#install.packages("knitr")
Warning message:
Unknown or uninitialised column: `Country`.
#install.packages("grid")
#install.packages("MLRMPA")
#install.packages("dprep")
#install.packages("normalr")
#install.packages("ggcorrplot")
#install.packages("RColorBrewer")
#install.packages("rgdal")
#install.packages("jsonlite")
#install.packages("RColorBrewer")
#install.packages("readr")
Sys.setlocale("LC_ALL", "English")
[1] "LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252"
Sys.setenv("LANGUAGE"="En")
library(gridExtra)
There were 12 warnings (use warnings() to see them)
library(dplyr)
library(lubridate)
library(magrittr)
library(ggplot2)
library(tidyr)
library(knitr)
library(normalr)
library(ggcorrplot)
library(leaflet)
library(plotly)
library(RColorBrewer)
library(readr)
#??src_mysql
my_db <- src_mysql(
dbname = "covid19",
host = "cvd-database.cmxpnyylhkiw.us-east-2.rds.amazonaws.com",
user = "admin",
password = "RprojectCovid19"
)
my_db
src: mysql 8.0.21 [admin@cvd-database.cmxpnyylhkiw.us-east-2.rds.amazonaws.com:/covid19]
tbls: covid19_confirmed, covid19_confirmed_global, covid19_deaths, covid19_deaths_global, covid19_recovered,
covid19_recovered_global, covid19_thailand, covidus, data_ethnic, data_gender, data_health, data_lockdown,
data_population, data_testing, gdp19, healthranking, population, pornhub, sars_2003, us_gdp, us_homeless,
us_latlong, world_temp
##import data
df_conf <- tbl(my_db, sql("select * from covid19_confirmed"))
df_conf <- as.data.frame(df_conf)
df_conf
df_deaths <- tbl(my_db, sql("select * from covid19_deaths"))
df_deaths <- as.data.frame(df_deaths)
df_deaths
df_recover <- tbl(my_db, sql("select * from covid19_recovered"))
df_recover <- as.data.frame(df_recover)
df_recover
##check the time frame of the data
n.col <- ncol(df_conf)
dates <- names(df_conf)[5:n.col]%>% mdy()
range(dates)
[1] "2020-01-22" "2021-01-22"
min.date <- min(dates)
max.date <- max(dates)
min.date.txt <- min.date %>% format('%d %b %Y')
max.date.txt <- max.date %>% format('%d %b %Y')
#clean data
cleanData <- function(data) {
## remove some columns
data %<>% select(-c(Province.State, Lat, Long)) %>% rename(country=Country.Region)
## convert from wide to long format
data %<>% gather(key=date, value=count, -country)
## convert from character to date
data %<>% mutate(date = date %>% mdy())
## aggregate by country
data %<>% group_by(country, date) %>% summarise(count=sum(count, na.rm=T)) %>% as.data.frame()
return(data)
}
## clean the three data sets
data.confirmed <- df_conf %>% cleanData() %>% rename(confirmed=count)
`summarise()` regrouping output by 'country' (override with `.groups` argument)
data.deaths <- df_deaths %>% cleanData() %>% rename(deaths=count)
`summarise()` regrouping output by 'country' (override with `.groups` argument)
data.recovered <- df_recover %>% cleanData() %>% rename(recovered=count)
`summarise()` regrouping output by 'country' (override with `.groups` argument)
data <- data.confirmed %>% merge(data.deaths, all=T) %>% merge(data.recovered, all=T)
data
## countries/regions with confirmed cases, excl. cruise ships
countries <- data %>% pull(country) %>% setdiff('Cruise Ship')
data
data.world <- data %>% group_by(date) %>%
summarise(country='World',
confirmed = sum(confirmed, na.rm=T),
deaths = sum(deaths, na.rm=T),
recovered = sum(recovered, na.rm=T))
`summarise()` ungrouping output (override with `.groups` argument)
data %<>% rbind(data.world)
data
data %<>% mutate(current.confirmed = confirmed - deaths - recovered)
#View(data)
# World confirmed cases map
world.location <- df_conf
#View(world.location )
world.location$confirmed <- world.location [, ncol(world.location )]
world.location %<>% select(c(Country.Region, Province.State, Lat, Long, confirmed)) %>%
mutate(txt=paste0(Country.Region, ' - ', Province.State, ': ', confirmed))
map.confirmed <- leaflet(width=1200, height=800) %>% addTiles()
# circle marker (units in pixels)
map.confirmed %<>% addCircleMarkers(world.location$Long, world.location$Lat,
# radius=2+log2(x$confirmed),
radius=0.02*sqrt(world.location$confirmed),
stroke=F,
color='red', fillOpacity=0.3,
label=world.location$txt)
# world
map.confirmed
## China
map.confirmed %>% setView(95, 35, zoom=4)
## Australia and New Zealand
map.confirmed %>% setView(135, -27, zoom=4)
## US and Canada
map.confirmed %>% setView(-105, 40, zoom=4)
## Europe
map.confirmed %>% setView(10, 50, zoom=4)
# new.confirmed new.deaths new.recovered
data %<>% arrange(country, date)
n <- nrow(data)
day1 <- min(data$date)
data %<>% mutate(new.confirmed = ifelse(date == day1, NA, confirmed - lag(confirmed, n=1)),
new.deaths = ifelse(date == day1, NA, deaths - lag(deaths, n=1)),
new.recovered = ifelse(date == day1, NA, recovered - lag(recovered, n=1)))
data %<>% mutate(new.confirmed = ifelse(new.confirmed < 0, 0, new.confirmed),
new.deaths = ifelse(new.deaths < 0, 0, new.deaths),
new.recovered = ifelse(new.recovered < 0, 0, new.recovered))
## death rate based on total deaths and recovered cases
data %<>% mutate(rate.upper = (100 * deaths / (deaths + recovered)) %>% round(1),
rate.upper = ifelse(is.nan(rate.upper), 0, rate.upper))
## lower bound: death rate based on total confirmed cases
data %<>% mutate(rate.lower = (100 * deaths / confirmed) %>% round(1),
rate.lower = ifelse(is.nan(rate.lower), 0, rate.lower))
## death rate based on the number of death/recovered on every single day
data %<>% mutate(rate.daily = (100 * new.deaths / (new.deaths + new.recovered)) %>% round(1),
rate.daily = ifelse(is.nan(rate.daily), 0, rate.daily))
#View(data)
## convert from wide to long format
data.long <- data %>%
select(c(country, date, confirmed, current.confirmed, recovered, deaths)) %>%
gather(key=type, value=count, -c(country, date))
## set factor levels to show them in a desirable order
data.long %<>% mutate(type=recode_factor(type, confirmed='Total Confirmed',
current.confirmed='Current Confirmed',
recovered='Recovered',
deaths='Deaths'))
#View(data.long)
##Number of case World
There were 35 warnings (use warnings() to see them)
world <- filter(data.long,country == 'World')
world
plot_world.area <- world %>% filter(type != 'Total Confirmed') %>%
ggplot(aes(x=date, y=count)) +
geom_area(aes(fill=type), alpha=0.5) +
labs(title=paste0('Numbers of Cases Worldwide - ', max.date.txt)) +
scale_fill_manual(values=c('red', 'green', 'black')) +
theme(legend.title=element_blank(), legend.position='bottom',
plot.title = element_text(size=7),
axis.title.x=element_blank(),
axis.title.y=element_blank(),
legend.key.size=unit(0.2, 'cm'),
legend.text=element_text(size=6),
axis.text=element_text(size=7),
axis.text.x=element_text(angle=45, hjust=1))
plot_world.line <- world %>%
ggplot(aes(x=date, y=count,color = type)) +
geom_line(aes(color=type)) +
labs(title=paste0('Numbers of Cases Worldwide (log scale) - ', max.date.txt)) +
#scale_color_manual(values=c('blue', 'red', 'green', 'black')) +
theme(legend.title=element_blank(), legend.position='bottom',
plot.title = element_text(size=7),
axis.title.x=element_blank(),
axis.title.y=element_blank(),
legend.key.size=unit(0.2, 'cm'),
legend.text=element_text(size=6),
axis.text=element_text(size=7),
axis.text.x=element_text(angle=45, hjust=1)) +
scale_y_continuous(trans='log10')
## show two plots side by side
grid.arrange(plot_world.area, plot_world.line, ncol=2)

gly.plot_world.area <- ggplotly(plot_world.area)
gly.plot_world.area
gly.plot_world.line <- ggplotly(plot_world.line)
gly.plot_world.line
data.world.all <- data %>% filter(country == "World")
## a scatter plot with a smoothed line and vertical x-axis labels
plot_acc.deaths <- ggplot(data.world.all, aes(x=date, y=deaths)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='Accumulative Deaths') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_recovered <- ggplot(data.world.all, aes(x=date, y=recovered)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='Accumulative Recovered Cases') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_new.deaths <- ggplot(data.world.all, aes(x=date, y=new.deaths)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='New Deaths') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_new.recovered <- ggplot(data.world.all, aes(x=date, y=new.recovered)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='New Recovered Cases') +
theme(axis.text.x=element_text(angle=45, hjust=1))
## show four plots together, with 2 plots in each row
grid.arrange(plot_acc.deaths, plot_recovered, plot_new.deaths, plot_new.recovered, nrow=2)

## convert from wide to long format, for drawing area plots
rates.long <- data %>%
select(c(country, date, rate.upper, rate.lower, rate.daily)) %>%
gather(key=type, value=count, -c(country, date))
# set factor levels to show them in a desirable order
rates.long %<>% mutate(type=recode_factor(type, rate.daily='Daily',
rate.upper='Upper bound',
rate.lower = 'Lower bound'))
#View(rates.long)
plot_rates <- rates.long %>% filter(country == "World") %>%
ggplot(aes(x = date, y = count, color = type)) +
geom_line() +
labs(title = "World Death Rate (%)") +
xlab("") + ylab("Death Rate (%)")
plot_rates

gly.rates <- ggplotly(plot_rates)
gly.rates
## ranking by confirmed cases
data.latest.all <- data %>% filter(date == max(date)) %>%
select(country, date,confirmed, new.confirmed, current.confirmed,
recovered, deaths, new.deaths, death.rate=rate.lower) %>%
mutate(ranking = dense_rank(desc(confirmed))) %>%
arrange(ranking)
#View(data.latest.all)
k <- 20
## top 20 countries: 21 incl. 'World'
top.countries <- data.latest.all %>% filter(ranking <= k + 1) %>%
arrange(ranking) %>% pull(country) %>% as.character()
top.countries %>% setdiff('World') %>% print()
[1] "US" "India" "Brazil" "Russia" "United Kingdom" "France"
[7] "Spain" "Italy" "Turkey" "Germany" "Colombia" "Argentina"
[13] "Mexico" "Poland" "South Africa" "Iran" "Ukraine" "Peru"
[19] "Indonesia" "Netherlands"
df <- data.long %>% filter(country %in% top.countries) %>%
mutate(country=country %>% factor(levels=c(top.countries))) %>%
filter(country != 'World' & type != 'Total Confirmed') %>%
ggplot(aes(x=date, y=count, color=type)) +
geom_line()+
#geom_area(alpha=0.5) +
# xlab('') + ylab('') +
labs(title=paste0('Numbers of COVID-19 Cases in Top 20 Countries - ',
max.date.txt)) +
scale_color_manual(values=c('red', 'green', 'black')) +
theme(legend.title=element_blank(), legend.position='bottom',
plot.title = element_text(size=12),
axis.title.x=element_blank(),
axis.title.y=element_blank(),
legend.key.size=unit(0.4, 'cm'),
legend.text=element_text(size=12),
strip.text.x=element_text(size=12),
axis.text=element_text(size=12),
axis.text.x=element_text(angle=45, hjust=1)) +
facet_wrap(~country, ncol=4, scales='free_y') +
scale_y_continuous(trans='log10')
df

data.latest <- data.latest.all %>% filter(!is.na(country)) %>%
mutate(country=ifelse(ranking <= k + 1, as.character(country), 'Others')) %>%
mutate(country=country %>% factor(levels=c(top.countries, 'Others')))
data.latest %<>% group_by(country) %>%
summarise(confirmed=sum(confirmed), new.confirmed=sum(new.confirmed),
current.confirmed=sum(current.confirmed),
recovered=sum(recovered), deaths=sum(deaths), new.deaths=sum(new.deaths)) %>%
mutate(death.rate=(100 * deaths/confirmed) %>% round(1))
`summarise()` ungrouping output (override with `.groups` argument)
data.latest
data.latest %<>% select(c(country, confirmed, deaths, death.rate,
new.confirmed, new.deaths, current.confirmed,recovered)) %>%
mutate(recover.rate=(100 * recovered/confirmed) %>% round(1))
data.latest
df_pop <- tbl(my_db, sql("select * from population "))
df_pop <- as.data.frame(df_pop)
df_pop <- rename(df_pop,"country"="Country")
# Add World Population
world_pop <- sum(df_pop$`Population (2020)`)
df_pop[nrow(df_pop) + 1,] = c("World", world_pop)
# Add Other Countries Population
top_pop <- filter(df_pop, df_pop$country %in% top.countries & df_pop$country != "World")
top_pop <- sum(top_pop$`Population (2020)` %>% as.numeric())
others_pop <- (world_pop - top_pop)
df_pop[nrow(df_pop) + 1,] = c("Others", others_pop)
#View(df_pop)
data.latest <- merge(x = data.latest, y = df_pop, by = "country", all.x = TRUE)
data.latest
data.latest <- rename(data.latest,"population" = "Population (2020)")
data.latest$population <- data.latest$population %>% as.numeric()
data.latest <- data.latest %>%
select(c(country, confirmed, deaths, death.rate,
new.confirmed, new.deaths,
current.confirmed, recovered, recover.rate, population)) %>%
mutate(confirm.rate = (100 * confirmed / population) %>% round(1))
data.latest
## % of death
data.latest %>% mutate(death.rate=death.rate %>% format(nsmall=1) %>% paste0('%'))
## convert from wide to long format, for drawing area plots
data.latest.long <- data.latest %>% filter(country!='World') %>%
gather(key=type, value=count, -country)
## set factor levels to show them with proper text and in a desirable order
data.latest.long %<>% mutate(type=recode_factor(type,
confirmed='Total Confirmed',
deaths='Total Deaths',
death.rate='Death Rate (%)',
new.confirmed='New Confirmed (compared with one day before)',
new.deaths='New Deaths (compared with one day before)',
current.confirmed='Current Confirmed',
recover.rate = 'Recovered Rate(%)',
confirm.rate = 'Confirmed Rate(%)'))
#View(data.latest.long)
data.one.dem <- filter(data.latest.long,type=='Total Confirmed'
| type=='Total Deaths'
| type=='Current Confirmed')
data.two.dem <- filter(data.latest.long,type=='Death Rate (%)'
# | type=='New Confirmed (compared with one day before)'
# | type=='New Deaths (compared with one day before)'
| type=='Recovered Rate(%)'
| type=='Confirmed Rate(%)')
data.two.dem
## bar chart
data.one.dem %>% ggplot(aes(x=country, y=count, fill=country, group=country)) +
geom_bar(stat='identity') +
geom_text(aes(label=count, y=count), size=3, vjust=0) +
xlab('') + ylab('') +
labs(title=paste0('Top 20 Countries with Most Confirmed Cases - ', max.date.txt)) +
scale_fill_discrete(name='Country', labels=aes(count)) +
theme(legend.title=element_blank(),
legend.position='none',
plot.title=element_text(size=13),
axis.text=element_text(size=8),
axis.text.x=element_text(angle=45, hjust=1)) +
facet_wrap(~type, ncol=1, scales='free_y')

data.two.dem$facet <- factor(data.two.dem$type, levels = c('Confirmed Rate(%)', 'Recovered Rate(%)','Death Rate (%)'))
data.two.dem %>%
ggplot(aes(x=country, y=count, fill=country, group=country)) +
geom_bar(stat='identity') +
geom_text(aes(label=count, y=count), size=4, vjust=0) +
xlab('') + ylab('') +
labs(title=paste0('Top 20 Countries with Most Confirmed Cases - ', max.date.txt)) +
scale_fill_discrete(name='Country', labels=aes(count)) +
theme(legend.title=element_blank(),
legend.position='none',
plot.title=element_text(size=13),
axis.text=element_text(size=10),
axis.text.x=element_text(size=10,angle=45, hjust=1)) +
facet_wrap(~facet, ncol=1, scales='free_y')

##GDP
df_gdp2019 <- tbl(my_db, sql("select * from gdp19"))
df_gdp2019 <- as.data.frame(df_gdp2019)
df_gdp2019
#healthranking data
df_healt <- tbl(my_db, sql("select * from healthranking"))
df_healt <- as.data.frame(df_healt)
df_healt <- select(df_healt,c("country","healthCareIndex"))
#View(df_healt)
#Top20Pornhub data
#df_pornhub <- tbl(my_db, sql("select * from Pornhub"))
#df_pornhub <- as.data.frame(df_pornhub)
#df_pornhub
# temp data
df_temp <- tbl(my_db, sql("select * from world_temp"))
df_temp <- as.data.frame(df_temp)
df_temp$Country[df_temp$Country == "United States"] <- "US"
df_city <- select(df_temp,c("Country","City")) %>%
rename(country=Country) %>%
rename(city=City)
numofcity <- aggregate(city ~ country, data = df_city, length)
#clean data
df_temp <- select(df_temp,c("Country","Avg_Year")) %>%
rename(country=Country)
#View(df_temp)
#df_temp <- data.frame(country=df_temp[,1],avg=rowMeans(df_temp[,-1]))
df_temp <- df_temp %<>% group_by(country) %>% summarise(avg_temp = mean(Avg_Year,na.rm = TRUE)%>% round(1))
`summarise()` ungrouping output (override with `.groups` argument)
df_temp
##temp bar
df_temp.all <- df_temp %>% merge(data.latest.all)
#View(df_temp.all)
df_temp_top.all <- df_temp.all %>% filter(country %in% top.countries) %>%
mutate(ranking = ranking - 1) %>%
arrange(ranking)
#View(df_temp_top.all)
g_temp_top <- df_temp_top.all %>%
ggplot(aes(x = reorder(country, ranking), y = avg_temp, fill = avg_temp)) +
labs(title=paste0("Temperature in Top 20 countries"), subtitle = "Average Temperature in Top 20 countries with most confirmed cases (°C) (2020)") +
scale_color_gradient(low = "#93DBFF", high = "#FF7771") +
geom_text(aes(label=avg_temp, y=avg_temp), size=4, vjust=-0.5) +
geom_bar(stat = "identity", position = "dodge") +
theme(
legend.title=element_blank(),
legend.position='none',
plot.title=element_text(size = 15, hjust = 0.5),
plot.subtitle = element_text(size = 12, hjust = 0.5),
axis.text=element_text(size=8),
axis.text.x=element_text(size = 9, angle=45, hjust=1)) +
scale_x_discrete(name = "Country") +
scale_y_discrete(name = "Average Temperature")
#labs(title = "Temperature in Top 20 countries", subtitle = "Temperature in Top 20 countries with most confirmed cases (°C)")
#g_temp_top
g_temp_top

#data.latest.all
lat.long <- rename(df_conf, "country" = "Country.Region", "city" = "Province.State") %>%
select("country", "Lat", "Long") %>%
merge(df_temp.all[c("country","confirmed", "recovered", "deaths", "avg_temp", "ranking")], by = "country") %>%
distinct(country, .keep_all = TRUE) %>%
mutate(ranking = ranking - 1) %>%
arrange(ranking)
#View(lat.long)
label_world <- lat.long
label_world$avg_temp <- as.numeric(label_world[, names(label_world) %in% c("avg_temp")])
label_world <- label_world %>%
mutate(txt=paste0('<b>',ranking, '</b>',
'<br/>','<b>',country, '</b>',
'<br/>', "Temperature: ",avg_temp, ' °C',
'<br/>', "Confirmed: ", confirmed,
'<br/>', "Deaths: ", deaths,
'<br/>', "Recovered: ", recovered
))
label_world$txt <- label_world$txt %>% lapply(htmltools::HTML)
label_world
label_top <- label_world %>% filter(ranking < 21)
label_top
# Temperature Map
wpal <- colorNumeric("YlOrRd", label_world$avg_temp, n = 4)
topIcon <- makeIcon("star.png",
#iconUrl = "https://static.vecteezy.com/system/resources/previews/001/189/063/non_2x/star-rounded-png.png",
iconWidth = 10, iconHeight = 10
#iconAnchorX = 20, iconAnchorY = 20
)
label_world <- label_world %>% filter(ranking > 20)
m <- leaflet(width=1200, height=800) %>% addTiles()
m %<>% addCircleMarkers(label_world$Long, label_world$Lat,
# radius=2+log2(x$confirmed),
radius=10,#*log2(m.world$avg.temp),
stroke=F,
#color='red',
color = wpal(label_world$avg_temp),
fillOpacity=0.5,
#popup=label.top$txt
label= label_world$txt,
group = "World"
) %>%
addCircleMarkers(label_top$Long, label_top$Lat,
# radius=2+log2(x$confirmed),
radius=10,#*log2(m.world$avg.temp),
stroke=F,
#color='red',
color = wpal(label_top$avg_temp),
fillOpacity=0.5,
#popup=label.top$txt
label= label_top$txt,
group = "Top 20 Countries"
) %>%
addLabelOnlyMarkers(label_top$Long, label_top$Lat, label = label_top$ranking,
labelOptions = labelOptions(noHide = TRUE, textOnly = TRUE,
direction = "head",
offset = c(5,4)),
group = "Top 20 Countries") %>%
addLegend("bottomright", pal = wpal, values = label_world$avg_temp, opacity = 1,
labFormat = labelFormat(suffix = " °C"),
title = "Temperature") %>%
addLayersControl(
#baseGroups = c("OSM (default)", "Toner", "Toner Lite"),
overlayGroups = c("Top 20 Countries", "World"),
options = layersControlOptions(collapsed = FALSE)
)
m
gdp.top20 <- df_gdp2019 %>%
select(c("rank", "country", "GDP (millions of US dollars)")) %>%
merge(data.latest.all %>%
select(country, ranking, confirmed, recovered, deaths) %>%
filter(country %in% top.countries & country != "World"), by = "country") %>%
arrange(ranking) %>%
mutate(ranking = ranking - 1)
gdp.top20 %<>% rename("GDP" = "GDP (millions of US dollars)")
gdp.top20
g <- ggplot(gdp.top20, aes(x = GDP, y = reorder(country, -ranking))) +
geom_histogram(stat = "identity", aes(fill = GDP))+
scale_fill_gradient("GDP", low = "#FF4038", high = "#50E952") +
labs(title=paste0("GDP of Top 20 Countries in 2019 (millions of US dollars)")) +
geom_text(aes(label=GDP, x = GDP), size=3.5, hjust=-0.2) +
xlab("GDP (millions of US dollars)") +
ylab("") +
theme(legend.title=element_blank())
Ignoring unknown parameters: binwidth, bins, pad
g

gly.top.gdp <- ggplotly(g)
gly.top.gdp
# Pornhub
df_pornhub <- tbl(my_db, sql("select * from pornhub"))
df_pornhub <- as.data.frame(df_pornhub)
df_pornhub
# Sars
df_sars <- tbl(my_db, sql("select * from sars_2003"))
df_sars <- as.data.frame(df_sars)
#View(df_sars)
dates.s <- df_sars[,1]%>% mdy()
range(dates.s)
[1] "2003-03-17" "2003-07-11"
min.date.s <- min(dates.s)
max.date.s <- max(dates.s)
min.date.txt.s <- min.date.s %>% format('%d %b %Y')
max.date.txt.s <- max.date.s %>% format('%d %b %Y')
day1.sars <- min.date.s
# clean data Sars
data.sars <- df_sars %>% rename(c("date" = "Date", "confirmed" = "Cumulative_number" ,"deaths" = "Number_deaths", "recovered" = "Number_recovered")) %>%
mutate(date = date %>% mdy()) %>%
group_by(country, date) %>% as.data.frame()
# Add World's Sars cases
world.sars <- data.sars %>% group_by(date) %>%
summarise(country='World',
confirmed = sum(confirmed, na.rm=T),
deaths = sum(deaths, na.rm=T),
recovered = sum(recovered, na.rm=T))
`summarise()` ungrouping output (override with `.groups` argument)
data.sars %<>% rbind(world.sars)
data.sars %<>% mutate(current.confirmed = confirmed - deaths - recovered)
#View(world.sars)
data.sars
#rate
data.sars %<>% arrange(country, date)
n <- nrow(data.sars)
day1.sars <- min(data.sars$date)
data.sars %<>% mutate(new.confirmed = ifelse(date == day1.sars, NA, confirmed - lag(confirmed, n=1)),
new.deaths = ifelse(date == day1.sars, NA, deaths - lag(deaths, n=1)),
new.recovered = ifelse(date == day1.sars, NA, recovered - lag(recovered, n=1)))
data.sars %<>% mutate(new.confirmed = ifelse(new.confirmed < 0, 0, new.confirmed),
new.deaths = ifelse(new.deaths < 0, 0, new.deaths),
new.recovered = ifelse(new.recovered < 0, 0, new.recovered))
## death rate based on total deaths and recovered cases
data.sars %<>% mutate(rate.upper = (100 * deaths / (deaths + recovered)) %>% round(1),
rate.upper = ifelse(is.nan(rate.upper), 0, rate.upper))
## lower bound: death rate based on total confirmed cases
data.sars %<>% mutate(rate.lower = (100 * deaths / confirmed) %>% round(1),
rate.lower = ifelse(is.nan(rate.lower), 0, rate.lower))
## death rate based on the number of death/recovered on every single day
data.sars %<>% mutate(rate.daily = (100 * new.deaths / (new.deaths + new.recovered)) %>% round(1),
rate.daily = ifelse(is.nan(rate.daily), 0, rate.daily))
#View(data.sars)
## convert from wide to long format
data.sars.long <- data.sars %>%
select(c(country, date, confirmed, current.confirmed, recovered, deaths)) %>%
gather(key=type, value=count, -c(country, date))
## set factor levels to show them in a desirable order
data.sars.long %<>% mutate(type=recode_factor(type, confirmed='Total Confirmed',
current.confirmed='Current Confirmed',
recovered='Recovered',
deaths='Deaths'))
#View(data.sars.long)
# World sars' long data
world.sars.long <- data.sars.long %>%
filter(country == "World")
#View(world.sars.long)
g <- ggplot(world.sars.long, aes(date, count, color = type)) +
geom_line()+
labs(title = "Number of Cases Worldwide: SARs")+
xlab("")+
ylab("")
g

gly.g <- ggplotly(g)
gly.g
gly.plot_world.line
df_sars_lastdate <- data.sars %>%
filter(date == max.date.s)
df_sars_lastdate
## Current Confirmed Cases
data.sars.world <- data.sars %>% filter(country=='World')
#View(data.sars.world)
n <- nrow(data.sars.world)
#View(data.sars.world)
plot_sars.currconf <- ggplot(data.sars.world, aes(x=date, y=current.confirmed)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='Current Confirmed Cases') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_sars.newconf <- ggplot(data.sars.world, aes(x=date, y=new.confirmed)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='Daily New Confirmed Cases') +
theme(axis.text.x=element_text(angle=45, hjust=1))
## show two plots side by side
grid.arrange(plot_sars.currconf, plot_sars.newconf, ncol=2)

## a scatter plot with a smoothed line and vertical x-axis labels
plot_sars.accdeaths <- ggplot(data.sars.world, aes(x=date, y=deaths)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='Accumulative Deaths') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_sars.accrecov <- ggplot(data.sars.world, aes(x=date, y=recovered)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='Accumulative Recovered Cases') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_sars.newdeaths <- ggplot(data.sars.world, aes(x=date, y=new.deaths)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='New Deaths') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot_sars.newrecov <- ggplot(data.sars.world, aes(x=date, y=new.recovered)) +
geom_point() + geom_smooth() +
xlab('') + ylab('Count') + labs(title='New Recovered Cases') +
theme(axis.text.x=element_text(angle=45, hjust=1))
## show four plots together, with 2 plots in each row
grid.arrange(plot_sars.accdeaths, plot_sars.accrecov, plot_sars.newdeaths, plot_sars.newrecov, nrow=2)

heatmap.2(as.matrix(norm_data_plot),
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
scale="none",
col = colorRampPalette(c("#6D9EC1","white","#E46726"))(n = 200),
margins=c(10,6),trace="column")

#correlation
corr_data %<>% select(c(GDP,confirm.rate,death.rate,recover.rate,healthcare,avg_temp,Pornhub, sars))
head(corr_data)
cor(corr_data)
GDP confirm.rate death.rate recover.rate healthcare avg_temp Pornhub
GDP 1.00000000 0.4342894 -0.23113988 -0.3937431 0.25183938 -0.05173909 0.71560225
confirm.rate 0.43428937 1.0000000 -0.34422360 -0.7182585 0.48153104 -0.36431596 0.43162144
death.rate -0.23113988 -0.3442236 1.00000000 0.1733012 -0.17258588 0.34643395 -0.04353792
recover.rate -0.39374308 -0.7182585 0.17330123 1.0000000 -0.76240251 0.26074521 -0.40887125
healthcare 0.25183938 0.4815310 -0.17258588 -0.7624025 1.00000000 -0.06716339 0.37579971
avg_temp -0.05173909 -0.3643160 0.34643395 0.2607452 -0.06716339 1.00000000 -0.11100774
Pornhub 0.71560225 0.4316214 -0.04353792 -0.4088713 0.37579971 -0.11100774 1.00000000
sars 0.97481805 0.4506609 -0.21325881 -0.3582828 0.16579089 -0.05509975 0.68006901
sars
GDP 0.97481805
confirm.rate 0.45066089
death.rate -0.21325881
recover.rate -0.35828282
healthcare 0.16579089
avg_temp -0.05509975
Pornhub 0.68006901
sars 1.00000000
plotcorr <- ggcorrplot(cor(corr_data),hc.order = TRUE,
outline.color = "white",
colors = c("#6D9EC1","white","#E46726"),
lab = TRUE)
plotcorr

## Data in US
df_us <- tbl(my_db, sql("select * from covidus"))
df_us <- as.data.frame(df_us)
df_us <- select(df_us, date, state, cases, deaths)
df_us_pop <- tbl(my_db, sql("select * from data_population"))
df_us_pop <- as.data.frame(df_us_pop)
df_us_gender <- tbl(my_db, sql("select * from data_gender"))
df_us_gender <- as.data.frame(df_us_gender)
df_us_ethnic <- tbl(my_db, sql("select * from data_ethnic"))
df_us_ethnic <- as.data.frame(df_us_ethnic)
df_us_lockdown <- tbl(my_db, sql("select * from data_lockdown"))
df_us_lockdown <- as.data.frame(df_us_lockdown)
df_us_health <- tbl(my_db, sql("select * from data_health"))
df_us_health <- as.data.frame(df_us_health)
df_us_testing <- tbl(my_db, sql("select * from data_testing"))
df_us_testing <- as.data.frame(df_us_testing)
df_us <- df_us %>%
mutate(date = date %>% mdy()) %>%
rename("confirmed" = "cases")
#data.us
dates.us <- df_us[,1]
range(dates.us)
[1] "2020-01-21" "2020-12-24"
min.date.us <- min(dates.us)
max.date.us <- max(dates.us)
min.date.txt.us <- min.date.us %>% format('%d %b %Y')
max.date.txt.us <- max.date.us %>% format('%d %b %Y')
day1.us <- min(df_us$date)
data.us.total <- df_us %>% group_by(date) %>%
summarise(state='US',
confirmed = sum(confirmed, na.rm=T),
deaths = sum(deaths, na.rm=T))
`summarise()` ungrouping output (override with `.groups` argument)
#View(data.us.total)
data.us <- df_us
data.us %<>% rbind(data.us.total)
View(data.us)
data.us.long <- data.us %>%
gather(key = type, value = count, -c(date, state))
View(data.us.long)
us.total <- data.us.total %>%
mutate(new.confirmed = ifelse(date == day1, 0, confirmed - lag(confirmed, n=1)),
new.deaths = ifelse(date == day1, 0, deaths - lag(deaths, n=1)))
us.total %<>% mutate(new.confirmed = ifelse(new.confirmed < 0, 0, new.confirmed),
new.deaths = ifelse(new.deaths < 0, 0, new.deaths))
us.total
us.total.long <- us.total %>%
gather(key = type, value = count, -c(date, state))
#View(us.total.long)
plot_us.cases <- data.us.long %>% filter(state == "US") %>%
ggplot(aes(x = date, y = count)) +
geom_area(aes(fill=type), alpha=0.5) +
# labs(title = paste0("Cumulative cases in US : ", min.date.txt.us, '-', max.date.txt.us, " (Stack, Log Scale)")) +
scale_y_continuous(trans='log10') +
scale_fill_manual(values=c('red', 'black'))+
ylab("") +
theme(axis.title = element_text(size=8))
#View(us.total.long )
plot_us.newconf <- us.total.long %>%
filter(type %in% c("new.confirmed","new.deaths")) %>%
ggplot(aes(x = date, y = count, color = type)) +
geom_line() +
#labs(title = paste0("Daily confirmed cases in US : ", min.date.txt.us, '-', max.date.txt.us)) +
xlab("Date") +
ylab("Confirmed cases")
gly.plot_us.cases <- ggplotly(plot_us.cases)
Transformation introduced infinite values in continuous y-axis
gly.plot_us.cases <- gly.plot_us.cases %>%
layout(title = paste0("<b>Cumulative cases in US : ", min.date.txt.us, '-', max.date.txt.us, " (Stack, Log Scale)"),
font=list(size = 10))
gly.plot_us.newconf <- ggplotly(plot_us.newconf) %>%
layout(title = paste0("<b>Daily confirmed cases in US : ", min.date.txt.us, '-', max.date.txt.us),
font=list(size = 10))
gly.plot_us.cases
gly.plot_us.newconf
NA
data.us %<>% mutate(new.confirmed = ifelse(date == day1, NA, confirmed - lag(confirmed, n=1)),
new.deaths = ifelse(date == day1, NA, deaths - lag(deaths, n=1)))
data.us %<>% mutate(new.confirmed = ifelse(new.confirmed < 0, 0, new.confirmed),
new.deaths = ifelse(new.deaths < 0, 0, new.deaths))
View(data.us)
#data.us.daily <- data.us %>% filter(state == "US")
#View(data.us.daily)
data.us.pop <- df_us_pop %>% select(State, Population) %>%
rename("state" = "State")
#data.us.pop
data.us.gender <- df_us_gender %>% select(State, Male, Female) %>%
rename("state" = "State")
data.us.gender
data.us.latest <- data.us %>%
filter(date == max.date.us) %>%
merge(data.us.pop, by = "state", all.x = T) %>%
merge(data.us.gender, by = "state", all.x = T)
#View(data.us.latest)
data.us.latest$Population[data.us.latest$state == "US"] <- sum(data.us.pop$Population)
data.us.latest %<>% mutate(ranking = dense_rank(desc(confirmed)),
confirmed.rate = (100 * confirmed / Population) %>% round(2),
death.rate = (100 * deaths / confirmed) %>% round(2),
Male.confirmed = (Male * confirmed) %>% round(0),
Female.confirmed = Female * confirmed %>% round(0),
Male.deaths = Male * deaths,
Female.deaths = Female * deaths) %>%
arrange(ranking)
top.us <- data.us.latest[,1]
top.us
[1] "US" "California" "Texas"
[4] "Florida" "Illinois" "New York"
[7] "Ohio" "Georgia" "Pennsylvania"
[10] "Tennessee" "Michigan" "Wisconsin"
[13] "North Carolina" "Indiana" "Arizona"
[16] "New Jersey" "Minnesota" "Missouri"
[19] "Massachusetts" "Alabama" "Virginia"
[22] "Colorado" "Louisiana" "South Carolina"
[25] "Iowa" "Oklahoma" "Maryland"
[28] "Utah" "Kentucky" "Washington"
[31] "Kansas" "Nevada" "Arkansas"
[34] "Mississippi" "Connecticut" "Nebraska"
[37] "Idaho" "New Mexico" "Oregon"
[40] "Puerto Rico" "South Dakota" "North Dakota"
[43] "Rhode Island" "Montana" "West Virginia"
[46] "Delaware" "Alaska" "Wyoming"
[49] "New Hampshire" "District of Columbia" "Maine"
[52] "Hawaii" "Guam" "Vermont"
[55] "Virgin Islands" "Northern Mariana Islands"
#View(top.us)
data.us.latest.show <- data.us.latest %>% select(state,date,confirmed,deaths, confirmed.rate, death.rate)
#View(data.us.latest.show)
# List of top 20 state
k <- 20
data.us.top <- data.us.latest %>%
filter(ranking <= k+1) %>%
arrange(ranking)
#View(data.us.top)
us.state.top <- data.us.top %>% pull(state) %>% as.character()
us.state.top %>% setdiff('US') %>% print()
[1] "California" "Texas" "Florida" "Illinois" "New York" "Ohio"
[7] "Georgia" "Pennsylvania" "Tennessee" "Michigan" "Wisconsin" "North Carolina"
[13] "Indiana" "Arizona" "New Jersey" "Minnesota" "Missouri" "Massachusetts"
[19] "Alabama" "Virginia"
# confirmed rate & death rate of top 20 state
g.rate <- data.us.latest %>% filter(state %in% us.state.top & state != "US") %>%
select(state, confirmed.rate, death.rate, ranking) %>%
gather(key = Type, value = Percent, -c(state, ranking)) %>%
ggplot(aes(x=reorder(state, -desc(ranking)), y=Percent, fill = Percent)) +
geom_bar(stat='identity') +
scale_fill_gradient(low = "#ebbc62", high = "#b42006") +
geom_text(aes(label=Percent, y=Percent), size=3, vjust=0) +
xlab('') + ylab('') +
labs(title=paste0('Confirmed Rate & Death Rate of Top 20 State in US')) +
#scale_fill_continuous(name='State', labels=aes(Percent)) +
theme(legend.title=element_blank(),
legend.position='none',
plot.title=element_text(size=13),
axis.text=element_text(size=8),
axis.text.x=element_text(angle=45, hjust=1)) +
facet_wrap(~Type, ncol=1, scales='free_y')
g.rate

us.confirmed.num <- data.us.top$confirmed[data.us.top$state == "US"] %>% as.numeric()
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
3: Unknown or uninitialised column: `Country`.
us.confirmed.num
[1] 18706036
data.us.top %<>% mutate(confirmed.per.us = (confirmed * 100 / us.confirmed.num) %>% round(1))
data.us.top
g.us.top1 <- data.us.top %>%
filter(state != "US") %>%
select(state, confirmed, confirmed.per.us, ranking) %>%
gather(key = Type, value = count, -c(state, ranking, confirmed.per.us)) %>%
arrange(ranking) %>%
ggplot(aes(x = reorder(state, ranking), y = count, fill = count)) +
geom_bar(stat = "identity") +
#labs(title = "20 state in US with most confirmed cases") +
scale_color_gradient(low = "#93DBFF", high = "#FF7771") +
xlab("") +
ylab("Confirmed Cases") +
geom_text(aes(label=paste0(confirmed.per.us, "%")), size=3, vjust=-0.5) +
theme(axis.text.x = element_text(angle = 90,vjust = 0.5))+
theme(
axis.line = element_blank(),
axis.ticks = element_blank(),
panel.grid.minor = element_blank(),
panel.grid.major = element_blank(),
panel.border = element_blank(),
panel.background = element_blank(),
#legend.position = "none"
)
g.us.top1

gly.us.top1 <- ggplotly(g.us.top1)
gly.us.top1
data.us.latest.corr <- df_us %>%filter(date==max.date.us)
#View(data.us.latest.corr)
data.us.latest.corr %<>% mutate(ranking = dense_rank(desc(confirmed))) %>% arrange(ranking)
top.us <- data.us.latest.corr[,2]
#top.us
# List of top 20 state
k <- 19
data.us.latest.top <- data.us.latest.corr %>%
filter(ranking <= k+1) %>%
arrange(ranking)
#View(data.us.latest.top)
us.state.top <- data.us.latest.top %>% pull(state) %>% as.character()
us.state.top %>% setdiff('US') %>% print()
[1] "California" "Texas" "Florida" "Illinois" "New York" "Ohio"
[7] "Georgia" "Pennsylvania" "Tennessee" "Michigan" "Wisconsin" "North Carolina"
[13] "Indiana" "Arizona" "New Jersey" "Minnesota" "Missouri" "Massachusetts"
[19] "Alabama" "Virginia"
#View(us.state.top)
#gender in us
df_us_gender <- tbl(my_db, sql("select * from data_gender"))
df_us_gender <- as.data.frame(df_us_gender)
df_us_gender <- select(df_us_gender,c("State","Male","Female"))
df_us_gender <- rename(df_us_gender,"state"="State")
df_us_gender
#population in us
df_us_pop <- tbl(my_db, sql("select * from data_population"))
df_us_pop <- as.data.frame(df_us_pop)
df_us_pop <- select(df_us_pop,c("State","Population"))
df_us_pop <- rename(df_us_pop,"state"="State")
df_us_pop
#lockdown in us
df_us_lockdown <- tbl(my_db, sql("select * from data_lockdown"))
df_us_lockdown <- as.data.frame(df_us_lockdown)
df_us_lockdown <- select(df_us_lockdown,c("State","Day lockdown"))
df_us_lockdown <- rename(df_us_lockdown,"state"="State")
df_us_lockdown
#GDP in us
df_us_gdp <- tbl(my_db, sql("select * from us_gdp"))
df_us_gdp <- as.data.frame(df_us_gdp )
df_us_gdp <- select(df_us_gdp ,c("State","GDPs"))
df_us_gdp <- rename(df_us_gdp ,"state"="State")
df_us_gdp
#homeless in us
df_us_homeless <- tbl(my_db, sql("select * from us_homeless"))
df_us_homeless <- as.data.frame(df_us_homeless)
df_us_homeless <- select(df_us_homeless,c("State","Homeless"))
df_us_homeless <- rename(df_us_homeless,"state"="State")
df_us_homeless
#View(df_us)
#merge
mergcountry = function(data1,data2){
data <- merge(x = data1, y = data2, by = "state", all.x = TRUE)
return(data)
}
data.top.state <- data.us.latest.top %>% filter(state != "US") %>%
select(state, confirmed, deaths)
#View(data.top.state)
df_Allus <- mergcountry(data.top.state, df_us_gender)
df_Allus <- mergcountry(df_Allus, df_us_pop)
df_Allus <- mergcountry(df_Allus, df_us_lockdown)
df_Allus <- mergcountry(df_Allus, df_us_gdp)
df_Allus <- mergcountry(df_Allus, df_us_homeless)
#View(df_Allus)
index <- is.na(df_Allus)
df_Allus[index] <- 0
#View(df_Allus)
normalize = function(data){
#return ((data - min(data,na.rm = TRUE))/(max(data,na.rm = TRUE) - min(data,na.rm = TRUE)))
z <- scale(data);
tanh(z/2)
}
Allus = as.data.frame(apply(df_Allus[,2:9],2,normalize))
corr_dataUS <- Allus
#View(corr_dataUS)
Allus$state <- c(df_Allus$state)
#View(Allus)
corr_dataUS <- rename(corr_dataUS ,"Daylockdown" = "Day lockdown")
heatmap.2(as.matrix(Allus_Heat),
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
3: Unknown or uninitialised column: `Country`.
4: Unknown or uninitialised column: `Country`.
5: Unknown or uninitialised column: `Country`.
6: Unknown or uninitialised column: `Country`.
scale="none",
col = colorRampPalette(c("#6D9EC1","white","#E46726"))(n = 200),
margins=c(10,6),trace="column")

#correlation
corr_dataUS %<>% select(c(confirmed,deaths,Male,Female,Population,Daylockdown,GDPs,Homeless))
head(corr_dataUS)
cor(corr_dataUS)
confirmed deaths Male Female Population Daylockdown GDPs
confirmed 1.000000000 0.792654639 0.105367061 -0.005945092 0.959941301 -0.2058259 0.91047827
deaths 0.792654639 1.000000000 0.008041947 0.017335533 0.852691131 -0.3797468 0.88497061
Male 0.105367061 0.008041947 1.000000000 -0.877329594 -0.009961687 -0.2898630 -0.02437981
Female -0.005945092 0.017335533 -0.877329594 1.000000000 0.069101784 0.3629328 0.02235864
Population 0.959941301 0.852691131 -0.009961687 0.069101784 1.000000000 -0.2155126 0.97052821
Daylockdown -0.205825865 -0.379746786 -0.289863040 0.362932815 -0.215512639 1.0000000 -0.36113782
GDPs 0.910478267 0.884970609 -0.024379811 0.022358644 0.970528214 -0.3611378 1.00000000
Homeless 0.791615898 0.796426563 -0.080155515 0.151488777 0.863178014 -0.3573560 0.89961070
Homeless
confirmed 0.79161590
deaths 0.79642656
Male -0.08015552
Female 0.15148878
Population 0.86317801
Daylockdown -0.35735597
GDPs 0.89961070
Homeless 1.00000000
ggcorrplot(cor(corr_dataUS),hc.order = TRUE,
outline.color = "white",
colors = c("#6D9EC1","white","#E46726"),
lab = TRUE)

NA
NA
NA
#Covid_Thailand
df_thai <- tbl(my_db, sql("select * from covid19_thailand"))
df_thai <- as.data.frame(df_thai)
#View(df_thai)
#clean Covid_Thailand
dates.th <- df_thai[,2]%>% mdy()
range(dates.th)
[1] "2020-01-12" "2021-01-15"
min.date.th <- min(dates.th)
max.date.th <- max(dates.th)
min.date.txt.th <- min.date.th %>% format('%d %b %Y')
max.date.txt.th <- max.date.th %>% format('%d %b %Y')
df_thai$announce_date <- mdy(df_thai$announce_date)
df_thai$notification_date <- mdy(df_thai$notification_date)
df_thai
df_thai <- df_thai %>% select(!No.) %>% select(!notification_date) %>%
group_by(announce_date)
df_thai
View(df_thai)
# Total confirmed cases in Thailand
There were 11 warnings (use warnings() to see them)
data.thai.count <- df_thai %>%
select(announce_date) %>%
summarise(confirmed = n()) %>% as.data.frame()
`summarise()` ungrouping output (override with `.groups` argument)
View(data.thai.count)
data.thai.count$cumulative_confirmed <- cumsum(data.thai.count[, 2])
data.thai.count <- data.thai.count %>% rename(new.confirmed=confirmed) %>% rename(confirmed=cumulative_confirmed)
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
thai.total.long <- data.thai.count %>%
gather(key = type, value = count, -c(announce_date))
thai.total.long
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
3: Unknown or uninitialised column: `Country`.
plot_thailand.conf <- thai.total.long %>%
filter(type %in% c("new.confirmed")) %>%
ggplot(aes(x = announce_date, y = count, color = type)) +
geom_line() +
xlab("Date") +
ylab("cases")
gly.plot_thai.cases <- ggplotly(plot_thailand.conf)
gly.plot_thai.cases <- gly.plot_thai.cases %>%
layout(title = paste0("<b>Daily Thai Confirmed Cases in Thailand : Jan 2020 - Jan 2021"),
font=list(size = 10))
gly.plot_thai.cases
plot_thailand.cumulative <- thai.total.long %>%
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
filter(type %in% c("confirmed")) %>%
ggplot(aes(x = announce_date, y = count)) +
geom_area(aes(fill=type), alpha=0.5) +
scale_y_continuous(trans='log10') +
scale_fill_manual(values=c('red'))+
ylab("") +
theme(axis.title = element_text(size=8))
gly.plot_thai.cumulative <- ggplotly(plot_thailand.cumulative)
gly.plot_thai.cumulative <- gly.plot_thai.cumulative %>%
layout(title = paste0("<b>Cumulative Cases in Thailand : Jan 2020 - Jan 2021 (Stack,Log scale)"),
font=list(size = 10))
gly.plot_thai.cumulative
## Thai Confirmed Cases (Jan 2020 - Jan 2021
plot1 <- ggplot(data.thai.count, aes(x=announce_date, y=cumulative_confirmed)) +
geom_point() + geom_smooth() +
xlab(" ") + ylab("Count") + labs(title='Thai Cumulative Confirmed Cases (Jan 2020 - Jan 2021)') +
theme(axis.text.x=element_text(angle=45, hjust=1))
plot2 <- ggplot(data.thai.count, aes(x=announce_date, y=comfirmed)) +
geom_point() + geom_smooth() +
xlab(" ") + ylab("Count")+ labs(title='Thai Confirmed Cases (Jan 2020 - Jan 2021 log scale)') +
theme(axis.text.x=element_text(angle=45, hjust=1))+scale_y_continuous(trans='log10')
## show two plots side by side
grid.arrange(plot1, plot2, ncol=1)
# Confirmed cases divided by sex (gender)
df_thai$sex[df_thai$sex == ""] <- "Unknown"
data.thai.gender <- df_thai %>%
group_by(sex) %>%
summarise(count = n()) %>%
mutate(percent = (count / sum(count) * 100) %>% round(2)) %>%
filter(percent>1)%>%
#mutate(pos = cumsum(percent) - 0.5*percent) %>%
arrange(desc(percent))
`summarise()` ungrouping output (override with `.groups` argument)
data.thai.gender
data.thai.gender$sex <- factor(data.thai.gender$sex, levels = as.character(data.thai.gender$sex))
data.thai.gender$sex
[1] Male Female
Levels: Male Female
g.th.gender <- data.thai.gender %>%
ggplot(aes(x = "", y = percent, fill = sex)) +
geom_bar(stat = "identity", width = 1) +
coord_polar("y") +
theme_void() +
labs(title='Gender of Thai Confirmed Cases (Jan 2020 - Jan 2021)')+
geom_text(aes(label = paste0(percent, "%")), color = "white", size = 5, position = position_stack(vjust = 0.5)) +
guides(fill = guide_legend(reverse = TRUE))
g.th.gender

# Confirmed cases divided by province_of_onset
data.thai.onset <- df_thai %>%
group_by(province_of_onset) %>%
summarise(count = n()) %>%
arrange(desc(count))%>%
rename(onset="count")%>%
rename(province="province_of_onset")
`summarise()` ungrouping output (override with `.groups` argument)
data.thai.onset$province[data.thai.onset$province == ""] <- "Unknown"
data.thai.onset
# Confirmed cases divided by province_of_isolation
data.thai.isolation <- df_thai %>%
group_by(province_of_isolation) %>%
summarise(count = n()) %>%
arrange(desc(count))%>%
rename(isolation ="count")%>%
rename(province="province_of_isolation")
`summarise()` ungrouping output (override with `.groups` argument)
data.thai.isolation$province[data.thai.isolation$province == ""] <- "Unknown"
data.thai.isolation
data.province <- merge(x = data.thai.onset, y = data.thai.isolation, by = "province", all.x = TRUE)
data.province <- data.province%>%filter(onset>125)
data.province
data.province.long <- data.province %>%
gather(key=type, value=count, -c(province))
data.province.long
# Confirmed cases divided by risk
data.thai.risk <- df_thai %>%
group_by(risk) %>%
summarise(count = n()) %>%
arrange(desc(count))
`summarise()` ungrouping output (override with `.groups` argument)
data.thai.risk
# Confirmed cases divided by risk
Warning messages:
1: Unknown or uninitialised column: `Country`.
2: Unknown or uninitialised column: `Country`.
3: Unknown or uninitialised column: `Country`.
data.thai.risk <- df_thai %>%
group_by(risk) %>%
summarise(count = n()) %>%
mutate(percent = (count / sum(count) * 100) %>% round(2)) %>%
filter(percent>4) %>%
arrange(desc(percent))
`summarise()` ungrouping output (override with `.groups` argument)
data.thai.risk
data.thai.risk$risk[data.thai.risk$risk == "C"] <- "Close contact with the patient"
data.thai.risk$risk[data.thai.risk$risk == "CGustDr Samut SaGGAF"] <- "Cluster Samut Sakhon"
data.thai.risk$risk[data.thai.risk$risk == "F"] <- "State Quarantine"
data.thai.risk$risk[data.thai.risk$risk == "G"] <- "Go to a crowded place"
data.thai.risk$risk[data.thai.risk$risk == "A"] <- "People travelling from abroad"
data.thai.risk$risk[data.thai.risk$risk == "UFGFAwF"] <- "Unknown"
data.thai.risk$risk[data.thai.risk$risk == "D"] <- "Career at risk"
data.thai.risk$risk[data.thai.risk$risk == "CGustDr RayAFg"] <- "Cluster Rayong"
data.thai.risk$risk[data.thai.risk$risk == "CGustDr Pattaya CasCFA"] <- "Cluster Pattaya Casino"
data.thai.risk$risk[data.thai.risk$risk == "CGustDr AFg TGAFg"] <- "Cluster Ang Thong"
data.thai.risk
data.thai.risk$risk <- factor(data.thai.risk$risk, levels = as.character(data.thai.risk$risk))
data.thai.risk$risk
[1] Close contact with the patient Cluster Samut Sakhon State Quarantine
[4] Go to a crowded place People travelling from abroad Unknown
[7] Career at risk Cluster Rayong
8 Levels: Close contact with the patient Cluster Samut Sakhon State Quarantine ... Cluster Rayong
g.th.risk <- data.thai.risk %>%
ggplot(aes(x = "", y = percent, fill = risk)) +
geom_bar(stat = "identity", width = 0.5) +
coord_polar("y") +
theme_void() +
labs(title='Risk of Thai Confirmed Cases(Jan 2020 - Jan 2021)')+
geom_text(aes(label = paste0(percent, "%")), color = "Black", size = 3, position = position_stack(vjust = 0.5)) +
guides(fill = guide_legend(reverse = TRUE))
g.th.risk

# Confirmed cases divided by age
data.thai.age <- df_thai %>%
group_by(age,sex) %>%
filter(age != 0)%>%
summarise(count = n()) %>%
arrange(desc(count))
`summarise()` regrouping output by 'age' (override with `.groups` argument)
data.thai.age
data.thai.age <- data.thai.age# %>% filter(sex != "Unknown")
ggplot(data.thai.age,aes(x=age,y=count,fill=sex))+geom_bar(stat = "identity")+
labs(title='Age of Thai Confirmed Cases (Stack)')+guides(fill=guide_legend(reverse = T))

# Confirmed cases divided by nationality
data.thai.nationality <- df_thai %>%
group_by(nationality) %>%
summarise(count = n()) %>%
filter(count > 11)%>%
arrange(desc(count))
`summarise()` ungrouping output (override with `.groups` argument)
data.thai.nationality$nationality[data.thai.nationality$nationality == "????????"] <- "Unknown"
data.thai.nationality$nationality[data.thai.nationality$nationality == ""] <- "Unknown"
data.thai.nationality$nationality[data.thai.nationality$nationality == "Burma"] <- "Myanmar"
data.thai.nationality
ggplot(data.thai.nationality)+
geom_bar(aes(x=nationality,y=count,fill = count),stat = "identity",show.legend = F)+
ggtitle('Nationality of Thai Confirmed Cases (log scale)')+
theme(plot.title = element_text(hjust = 0.5))+
theme(axis.text = element_text(size = 10),
axis.title = element_text(size = 10),
plot.title = element_text(size = 10))+
xlab("")+ylab("Number of Nationality")+
scale_y_continuous(trans='log10')+coord_flip()+scale_fill_gradient(low="blue",high = "red")

LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpgYGB7cn0NCiNpbnN0YWxsLnBhY2thZ2VzKCJrbml0ciIpDQojaW5zdGFsbC5wYWNrYWdlcygiZ3JpZCIpDQojaW5zdGFsbC5wYWNrYWdlcygiTUxSTVBBIikNCiNpbnN0YWxsLnBhY2thZ2VzKCJkcHJlcCIpDQojaW5zdGFsbC5wYWNrYWdlcygibm9ybWFsciIpDQojaW5zdGFsbC5wYWNrYWdlcygiZ2djb3JycGxvdCIpDQojaW5zdGFsbC5wYWNrYWdlcygiUkNvbG9yQnJld2VyIikNCiNpbnN0YWxsLnBhY2thZ2VzKCJyZ2RhbCIpDQojaW5zdGFsbC5wYWNrYWdlcygianNvbmxpdGUiKQ0KI2luc3RhbGwucGFja2FnZXMoIlJDb2xvckJyZXdlciIpDQojaW5zdGFsbC5wYWNrYWdlcygicmVhZHIiKQ0KU3lzLnNldGxvY2FsZSgiTENfQUxMIiwgIkVuZ2xpc2giKQ0KU3lzLnNldGVudigiTEFOR1VBR0UiPSJFbiIpDQoNCmBgYA0KDQpgYGB7cn0NCmxpYnJhcnkoZ3JpZEV4dHJhKQ0KbGlicmFyeShkcGx5cikNCmxpYnJhcnkobHVicmlkYXRlKQ0KbGlicmFyeShtYWdyaXR0cikNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkodGlkeXIpDQpsaWJyYXJ5KGtuaXRyKQ0KbGlicmFyeShub3JtYWxyKQ0KbGlicmFyeShnZ2NvcnJwbG90KQ0KbGlicmFyeShsZWFmbGV0KQ0KbGlicmFyeShwbG90bHkpDQpsaWJyYXJ5KFJDb2xvckJyZXdlcikNCmxpYnJhcnkocmVhZHIpDQoNCiM/P3NyY19teXNxbA0KbXlfZGIgPC0gc3JjX215c3FsKA0KICBkYm5hbWUgPSAiY292aWQxOSIsDQogIGhvc3QgPSAiY3ZkLWRhdGFiYXNlLmNteHBueXlsaGtpdy51cy1lYXN0LTIucmRzLmFtYXpvbmF3cy5jb20iLA0KICB1c2VyID0gImFkbWluIiwNCiAgcGFzc3dvcmQgPSAiUnByb2plY3RDb3ZpZDE5Ig0KKQ0KbXlfZGINCg0KIyNpbXBvcnQgZGF0YQ0KZGZfY29uZiA8LSB0YmwobXlfZGIsIHNxbCgic2VsZWN0ICogZnJvbSBjb3ZpZDE5X2NvbmZpcm1lZCIpKQ0KZGZfY29uZiA8LSBhcy5kYXRhLmZyYW1lKGRmX2NvbmYpDQpkZl9jb25mDQpkZl9kZWF0aHMgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gY292aWQxOV9kZWF0aHMiKSkNCmRmX2RlYXRocyA8LSBhcy5kYXRhLmZyYW1lKGRmX2RlYXRocykNCmRmX2RlYXRocw0KZGZfcmVjb3ZlciA8LSB0YmwobXlfZGIsIHNxbCgic2VsZWN0ICogZnJvbSBjb3ZpZDE5X3JlY292ZXJlZCIpKQ0KZGZfcmVjb3ZlciA8LSBhcy5kYXRhLmZyYW1lKGRmX3JlY292ZXIpDQpkZl9yZWNvdmVyDQpgYGANCg0KYGBge3J9DQojI2NoZWNrIHRoZSB0aW1lIGZyYW1lIG9mIHRoZSBkYXRhDQpuLmNvbCA8LSBuY29sKGRmX2NvbmYpDQpkYXRlcyA8LSBuYW1lcyhkZl9jb25mKVs1Om4uY29sXSU+JSBtZHkoKQ0KcmFuZ2UoZGF0ZXMpDQptaW4uZGF0ZSA8LSBtaW4oZGF0ZXMpDQptYXguZGF0ZSA8LSBtYXgoZGF0ZXMpDQptaW4uZGF0ZS50eHQgPC0gbWluLmRhdGUgJT4lIGZvcm1hdCgnJWQgJWIgJVknKQ0KbWF4LmRhdGUudHh0IDwtIG1heC5kYXRlICU+JSBmb3JtYXQoJyVkICViICVZJykNCmBgYA0KYGBge3J9DQojY2xlYW4gZGF0YQ0KY2xlYW5EYXRhIDwtIGZ1bmN0aW9uKGRhdGEpIHsNCiAgIyMgcmVtb3ZlIHNvbWUgY29sdW1ucw0KICBkYXRhICU8PiUgc2VsZWN0KC1jKFByb3ZpbmNlLlN0YXRlLCBMYXQsIExvbmcpKSAlPiUgcmVuYW1lKGNvdW50cnk9Q291bnRyeS5SZWdpb24pDQogICMjIGNvbnZlcnQgZnJvbSB3aWRlIHRvIGxvbmcgZm9ybWF0DQogIGRhdGEgJTw+JSBnYXRoZXIoa2V5PWRhdGUsIHZhbHVlPWNvdW50LCAtY291bnRyeSkNCiAgIyMgY29udmVydCBmcm9tIGNoYXJhY3RlciB0byBkYXRlDQogIGRhdGEgJTw+JSBtdXRhdGUoZGF0ZSA9IGRhdGUgJT4lIG1keSgpKQ0KICAjIyBhZ2dyZWdhdGUgYnkgY291bnRyeQ0KICBkYXRhICU8PiUgZ3JvdXBfYnkoY291bnRyeSwgZGF0ZSkgJT4lIHN1bW1hcmlzZShjb3VudD1zdW0oY291bnQsIG5hLnJtPVQpKSAlPiUgYXMuZGF0YS5mcmFtZSgpDQogIHJldHVybihkYXRhKQ0KfQ0KIyMgY2xlYW4gdGhlIHRocmVlIGRhdGEgc2V0cw0KZGF0YS5jb25maXJtZWQgPC0gZGZfY29uZiAlPiUgY2xlYW5EYXRhKCkgJT4lIHJlbmFtZShjb25maXJtZWQ9Y291bnQpDQpkYXRhLmRlYXRocyA8LSBkZl9kZWF0aHMgJT4lIGNsZWFuRGF0YSgpICU+JSByZW5hbWUoZGVhdGhzPWNvdW50KQ0KZGF0YS5yZWNvdmVyZWQgPC0gZGZfcmVjb3ZlciAlPiUgY2xlYW5EYXRhKCkgJT4lIHJlbmFtZShyZWNvdmVyZWQ9Y291bnQpDQpkYXRhIDwtIGRhdGEuY29uZmlybWVkICU+JSBtZXJnZShkYXRhLmRlYXRocywgYWxsPVQpICU+JSBtZXJnZShkYXRhLnJlY292ZXJlZCwgYWxsPVQpDQpkYXRhDQojIyBjb3VudHJpZXMvcmVnaW9ucyB3aXRoIGNvbmZpcm1lZCBjYXNlcywgZXhjbC4gY3J1aXNlIHNoaXBzDQpjb3VudHJpZXMgPC0gZGF0YSAlPiUgcHVsbChjb3VudHJ5KSAlPiUgc2V0ZGlmZignQ3J1aXNlIFNoaXAnKQ0KZGF0YQ0KYGBgDQoNCmBgYHtyfQ0KZGF0YS53b3JsZCA8LSBkYXRhICU+JSBncm91cF9ieShkYXRlKSAlPiUNCiAgc3VtbWFyaXNlKGNvdW50cnk9J1dvcmxkJywNCiAgICAgICAgICAgIGNvbmZpcm1lZCA9IHN1bShjb25maXJtZWQsIG5hLnJtPVQpLA0KICAgICAgICAgICAgZGVhdGhzID0gc3VtKGRlYXRocywgbmEucm09VCksDQogICAgICAgICAgICByZWNvdmVyZWQgPSBzdW0ocmVjb3ZlcmVkLCBuYS5ybT1UKSkNCmRhdGEgJTw+JSByYmluZChkYXRhLndvcmxkKQ0KZGF0YQ0KZGF0YSAlPD4lIG11dGF0ZShjdXJyZW50LmNvbmZpcm1lZCA9IGNvbmZpcm1lZCAtIGRlYXRocyAtIHJlY292ZXJlZCkNCiNWaWV3KGRhdGEpDQpgYGANCg0KDQoNCmBgYHtyfQ0KIyBXb3JsZCBjb25maXJtZWQgY2FzZXMgbWFwDQp3b3JsZC5sb2NhdGlvbiA8LSBkZl9jb25mDQojVmlldyh3b3JsZC5sb2NhdGlvbiApDQp3b3JsZC5sb2NhdGlvbiRjb25maXJtZWQgPC0gd29ybGQubG9jYXRpb24gWywgbmNvbCh3b3JsZC5sb2NhdGlvbiApXQ0Kd29ybGQubG9jYXRpb24gICU8PiUgc2VsZWN0KGMoQ291bnRyeS5SZWdpb24sIFByb3ZpbmNlLlN0YXRlLCBMYXQsIExvbmcsIGNvbmZpcm1lZCkpICU+JQ0KICBtdXRhdGUodHh0PXBhc3RlMChDb3VudHJ5LlJlZ2lvbiwgJyAtICcsIFByb3ZpbmNlLlN0YXRlLCAnOiAnLCBjb25maXJtZWQpKQ0KbWFwLmNvbmZpcm1lZCA8LSBsZWFmbGV0KHdpZHRoPTEyMDAsIGhlaWdodD04MDApICU+JSBhZGRUaWxlcygpDQojIGNpcmNsZSBtYXJrZXIgKHVuaXRzIGluIHBpeGVscykNCm1hcC5jb25maXJtZWQgJTw+JSBhZGRDaXJjbGVNYXJrZXJzKHdvcmxkLmxvY2F0aW9uJExvbmcsIHdvcmxkLmxvY2F0aW9uJExhdCwNCiAgICAgICAgICAgICAgICAgICAgICAgICMgcmFkaXVzPTIrbG9nMih4JGNvbmZpcm1lZCksDQogICAgICAgICAgICAgICAgICAgICAgICByYWRpdXM9MC4wMipzcXJ0KHdvcmxkLmxvY2F0aW9uJGNvbmZpcm1lZCksDQogICAgICAgICAgICAgICAgICAgICAgICBzdHJva2U9RiwNCiAgICAgICAgICAgICAgICAgICAgICAgIGNvbG9yPSdyZWQnLCBmaWxsT3BhY2l0eT0wLjMsDQogICAgICAgICAgICAgICAgICAgICAgICBsYWJlbD13b3JsZC5sb2NhdGlvbiR0eHQpDQojIHdvcmxkDQptYXAuY29uZmlybWVkDQpgYGANCg0KYGBge3J9DQojIyBDaGluYQ0KbWFwLmNvbmZpcm1lZCAlPiUgc2V0Vmlldyg5NSwgMzUsIHpvb209NCkNCiMjIEF1c3RyYWxpYSBhbmQgTmV3IFplYWxhbmQNCm1hcC5jb25maXJtZWQgJT4lIHNldFZpZXcoMTM1LCAtMjcsIHpvb209NCkNCiMjIFVTIGFuZCBDYW5hZGENCm1hcC5jb25maXJtZWQgJT4lIHNldFZpZXcoLTEwNSwgNDAsIHpvb209NCkNCiMjIEV1cm9wZQ0KbWFwLmNvbmZpcm1lZCAlPiUgc2V0VmlldygxMCwgNTAsIHpvb209NCkNCmBgYA0KDQpgYGB7cn0NCiMgbmV3LmNvbmZpcm1lZCBuZXcuZGVhdGhzIG5ldy5yZWNvdmVyZWQNCmRhdGEgJTw+JSBhcnJhbmdlKGNvdW50cnksIGRhdGUpDQpuIDwtIG5yb3coZGF0YSkNCmRheTEgPC0gbWluKGRhdGEkZGF0ZSkNCmRhdGEgJTw+JSBtdXRhdGUobmV3LmNvbmZpcm1lZCA9IGlmZWxzZShkYXRlID09IGRheTEsIE5BLCBjb25maXJtZWQgLSBsYWcoY29uZmlybWVkLCBuPTEpKSwNCiAgICAgICAgICAgICAgICAgbmV3LmRlYXRocyA9IGlmZWxzZShkYXRlID09IGRheTEsIE5BLCBkZWF0aHMgLSBsYWcoZGVhdGhzLCBuPTEpKSwNCiAgICAgICAgICAgICAgICAgbmV3LnJlY292ZXJlZCA9IGlmZWxzZShkYXRlID09IGRheTEsIE5BLCByZWNvdmVyZWQgLSBsYWcocmVjb3ZlcmVkLCBuPTEpKSkNCmRhdGEgJTw+JSBtdXRhdGUobmV3LmNvbmZpcm1lZCA9IGlmZWxzZShuZXcuY29uZmlybWVkIDwgMCwgMCwgbmV3LmNvbmZpcm1lZCksDQogICAgICAgICAgICAgICAgIG5ldy5kZWF0aHMgPSBpZmVsc2UobmV3LmRlYXRocyA8IDAsIDAsIG5ldy5kZWF0aHMpLA0KICAgICAgICAgICAgICAgICBuZXcucmVjb3ZlcmVkID0gaWZlbHNlKG5ldy5yZWNvdmVyZWQgPCAwLCAwLCBuZXcucmVjb3ZlcmVkKSkNCiMjIGRlYXRoIHJhdGUgYmFzZWQgb24gdG90YWwgZGVhdGhzIGFuZCByZWNvdmVyZWQgY2FzZXMNCmRhdGEgJTw+JSBtdXRhdGUocmF0ZS51cHBlciA9ICgxMDAgKiBkZWF0aHMgLyAoZGVhdGhzICsgcmVjb3ZlcmVkKSkgJT4lIHJvdW5kKDEpLA0KICAgICAgICAgICAgICAgICByYXRlLnVwcGVyID0gaWZlbHNlKGlzLm5hbihyYXRlLnVwcGVyKSwgMCwgcmF0ZS51cHBlcikpDQoNCiMjIGxvd2VyIGJvdW5kOiBkZWF0aCByYXRlIGJhc2VkIG9uIHRvdGFsIGNvbmZpcm1lZCBjYXNlcw0KZGF0YSAlPD4lIG11dGF0ZShyYXRlLmxvd2VyID0gKDEwMCAqIGRlYXRocyAvIGNvbmZpcm1lZCkgJT4lIHJvdW5kKDEpLA0KICAgICAgICAgICAgICAgICByYXRlLmxvd2VyID0gaWZlbHNlKGlzLm5hbihyYXRlLmxvd2VyKSwgMCwgcmF0ZS5sb3dlcikpDQoNCiMjIGRlYXRoIHJhdGUgYmFzZWQgb24gdGhlIG51bWJlciBvZiBkZWF0aC9yZWNvdmVyZWQgb24gZXZlcnkgc2luZ2xlIGRheQ0KZGF0YSAlPD4lIG11dGF0ZShyYXRlLmRhaWx5ID0gKDEwMCAqIG5ldy5kZWF0aHMgLyAobmV3LmRlYXRocyArIG5ldy5yZWNvdmVyZWQpKSAlPiUgcm91bmQoMSksDQogICAgICAgICAgICAgICAgIHJhdGUuZGFpbHkgPSBpZmVsc2UoaXMubmFuKHJhdGUuZGFpbHkpLCAwLCByYXRlLmRhaWx5KSkNCg0KI1ZpZXcoZGF0YSkNCmBgYA0KDQpgYGB7cn0NCiMjIGNvbnZlcnQgZnJvbSB3aWRlIHRvIGxvbmcgZm9ybWF0DQpkYXRhLmxvbmcgPC0gZGF0YSAlPiUNCiAgc2VsZWN0KGMoY291bnRyeSwgZGF0ZSwgY29uZmlybWVkLCBjdXJyZW50LmNvbmZpcm1lZCwgcmVjb3ZlcmVkLCBkZWF0aHMpKSAlPiUNCiAgZ2F0aGVyKGtleT10eXBlLCB2YWx1ZT1jb3VudCwgLWMoY291bnRyeSwgZGF0ZSkpDQojIyBzZXQgZmFjdG9yIGxldmVscyB0byBzaG93IHRoZW0gaW4gYSBkZXNpcmFibGUgb3JkZXINCmRhdGEubG9uZyAlPD4lIG11dGF0ZSh0eXBlPXJlY29kZV9mYWN0b3IodHlwZSwgY29uZmlybWVkPSdUb3RhbCBDb25maXJtZWQnLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjdXJyZW50LmNvbmZpcm1lZD0nQ3VycmVudCBDb25maXJtZWQnLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICByZWNvdmVyZWQ9J1JlY292ZXJlZCcsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRlYXRocz0nRGVhdGhzJykpDQojVmlldyhkYXRhLmxvbmcpDQpgYGANCg0KYGBge3J9DQojI051bWJlciBvZiBjYXNlIFdvcmxkDQp3b3JsZCA8LSBmaWx0ZXIoZGF0YS5sb25nLGNvdW50cnkgPT0gJ1dvcmxkJykNCndvcmxkDQpwbG90X3dvcmxkLmFyZWEgPC0gd29ybGQgJT4lIGZpbHRlcih0eXBlICE9ICdUb3RhbCBDb25maXJtZWQnKSAlPiUNCiAgZ2dwbG90KGFlcyh4PWRhdGUsIHk9Y291bnQpKSArDQogIGdlb21fYXJlYShhZXMoZmlsbD10eXBlKSwgYWxwaGE9MC41KSArDQogIGxhYnModGl0bGU9cGFzdGUwKCdOdW1iZXJzIG9mIENhc2VzIFdvcmxkd2lkZSAtICcsIG1heC5kYXRlLnR4dCkpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoJ3JlZCcsICdncmVlbicsICdibGFjaycpKSArDQogIHRoZW1lKGxlZ2VuZC50aXRsZT1lbGVtZW50X2JsYW5rKCksIGxlZ2VuZC5wb3NpdGlvbj0nYm90dG9tJywNCiAgICAgICAgcGxvdC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplPTcpLA0KICAgICAgICBheGlzLnRpdGxlLng9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBheGlzLnRpdGxlLnk9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBsZWdlbmQua2V5LnNpemU9dW5pdCgwLjIsICdjbScpLA0KICAgICAgICBsZWdlbmQudGV4dD1lbGVtZW50X3RleHQoc2l6ZT02KSwNCiAgICAgICAgYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTcpLA0KICAgICAgICBheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKQ0KcGxvdF93b3JsZC5saW5lIDwtIHdvcmxkICU+JQ0KICBnZ3Bsb3QoYWVzKHg9ZGF0ZSwgeT1jb3VudCxjb2xvciA9IHR5cGUpKSArDQogIGdlb21fbGluZShhZXMoY29sb3I9dHlwZSkpICsNCiAgbGFicyh0aXRsZT1wYXN0ZTAoJ051bWJlcnMgb2YgQ2FzZXMgV29ybGR3aWRlIChsb2cgc2NhbGUpIC0gJywgbWF4LmRhdGUudHh0KSkgKw0KICAjc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcz1jKCdibHVlJywgJ3JlZCcsICdncmVlbicsICdibGFjaycpKSArDQogIHRoZW1lKGxlZ2VuZC50aXRsZT1lbGVtZW50X2JsYW5rKCksIGxlZ2VuZC5wb3NpdGlvbj0nYm90dG9tJywNCiAgICAgICAgcGxvdC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplPTcpLA0KICAgICAgICBheGlzLnRpdGxlLng9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBheGlzLnRpdGxlLnk9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBsZWdlbmQua2V5LnNpemU9dW5pdCgwLjIsICdjbScpLA0KICAgICAgICBsZWdlbmQudGV4dD1lbGVtZW50X3RleHQoc2l6ZT02KSwNCiAgICAgICAgYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTcpLA0KICAgICAgICBheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKSArDQogIHNjYWxlX3lfY29udGludW91cyh0cmFucz0nbG9nMTAnKQ0KIyMgc2hvdyB0d28gcGxvdHMgc2lkZSBieSBzaWRlDQpncmlkLmFycmFuZ2UocGxvdF93b3JsZC5hcmVhLCBwbG90X3dvcmxkLmxpbmUsIG5jb2w9MikNCg0KZ2x5LnBsb3Rfd29ybGQuYXJlYSA8LSBnZ3Bsb3RseShwbG90X3dvcmxkLmFyZWEpDQpnbHkucGxvdF93b3JsZC5hcmVhDQoNCmdseS5wbG90X3dvcmxkLmxpbmUgPC0gZ2dwbG90bHkocGxvdF93b3JsZC5saW5lKQ0KZ2x5LnBsb3Rfd29ybGQubGluZQ0KYGBgDQoNCg0KDQpgYGB7cn0NCmRhdGEud29ybGQuYWxsIDwtIGRhdGEgJT4lIGZpbHRlcihjb3VudHJ5ID09ICJXb3JsZCIpDQojIyBhIHNjYXR0ZXIgcGxvdCB3aXRoIGEgc21vb3RoZWQgbGluZSBhbmQgdmVydGljYWwgeC1heGlzIGxhYmVscw0KcGxvdF9hY2MuZGVhdGhzIDwtIGdncGxvdChkYXRhLndvcmxkLmFsbCwgYWVzKHg9ZGF0ZSwgeT1kZWF0aHMpKSArDQogIGdlb21fcG9pbnQoKSArIGdlb21fc21vb3RoKCkgKw0KICB4bGFiKCcnKSArIHlsYWIoJ0NvdW50JykgKyBsYWJzKHRpdGxlPSdBY2N1bXVsYXRpdmUgRGVhdGhzJykgKw0KICB0aGVtZShheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKQ0KcGxvdF9yZWNvdmVyZWQgPC0gZ2dwbG90KGRhdGEud29ybGQuYWxsLCBhZXMoeD1kYXRlLCB5PXJlY292ZXJlZCkpICsNCiAgZ2VvbV9wb2ludCgpICsgZ2VvbV9zbW9vdGgoKSArDQogIHhsYWIoJycpICsgeWxhYignQ291bnQnKSArIGxhYnModGl0bGU9J0FjY3VtdWxhdGl2ZSBSZWNvdmVyZWQgQ2FzZXMnKSArDQogIHRoZW1lKGF4aXMudGV4dC54PWVsZW1lbnRfdGV4dChhbmdsZT00NSwgaGp1c3Q9MSkpDQpwbG90X25ldy5kZWF0aHMgPC0gZ2dwbG90KGRhdGEud29ybGQuYWxsLCBhZXMoeD1kYXRlLCB5PW5ldy5kZWF0aHMpKSArDQogIGdlb21fcG9pbnQoKSArIGdlb21fc21vb3RoKCkgKw0KICB4bGFiKCcnKSArIHlsYWIoJ0NvdW50JykgKyBsYWJzKHRpdGxlPSdOZXcgRGVhdGhzJykgKw0KICB0aGVtZShheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKQ0KcGxvdF9uZXcucmVjb3ZlcmVkIDwtIGdncGxvdChkYXRhLndvcmxkLmFsbCwgYWVzKHg9ZGF0ZSwgeT1uZXcucmVjb3ZlcmVkKSkgKw0KICBnZW9tX3BvaW50KCkgKyBnZW9tX3Ntb290aCgpICsNCiAgeGxhYignJykgKyB5bGFiKCdDb3VudCcpICsgbGFicyh0aXRsZT0nTmV3IFJlY292ZXJlZCBDYXNlcycpICsNCiAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdD0xKSkNCiMjIHNob3cgZm91ciBwbG90cyB0b2dldGhlciwgd2l0aCAyIHBsb3RzIGluIGVhY2ggcm93DQpncmlkLmFycmFuZ2UocGxvdF9hY2MuZGVhdGhzLCBwbG90X3JlY292ZXJlZCwgcGxvdF9uZXcuZGVhdGhzLCBwbG90X25ldy5yZWNvdmVyZWQsIG5yb3c9MikNCg0KYGBgDQoNCmBgYHtyfQ0KIyMgY29udmVydCBmcm9tIHdpZGUgdG8gbG9uZyBmb3JtYXQsIGZvciBkcmF3aW5nIGFyZWEgcGxvdHMNCnJhdGVzLmxvbmcgPC0gZGF0YSAlPiUNCiAgc2VsZWN0KGMoY291bnRyeSwgZGF0ZSwgcmF0ZS51cHBlciwgcmF0ZS5sb3dlciwgcmF0ZS5kYWlseSkpICU+JQ0KICBnYXRoZXIoa2V5PXR5cGUsIHZhbHVlPWNvdW50LCAtYyhjb3VudHJ5LCBkYXRlKSkNCiMgc2V0IGZhY3RvciBsZXZlbHMgdG8gc2hvdyB0aGVtIGluIGEgZGVzaXJhYmxlIG9yZGVyDQpyYXRlcy5sb25nICU8PiUgbXV0YXRlKHR5cGU9cmVjb2RlX2ZhY3Rvcih0eXBlLCByYXRlLmRhaWx5PSdEYWlseScsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgIHJhdGUudXBwZXI9J1VwcGVyIGJvdW5kJywNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcmF0ZS5sb3dlciA9ICdMb3dlciBib3VuZCcpKSANCiNWaWV3KHJhdGVzLmxvbmcpIA0KDQpwbG90X3JhdGVzIDwtIHJhdGVzLmxvbmcgJT4lIGZpbHRlcihjb3VudHJ5ID09ICJXb3JsZCIpICU+JSANCiAgZ2dwbG90KGFlcyh4ID0gZGF0ZSwgeSA9IGNvdW50LCBjb2xvciA9IHR5cGUpKSArIA0KICBnZW9tX2xpbmUoKSArDQogIGxhYnModGl0bGUgPSAiV29ybGQgRGVhdGggUmF0ZSAoJSkiKSArDQogIHhsYWIoIiIpICsgeWxhYigiRGVhdGggUmF0ZSAoJSkiKQ0KDQpwbG90X3JhdGVzDQoNCmdseS5yYXRlcyA8LSBnZ3Bsb3RseShwbG90X3JhdGVzKQ0KZ2x5LnJhdGVzDQpgYGANCg0KYGBge3J9DQojIyByYW5raW5nIGJ5IGNvbmZpcm1lZCBjYXNlcw0KZGF0YS5sYXRlc3QuYWxsIDwtIGRhdGEgJT4lIGZpbHRlcihkYXRlID09IG1heChkYXRlKSkgJT4lDQogIHNlbGVjdChjb3VudHJ5LCBkYXRlLGNvbmZpcm1lZCwgbmV3LmNvbmZpcm1lZCwgY3VycmVudC5jb25maXJtZWQsDQogICAgICAgICByZWNvdmVyZWQsIGRlYXRocywgbmV3LmRlYXRocywgZGVhdGgucmF0ZT1yYXRlLmxvd2VyKSAlPiUNCiAgbXV0YXRlKHJhbmtpbmcgPSBkZW5zZV9yYW5rKGRlc2MoY29uZmlybWVkKSkpICU+JQ0KICBhcnJhbmdlKHJhbmtpbmcpDQojVmlldyhkYXRhLmxhdGVzdC5hbGwpDQoNCmsgPC0gMjANCiMjIHRvcCAyMCBjb3VudHJpZXM6IDIxIGluY2wuICdXb3JsZCcNCnRvcC5jb3VudHJpZXMgPC0gZGF0YS5sYXRlc3QuYWxsICU+JSBmaWx0ZXIocmFua2luZyA8PSBrICsgMSkgJT4lDQogIGFycmFuZ2UocmFua2luZykgJT4lIHB1bGwoY291bnRyeSkgJT4lIGFzLmNoYXJhY3RlcigpDQp0b3AuY291bnRyaWVzICU+JSBzZXRkaWZmKCdXb3JsZCcpICU+JSBwcmludCgpDQpgYGANCg0KDQpgYGB7cn0NCmRmIDwtIGRhdGEubG9uZyAlPiUgZmlsdGVyKGNvdW50cnkgJWluJSB0b3AuY291bnRyaWVzKSAlPiUNCiAgbXV0YXRlKGNvdW50cnk9Y291bnRyeSAlPiUgZmFjdG9yKGxldmVscz1jKHRvcC5jb3VudHJpZXMpKSkgJT4lIA0KICBmaWx0ZXIoY291bnRyeSAhPSAnV29ybGQnICYgdHlwZSAhPSAnVG90YWwgQ29uZmlybWVkJykgJT4lDQogIGdncGxvdChhZXMoeD1kYXRlLCB5PWNvdW50LCBjb2xvcj10eXBlKSkgKw0KICBnZW9tX2xpbmUoKSsNCiAgI2dlb21fYXJlYShhbHBoYT0wLjUpICsNCiMgeGxhYignJykgKyB5bGFiKCcnKSArDQogIGxhYnModGl0bGU9cGFzdGUwKCdOdW1iZXJzIG9mIENPVklELTE5IENhc2VzIGluIFRvcCAyMCBDb3VudHJpZXMgLSAnLA0KICAgICAgICAgICAgICAgICAgICBtYXguZGF0ZS50eHQpKSArDQogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXM9YygncmVkJywgJ2dyZWVuJywgJ2JsYWNrJykpICsNCiAgdGhlbWUobGVnZW5kLnRpdGxlPWVsZW1lbnRfYmxhbmsoKSwgbGVnZW5kLnBvc2l0aW9uPSdib3R0b20nLA0KICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTIpLA0KICAgICAgICBheGlzLnRpdGxlLng9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBheGlzLnRpdGxlLnk9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBsZWdlbmQua2V5LnNpemU9dW5pdCgwLjQsICdjbScpLA0KICAgICAgICBsZWdlbmQudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMiksDQogICAgICAgIHN0cmlwLnRleHQueD1lbGVtZW50X3RleHQoc2l6ZT0xMiksDQogICAgICAgIGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMiksDQogICAgICAgIGF4aXMudGV4dC54PWVsZW1lbnRfdGV4dChhbmdsZT00NSwgaGp1c3Q9MSkpICsNCiAgZmFjZXRfd3JhcCh+Y291bnRyeSwgbmNvbD00LCBzY2FsZXM9J2ZyZWVfeScpICsNCiAgc2NhbGVfeV9jb250aW51b3VzKHRyYW5zPSdsb2cxMCcpDQpkZg0KYGBgDQoNCmBgYHtyfQ0KZGF0YS5sYXRlc3QgPC0gZGF0YS5sYXRlc3QuYWxsICU+JSBmaWx0ZXIoIWlzLm5hKGNvdW50cnkpKSAlPiUNCiAgbXV0YXRlKGNvdW50cnk9aWZlbHNlKHJhbmtpbmcgPD0gayArIDEsIGFzLmNoYXJhY3Rlcihjb3VudHJ5KSwgJ090aGVycycpKSAlPiUNCiAgbXV0YXRlKGNvdW50cnk9Y291bnRyeSAlPiUgZmFjdG9yKGxldmVscz1jKHRvcC5jb3VudHJpZXMsICdPdGhlcnMnKSkpDQpkYXRhLmxhdGVzdCAlPD4lIGdyb3VwX2J5KGNvdW50cnkpICU+JQ0KICBzdW1tYXJpc2UoY29uZmlybWVkPXN1bShjb25maXJtZWQpLCBuZXcuY29uZmlybWVkPXN1bShuZXcuY29uZmlybWVkKSwNCiAgICAgICAgICAgIGN1cnJlbnQuY29uZmlybWVkPXN1bShjdXJyZW50LmNvbmZpcm1lZCksDQogICAgICAgICAgICByZWNvdmVyZWQ9c3VtKHJlY292ZXJlZCksIGRlYXRocz1zdW0oZGVhdGhzKSwgbmV3LmRlYXRocz1zdW0obmV3LmRlYXRocykpICU+JQ0KICBtdXRhdGUoZGVhdGgucmF0ZT0oMTAwICogZGVhdGhzL2NvbmZpcm1lZCkgJT4lIHJvdW5kKDEpKSANCmRhdGEubGF0ZXN0DQpkYXRhLmxhdGVzdCAlPD4lIHNlbGVjdChjKGNvdW50cnksIGNvbmZpcm1lZCwgZGVhdGhzLCBkZWF0aC5yYXRlLA0KICAgICAgICAgICAgICAgICAgICAgICAgICBuZXcuY29uZmlybWVkLCBuZXcuZGVhdGhzLCBjdXJyZW50LmNvbmZpcm1lZCxyZWNvdmVyZWQpKSAlPiUNCiAgbXV0YXRlKHJlY292ZXIucmF0ZT0oMTAwICogcmVjb3ZlcmVkL2NvbmZpcm1lZCkgJT4lIHJvdW5kKDEpKQ0KZGF0YS5sYXRlc3QNCg0KZGZfcG9wIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIHBvcHVsYXRpb24gIikpDQpkZl9wb3AgPC0gYXMuZGF0YS5mcmFtZShkZl9wb3ApDQpkZl9wb3AgPC0gcmVuYW1lKGRmX3BvcCwiY291bnRyeSI9IkNvdW50cnkiKQ0KDQojIEFkZCBXb3JsZCBQb3B1bGF0aW9uDQp3b3JsZF9wb3AgPC0gc3VtKGRmX3BvcCRgUG9wdWxhdGlvbiAoMjAyMClgKQ0KZGZfcG9wW25yb3coZGZfcG9wKSArIDEsXSA9IGMoIldvcmxkIiwgd29ybGRfcG9wKQ0KDQojIEFkZCBPdGhlciBDb3VudHJpZXMgUG9wdWxhdGlvbg0KdG9wX3BvcCA8LSBmaWx0ZXIoZGZfcG9wLCBkZl9wb3AkY291bnRyeSAlaW4lIHRvcC5jb3VudHJpZXMgJiBkZl9wb3AkY291bnRyeSAhPSAiV29ybGQiKQ0KdG9wX3BvcCA8LSBzdW0odG9wX3BvcCRgUG9wdWxhdGlvbiAoMjAyMClgICU+JSBhcy5udW1lcmljKCkpDQpvdGhlcnNfcG9wIDwtICh3b3JsZF9wb3AgLSB0b3BfcG9wKSANCmRmX3BvcFtucm93KGRmX3BvcCkgKyAxLF0gPSBjKCJPdGhlcnMiLCBvdGhlcnNfcG9wKQ0KI1ZpZXcoZGZfcG9wKQ0KDQpkYXRhLmxhdGVzdCA8LSBtZXJnZSh4ID0gZGF0YS5sYXRlc3QsIHkgPSBkZl9wb3AsIGJ5ID0gImNvdW50cnkiLCBhbGwueCA9IFRSVUUpIA0KZGF0YS5sYXRlc3QNCmRhdGEubGF0ZXN0IDwtIHJlbmFtZShkYXRhLmxhdGVzdCwicG9wdWxhdGlvbiIgPSAiUG9wdWxhdGlvbiAoMjAyMCkiKQ0KZGF0YS5sYXRlc3QkcG9wdWxhdGlvbiA8LSBkYXRhLmxhdGVzdCRwb3B1bGF0aW9uICU+JSBhcy5udW1lcmljKCkNCmRhdGEubGF0ZXN0ICA8LSBkYXRhLmxhdGVzdCAlPiUNCiAgc2VsZWN0KGMoY291bnRyeSwgY29uZmlybWVkLCBkZWF0aHMsIGRlYXRoLnJhdGUsDQogICAgICAgICAgICAgICAgICAgICAgICAgIG5ldy5jb25maXJtZWQsIG5ldy5kZWF0aHMsDQogICAgICAgICAgICAgICAgICAgICAgICAgIGN1cnJlbnQuY29uZmlybWVkLCByZWNvdmVyZWQsIHJlY292ZXIucmF0ZSwgcG9wdWxhdGlvbikpICU+JQ0KICBtdXRhdGUoY29uZmlybS5yYXRlID0gKDEwMCAqIGNvbmZpcm1lZCAvIHBvcHVsYXRpb24pICU+JSByb3VuZCgxKSkNCmRhdGEubGF0ZXN0DQpgYGANCg0KYGBge3J9DQojIyAlIG9mIGRlYXRoDQpkYXRhLmxhdGVzdCAlPiUgbXV0YXRlKGRlYXRoLnJhdGU9ZGVhdGgucmF0ZSAlPiUgZm9ybWF0KG5zbWFsbD0xKSAlPiUgcGFzdGUwKCclJykpDQoNCiMjIGNvbnZlcnQgZnJvbSB3aWRlIHRvIGxvbmcgZm9ybWF0LCBmb3IgZHJhd2luZyBhcmVhIHBsb3RzDQpkYXRhLmxhdGVzdC5sb25nIDwtIGRhdGEubGF0ZXN0ICU+JSBmaWx0ZXIoY291bnRyeSE9J1dvcmxkJykgJT4lDQogIGdhdGhlcihrZXk9dHlwZSwgdmFsdWU9Y291bnQsIC1jb3VudHJ5KQ0KIyMgc2V0IGZhY3RvciBsZXZlbHMgdG8gc2hvdyB0aGVtIHdpdGggcHJvcGVyIHRleHQgYW5kIGluIGEgZGVzaXJhYmxlIG9yZGVyDQpkYXRhLmxhdGVzdC5sb25nICU8PiUgbXV0YXRlKHR5cGU9cmVjb2RlX2ZhY3Rvcih0eXBlLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29uZmlybWVkPSdUb3RhbCBDb25maXJtZWQnLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZGVhdGhzPSdUb3RhbCBEZWF0aHMnLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZGVhdGgucmF0ZT0nRGVhdGggUmF0ZSAoJSknLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbmV3LmNvbmZpcm1lZD0nTmV3IENvbmZpcm1lZCAoY29tcGFyZWQgd2l0aCBvbmUgZGF5IGJlZm9yZSknLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbmV3LmRlYXRocz0nTmV3IERlYXRocyAoY29tcGFyZWQgd2l0aCBvbmUgZGF5IGJlZm9yZSknLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY3VycmVudC5jb25maXJtZWQ9J0N1cnJlbnQgQ29uZmlybWVkJywNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHJlY292ZXIucmF0ZSA9ICdSZWNvdmVyZWQgUmF0ZSglKScsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb25maXJtLnJhdGUgPSAnQ29uZmlybWVkIFJhdGUoJSknKSkNCiNWaWV3KGRhdGEubGF0ZXN0LmxvbmcpDQpkYXRhLm9uZS5kZW0gPC0gZmlsdGVyKGRhdGEubGF0ZXN0LmxvbmcsdHlwZT09J1RvdGFsIENvbmZpcm1lZCcNCiAgICAgICAgICAgICAgICAgICAgICAgfCB0eXBlPT0nVG90YWwgRGVhdGhzJw0KICAgICAgICAgICAgICAgICAgICAgICB8IHR5cGU9PSdDdXJyZW50IENvbmZpcm1lZCcpDQpkYXRhLnR3by5kZW0gPC0gZmlsdGVyKGRhdGEubGF0ZXN0LmxvbmcsdHlwZT09J0RlYXRoIFJhdGUgKCUpJw0KICAgICAgICAgICAgICAgICAgICAgIyAgfCB0eXBlPT0nTmV3IENvbmZpcm1lZCAoY29tcGFyZWQgd2l0aCBvbmUgZGF5IGJlZm9yZSknDQogICAgICAgICAgICAgICAgICAgICMgICB8IHR5cGU9PSdOZXcgRGVhdGhzIChjb21wYXJlZCB3aXRoIG9uZSBkYXkgYmVmb3JlKScNCiAgICAgICAgICAgICAgICAgICAgICAgfCB0eXBlPT0nUmVjb3ZlcmVkIFJhdGUoJSknDQogICAgICAgICAgICAgICAgICAgICAgIHwgdHlwZT09J0NvbmZpcm1lZCBSYXRlKCUpJykNCmRhdGEudHdvLmRlbQ0KYGBgDQoNCmBgYHtyfQ0KIyMgYmFyIGNoYXJ0DQpkYXRhLm9uZS5kZW0gJT4lIGdncGxvdChhZXMoeD1jb3VudHJ5LCB5PWNvdW50LCBmaWxsPWNvdW50cnksIGdyb3VwPWNvdW50cnkpKSArDQogIGdlb21fYmFyKHN0YXQ9J2lkZW50aXR5JykgKw0KICBnZW9tX3RleHQoYWVzKGxhYmVsPWNvdW50LCB5PWNvdW50KSwgc2l6ZT0zLCB2anVzdD0wKSArDQogIHhsYWIoJycpICsgeWxhYignJykgKw0KICBsYWJzKHRpdGxlPXBhc3RlMCgnVG9wIDIwIENvdW50cmllcyB3aXRoIE1vc3QgQ29uZmlybWVkIENhc2VzIC0gJywgbWF4LmRhdGUudHh0KSkgKw0KICBzY2FsZV9maWxsX2Rpc2NyZXRlKG5hbWU9J0NvdW50cnknLCBsYWJlbHM9YWVzKGNvdW50KSkgKw0KICB0aGVtZShsZWdlbmQudGl0bGU9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBsZWdlbmQucG9zaXRpb249J25vbmUnLA0KICAgICAgICBwbG90LnRpdGxlPWVsZW1lbnRfdGV4dChzaXplPTEzKSwNCiAgICAgICAgYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTgpLA0KICAgICAgICBheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKSArDQogIGZhY2V0X3dyYXAofnR5cGUsIG5jb2w9MSwgc2NhbGVzPSdmcmVlX3knKQ0KYGBgDQoNCmBgYHtyfQ0KZGF0YS50d28uZGVtJGZhY2V0IDwtIGZhY3RvcihkYXRhLnR3by5kZW0kdHlwZSwgbGV2ZWxzID0gYygnQ29uZmlybWVkIFJhdGUoJSknLCAnUmVjb3ZlcmVkIFJhdGUoJSknLCdEZWF0aCBSYXRlICglKScpKQ0KZGF0YS50d28uZGVtICU+JSANCiAgZ2dwbG90KGFlcyh4PWNvdW50cnksIHk9Y291bnQsIGZpbGw9Y291bnRyeSwgZ3JvdXA9Y291bnRyeSkpICsNCiAgZ2VvbV9iYXIoc3RhdD0naWRlbnRpdHknKSArDQogIGdlb21fdGV4dChhZXMobGFiZWw9Y291bnQsIHk9Y291bnQpLCBzaXplPTQsIHZqdXN0PTApICsNCiAgeGxhYignJykgKyB5bGFiKCcnKSArDQogIGxhYnModGl0bGU9cGFzdGUwKCdUb3AgMjAgQ291bnRyaWVzIHdpdGggTW9zdCBDb25maXJtZWQgQ2FzZXMgLSAnLCBtYXguZGF0ZS50eHQpKSArDQogIHNjYWxlX2ZpbGxfZGlzY3JldGUobmFtZT0nQ291bnRyeScsIGxhYmVscz1hZXMoY291bnQpKSArDQogIHRoZW1lKGxlZ2VuZC50aXRsZT1lbGVtZW50X2JsYW5rKCksDQogICAgICAgIGxlZ2VuZC5wb3NpdGlvbj0nbm9uZScsDQogICAgICAgIHBsb3QudGl0bGU9ZWxlbWVudF90ZXh0KHNpemU9MTMpLA0KICAgICAgICBheGlzLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9MTApLA0KICAgICAgICBheGlzLnRleHQueD1lbGVtZW50X3RleHQoc2l6ZT0xMCxhbmdsZT00NSwgaGp1c3Q9MSkpICsNCiAgZmFjZXRfd3JhcCh+ZmFjZXQsIG5jb2w9MSwgc2NhbGVzPSdmcmVlX3knKQ0KYGBgDQoNCmBgYHtyfQ0KIyNHRFANCmRmX2dkcDIwMTkgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gZ2RwMTkiKSkNCmRmX2dkcDIwMTkgPC0gYXMuZGF0YS5mcmFtZShkZl9nZHAyMDE5KQ0KZGZfZ2RwMjAxOQ0KYGBgDQoNCmBgYHtyfQ0KI2hlYWx0aHJhbmtpbmcgZGF0YQ0KZGZfaGVhbHQgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gaGVhbHRocmFua2luZyIpKQ0KZGZfaGVhbHQgPC0gYXMuZGF0YS5mcmFtZShkZl9oZWFsdCkNCmRmX2hlYWx0IDwtIHNlbGVjdChkZl9oZWFsdCxjKCJjb3VudHJ5IiwiaGVhbHRoQ2FyZUluZGV4IikpDQojVmlldyhkZl9oZWFsdCkNCmBgYA0KDQpgYGB7cn0NCiNUb3AyMFBvcm5odWIgZGF0YQ0KI2RmX3Bvcm5odWIgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gUG9ybmh1YiIpKQ0KI2RmX3Bvcm5odWIgPC0gYXMuZGF0YS5mcmFtZShkZl9wb3JuaHViKQ0KI2RmX3Bvcm5odWINCmBgYA0KDQpgYGB7cn0NCiMgdGVtcCBkYXRhDQpkZl90ZW1wIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIHdvcmxkX3RlbXAiKSkNCmRmX3RlbXAgPC0gYXMuZGF0YS5mcmFtZShkZl90ZW1wKSANCmRmX3RlbXAkQ291bnRyeVtkZl90ZW1wJENvdW50cnkgPT0gIlVuaXRlZCBTdGF0ZXMiXSA8LSAiVVMiDQoNCmRmX2NpdHkgPC0gc2VsZWN0KGRmX3RlbXAsYygiQ291bnRyeSIsIkNpdHkiKSkgJT4lDQogIHJlbmFtZShjb3VudHJ5PUNvdW50cnkpICU+JSANCiAgcmVuYW1lKGNpdHk9Q2l0eSkNCg0KbnVtb2ZjaXR5IDwtIGFnZ3JlZ2F0ZShjaXR5IH4gY291bnRyeSwgZGF0YSA9IGRmX2NpdHksIGxlbmd0aCkNCg0KI2NsZWFuIGRhdGENCmRmX3RlbXAgPC0gc2VsZWN0KGRmX3RlbXAsYygiQ291bnRyeSIsIkF2Z19ZZWFyIikpICU+JQ0KICByZW5hbWUoY291bnRyeT1Db3VudHJ5KQ0KI1ZpZXcoZGZfdGVtcCkNCg0KI2RmX3RlbXAgPC0gZGF0YS5mcmFtZShjb3VudHJ5PWRmX3RlbXBbLDFdLGF2Zz1yb3dNZWFucyhkZl90ZW1wWywtMV0pKQ0KZGZfdGVtcCA8LSBkZl90ZW1wICU8PiUgZ3JvdXBfYnkoY291bnRyeSkgJT4lIHN1bW1hcmlzZShhdmdfdGVtcCA9IG1lYW4oQXZnX1llYXIsbmEucm0gPSBUUlVFKSU+JSByb3VuZCgxKSkNCmRmX3RlbXANCmBgYA0KDQpgYGB7cn0NCiMjdGVtcCBiYXINCmRmX3RlbXAuYWxsIDwtIGRmX3RlbXAgJT4lIG1lcmdlKGRhdGEubGF0ZXN0LmFsbCkNCiNWaWV3KGRmX3RlbXAuYWxsKQ0KZGZfdGVtcF90b3AuYWxsIDwtIGRmX3RlbXAuYWxsICU+JSBmaWx0ZXIoY291bnRyeSAlaW4lIHRvcC5jb3VudHJpZXMpICU+JQ0KICBtdXRhdGUocmFua2luZyA9IHJhbmtpbmcgLSAxKSAlPiUNCiAgYXJyYW5nZShyYW5raW5nKQ0KI1ZpZXcoZGZfdGVtcF90b3AuYWxsKQ0KZ190ZW1wX3RvcCA8LSBkZl90ZW1wX3RvcC5hbGwgJT4lDQogIGdncGxvdChhZXMoeCA9IHJlb3JkZXIoY291bnRyeSwgcmFua2luZyksIHkgPSBhdmdfdGVtcCwgZmlsbCA9IGF2Z190ZW1wKSkgKw0KICBsYWJzKHRpdGxlPXBhc3RlMCgiVGVtcGVyYXR1cmUgaW4gVG9wICAyMCBjb3VudHJpZXMiKSwgc3VidGl0bGUgPSAiQXZlcmFnZSBUZW1wZXJhdHVyZSBpbiBUb3AgMjAgY291bnRyaWVzIHdpdGggbW9zdCBjb25maXJtZWQgY2FzZXMgKMKwQykgKDIwMjApIikgKw0KICBzY2FsZV9jb2xvcl9ncmFkaWVudChsb3cgPSAiIzkzREJGRiIsIGhpZ2ggPSAiI0ZGNzc3MSIpICsNCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbD1hdmdfdGVtcCwgeT1hdmdfdGVtcCksIHNpemU9NCwgdmp1c3Q9LTAuNSkgKw0KICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5IiwgcG9zaXRpb24gPSAiZG9kZ2UiKSArDQogIHRoZW1lKA0KICAgICAgICBsZWdlbmQudGl0bGU9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICBsZWdlbmQucG9zaXRpb249J25vbmUnLA0KICAgICAgICBwbG90LnRpdGxlPWVsZW1lbnRfdGV4dChzaXplID0gMTUsIGhqdXN0ID0gMC41KSwNCiAgICAgICAgcGxvdC5zdWJ0aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTIsIGhqdXN0ID0gMC41KSwNCiAgICAgICAgYXhpcy50ZXh0PWVsZW1lbnRfdGV4dChzaXplPTgpLA0KICAgICAgICBheGlzLnRleHQueD1lbGVtZW50X3RleHQoc2l6ZSA9IDksIGFuZ2xlPTQ1LCBoanVzdD0xKSkgKw0KICBzY2FsZV94X2Rpc2NyZXRlKG5hbWUgPSAiQ291bnRyeSIpICsNCiAgc2NhbGVfeV9kaXNjcmV0ZShuYW1lID0gIkF2ZXJhZ2UgVGVtcGVyYXR1cmUiKQ0KDQogICAgICAgICAgICNsYWJzKHRpdGxlID0gIlRlbXBlcmF0dXJlIGluIFRvcCAgMjAgY291bnRyaWVzIiwgc3VidGl0bGUgPSAiVGVtcGVyYXR1cmUgaW4gVG9wIDIwIGNvdW50cmllcyB3aXRoIG1vc3QgY29uZmlybWVkIGNhc2VzICjCsEMpIikNCiAgICAgICAgIA0KI2dfdGVtcF90b3AgICAgICAgICANCmdfdGVtcF90b3AgDQpgYGANCg0KYGBge3J9DQojZGF0YS5sYXRlc3QuYWxsDQpsYXQubG9uZyA8LSByZW5hbWUoZGZfY29uZiwgImNvdW50cnkiID0gIkNvdW50cnkuUmVnaW9uIiwgImNpdHkiID0gIlByb3ZpbmNlLlN0YXRlIikgJT4lIA0KICBzZWxlY3QoImNvdW50cnkiLCAiTGF0IiwgIkxvbmciKSAlPiUgDQogIG1lcmdlKGRmX3RlbXAuYWxsW2MoImNvdW50cnkiLCJjb25maXJtZWQiLCAicmVjb3ZlcmVkIiwgImRlYXRocyIsICJhdmdfdGVtcCIsICJyYW5raW5nIildLCBieSA9ICJjb3VudHJ5IikgJT4lDQogIGRpc3RpbmN0KGNvdW50cnksIC5rZWVwX2FsbCA9IFRSVUUpICU+JQ0KICBtdXRhdGUocmFua2luZyA9IHJhbmtpbmcgLSAxKSAlPiUNCiAgYXJyYW5nZShyYW5raW5nKQ0KI1ZpZXcobGF0LmxvbmcpDQpgYGANCg0KYGBge3J9DQpsYWJlbF93b3JsZCA8LSBsYXQubG9uZyANCmxhYmVsX3dvcmxkJGF2Z190ZW1wIDwtIGFzLm51bWVyaWMobGFiZWxfd29ybGRbLCBuYW1lcyhsYWJlbF93b3JsZCkgJWluJSBjKCJhdmdfdGVtcCIpXSkNCmxhYmVsX3dvcmxkIDwtIGxhYmVsX3dvcmxkICU+JSAgDQogIG11dGF0ZSh0eHQ9cGFzdGUwKCc8Yj4nLHJhbmtpbmcsICc8L2I+JywNCiAgICAgICAgICAgICAgICAgICAgJzxici8+JywnPGI+Jyxjb3VudHJ5LCAnPC9iPicsDQogICAgICAgICAgICAgICAgICAgICc8YnIvPicsICJUZW1wZXJhdHVyZTogICIsYXZnX3RlbXAsICcgwrBDJywNCiAgICAgICAgICAgICAgICAgICAgJzxici8+JywgIkNvbmZpcm1lZDogICIsIGNvbmZpcm1lZCwgDQogICAgICAgICAgICAgICAgICAgICc8YnIvPicsICJEZWF0aHM6ICIsIGRlYXRocywNCiAgICAgICAgICAgICAgICAgICAgJzxici8+JywgIlJlY292ZXJlZDogIiwgcmVjb3ZlcmVkDQogICAgICAgICAgICAgICAgICAgICkpIA0KDQpsYWJlbF93b3JsZCR0eHQgPC0gbGFiZWxfd29ybGQkdHh0ICU+JSBsYXBwbHkoaHRtbHRvb2xzOjpIVE1MKQ0KbGFiZWxfd29ybGQgDQoNCmxhYmVsX3RvcCA8LSBsYWJlbF93b3JsZCAlPiUgZmlsdGVyKHJhbmtpbmcgPCAyMSkNCmxhYmVsX3RvcA0KYGBgDQoNCg0KYGBge3J9DQojIFRlbXBlcmF0dXJlIE1hcA0Kd3BhbCA8LSBjb2xvck51bWVyaWMoIllsT3JSZCIsIGxhYmVsX3dvcmxkJGF2Z190ZW1wLCBuID0gNCkNCg0KdG9wSWNvbiA8LSBtYWtlSWNvbigic3Rhci5wbmciLA0KICAjaWNvblVybCA9ICJodHRwczovL3N0YXRpYy52ZWN0ZWV6eS5jb20vc3lzdGVtL3Jlc291cmNlcy9wcmV2aWV3cy8wMDEvMTg5LzA2My9ub25fMngvc3Rhci1yb3VuZGVkLXBuZy5wbmciLA0KICBpY29uV2lkdGggPSAxMCwgaWNvbkhlaWdodCA9IDEwDQogICNpY29uQW5jaG9yWCA9IDIwLCBpY29uQW5jaG9yWSA9IDIwDQogIA0KKQ0KDQpsYWJlbF93b3JsZCA8LSBsYWJlbF93b3JsZCAlPiUgZmlsdGVyKHJhbmtpbmcgPiAyMCkgDQogIA0KbSA8LSBsZWFmbGV0KHdpZHRoPTEyMDAsIGhlaWdodD04MDApICU+JSBhZGRUaWxlcygpICANCm0gJTw+JSAgYWRkQ2lyY2xlTWFya2VycyhsYWJlbF93b3JsZCRMb25nLCBsYWJlbF93b3JsZCRMYXQsDQogICAgICAgICAgICAgICAgICAgICAgICAjIHJhZGl1cz0yK2xvZzIoeCRjb25maXJtZWQpLA0KICAgICAgICAgICAgICAgICAgICAgICAgcmFkaXVzPTEwLCMqbG9nMihtLndvcmxkJGF2Zy50ZW1wKSwNCiAgICAgICAgICAgICAgICAgICAgICAgIHN0cm9rZT1GLA0KICAgICAgICAgICAgICAgICAgICAgICAgI2NvbG9yPSdyZWQnLA0KICAgICAgICAgICAgICAgICAgICAgICAgY29sb3IgPSB3cGFsKGxhYmVsX3dvcmxkJGF2Z190ZW1wKSwgDQogICAgICAgICAgICAgICAgICAgICAgICBmaWxsT3BhY2l0eT0wLjUsDQogICAgICAgICAgICAgICAgICAgICAgICAjcG9wdXA9bGFiZWwudG9wJHR4dA0KICAgICAgICAgICAgICAgICAgICAgICAgbGFiZWw9IGxhYmVsX3dvcmxkJHR4dCwNCiAgICAgICAgICAgICAgICAgICAgICAgIGdyb3VwID0gIldvcmxkIg0KICAgICAgICAgICAgICAgICAgICAgICAgKSAlPiUNCiAgDQogIGFkZENpcmNsZU1hcmtlcnMobGFiZWxfdG9wJExvbmcsIGxhYmVsX3RvcCRMYXQsDQogICAgICAgICAgICAgICAgICAgICAgICAjIHJhZGl1cz0yK2xvZzIoeCRjb25maXJtZWQpLA0KICAgICAgICAgICAgICAgICAgICAgICAgcmFkaXVzPTEwLCMqbG9nMihtLndvcmxkJGF2Zy50ZW1wKSwNCiAgICAgICAgICAgICAgICAgICAgICAgIHN0cm9rZT1GLA0KICAgICAgICAgICAgICAgICAgICAgICAgI2NvbG9yPSdyZWQnLA0KICAgICAgICAgICAgICAgICAgICAgICAgY29sb3IgPSB3cGFsKGxhYmVsX3RvcCRhdmdfdGVtcCksIA0KICAgICAgICAgICAgICAgICAgICAgICAgZmlsbE9wYWNpdHk9MC41LA0KICAgICAgICAgICAgICAgICAgICAgICAgI3BvcHVwPWxhYmVsLnRvcCR0eHQNCiAgICAgICAgICAgICAgICAgICAgICAgIGxhYmVsPSBsYWJlbF90b3AkdHh0LA0KICAgICAgICAgICAgICAgICAgICAgICAgZ3JvdXAgPSAiVG9wIDIwIENvdW50cmllcyINCiAgICAgICAgICAgICAgICAgICAgICAgICkgJT4lDQogIA0KICBhZGRMYWJlbE9ubHlNYXJrZXJzKGxhYmVsX3RvcCRMb25nLCBsYWJlbF90b3AkTGF0LCBsYWJlbCA9IGxhYmVsX3RvcCRyYW5raW5nLA0KICAgICAgICAgICAgICAgICAgICAgIGxhYmVsT3B0aW9ucyA9IGxhYmVsT3B0aW9ucyhub0hpZGUgPSBUUlVFLCB0ZXh0T25seSA9IFRSVUUsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkaXJlY3Rpb24gPSAiaGVhZCIsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvZmZzZXQgPSBjKDUsNCkpLA0KICAgICAgICAgICAgICAgICAgICAgIGdyb3VwID0gIlRvcCAyMCBDb3VudHJpZXMiKSAlPiUNCiAgDQogIGFkZExlZ2VuZCgiYm90dG9tcmlnaHQiLCBwYWwgPSB3cGFsLCB2YWx1ZXMgPSBsYWJlbF93b3JsZCRhdmdfdGVtcCwgb3BhY2l0eSA9IDEsDQogICAgICAgICAgICBsYWJGb3JtYXQgPSBsYWJlbEZvcm1hdChzdWZmaXggPSAiIMKwQyIpLA0KICAgICAgICAgICAgdGl0bGUgPSAiVGVtcGVyYXR1cmUiKSAlPiUgDQogIA0KICBhZGRMYXllcnNDb250cm9sKA0KICAgICNiYXNlR3JvdXBzID0gYygiT1NNIChkZWZhdWx0KSIsICJUb25lciIsICJUb25lciBMaXRlIiksDQogICAgb3ZlcmxheUdyb3VwcyA9IGMoIlRvcCAyMCBDb3VudHJpZXMiLCAiV29ybGQiKSwNCiAgICBvcHRpb25zID0gbGF5ZXJzQ29udHJvbE9wdGlvbnMoY29sbGFwc2VkID0gRkFMU0UpDQogICkNCm0NCmBgYA0KDQpgYGB7cn0NCmdkcC50b3AyMCA8LSBkZl9nZHAyMDE5ICU+JQ0KICBzZWxlY3QoYygicmFuayIsICJjb3VudHJ5IiwgIkdEUCAobWlsbGlvbnMgb2YgVVMgZG9sbGFycykiKSkgJT4lDQogIG1lcmdlKGRhdGEubGF0ZXN0LmFsbCAlPiUgDQogICAgICAgICAgc2VsZWN0KGNvdW50cnksIHJhbmtpbmcsIGNvbmZpcm1lZCwgcmVjb3ZlcmVkLCBkZWF0aHMpICU+JSANCiAgICAgICAgICBmaWx0ZXIoY291bnRyeSAlaW4lIHRvcC5jb3VudHJpZXMgJiBjb3VudHJ5ICE9ICJXb3JsZCIpLCBieSA9ICJjb3VudHJ5IikgJT4lDQogIGFycmFuZ2UocmFua2luZykgJT4lIA0KICBtdXRhdGUocmFua2luZyA9IHJhbmtpbmcgLSAxKSANCmdkcC50b3AyMCAlPD4lIHJlbmFtZSgiR0RQIiA9ICJHRFAgKG1pbGxpb25zIG9mIFVTIGRvbGxhcnMpIikNCmdkcC50b3AyMA0KDQpnIDwtIGdncGxvdChnZHAudG9wMjAsIGFlcyh4ID0gR0RQLCB5ID0gcmVvcmRlcihjb3VudHJ5LCAtcmFua2luZykpKSArDQogIGdlb21faGlzdG9ncmFtKHN0YXQgPSAiaWRlbnRpdHkiLCBhZXMoZmlsbCA9IEdEUCkpKyANCiAgc2NhbGVfZmlsbF9ncmFkaWVudCgiR0RQIiwgbG93ID0gIiNGRjQwMzgiLCBoaWdoID0gIiM1MEU5NTIiKSArIA0KICBsYWJzKHRpdGxlPXBhc3RlMCgiR0RQICBvZiBUb3AgMjAgQ291bnRyaWVzIGluIDIwMTkgKG1pbGxpb25zIG9mIFVTIGRvbGxhcnMpIikpICsNCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbD1HRFAsIHggPSBHRFApLCBzaXplPTMuNSwgaGp1c3Q9LTAuMikgKw0KICB4bGFiKCJHRFAgKG1pbGxpb25zIG9mIFVTIGRvbGxhcnMpIikgKw0KICB5bGFiKCIiKSArDQogIHRoZW1lKGxlZ2VuZC50aXRsZT1lbGVtZW50X2JsYW5rKCkpDQpnDQoNCmdseS50b3AuZ2RwIDwtIGdncGxvdGx5KGcpDQpnbHkudG9wLmdkcA0KYGBgDQoNCmBgYHtyfQ0KIyBQb3JuaHViDQpkZl9wb3JuaHViIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIHBvcm5odWIiKSkNCmRmX3Bvcm5odWIgPC0gYXMuZGF0YS5mcmFtZShkZl9wb3JuaHViKQ0KZGZfcG9ybmh1Yg0KYGBgDQoNCg0KYGBge3J9DQojIFNhcnMNCmRmX3NhcnMgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gc2Fyc18yMDAzIikpDQpkZl9zYXJzIDwtIGFzLmRhdGEuZnJhbWUoZGZfc2FycykNCiNWaWV3KGRmX3NhcnMpDQoNCmRhdGVzLnMgPC0gZGZfc2Fyc1ssMV0lPiUgbWR5KCkNCnJhbmdlKGRhdGVzLnMpDQptaW4uZGF0ZS5zIDwtIG1pbihkYXRlcy5zKQ0KbWF4LmRhdGUucyA8LSBtYXgoZGF0ZXMucykNCm1pbi5kYXRlLnR4dC5zIDwtIG1pbi5kYXRlLnMgJT4lIGZvcm1hdCgnJWQgJWIgJVknKQ0KbWF4LmRhdGUudHh0LnMgPC0gbWF4LmRhdGUucyAlPiUgZm9ybWF0KCclZCAlYiAlWScpDQpkYXkxLnNhcnMgPC0gbWluLmRhdGUucw0KYGBgDQoNCmBgYHtyfQ0KIyBjbGVhbiBkYXRhIFNhcnMgDQpkYXRhLnNhcnMgPC0gZGZfc2FycyAlPiUgcmVuYW1lKGMoImRhdGUiID0gIkRhdGUiLCAiY29uZmlybWVkIiA9ICJDdW11bGF0aXZlX251bWJlciIgICwiZGVhdGhzIiA9ICJOdW1iZXJfZGVhdGhzIiwgInJlY292ZXJlZCIgPSAiTnVtYmVyX3JlY292ZXJlZCIpKSAlPiUNCiAgbXV0YXRlKGRhdGUgPSBkYXRlICU+JSBtZHkoKSkgJT4lDQogIGdyb3VwX2J5KGNvdW50cnksIGRhdGUpICU+JSBhcy5kYXRhLmZyYW1lKCkgDQoNCiMgQWRkIFdvcmxkJ3MgU2FycyBjYXNlcyANCndvcmxkLnNhcnMgPC0gZGF0YS5zYXJzICU+JSBncm91cF9ieShkYXRlKSAlPiUgDQogIHN1bW1hcmlzZShjb3VudHJ5PSdXb3JsZCcsDQogICAgICAgICAgICBjb25maXJtZWQgPSBzdW0oY29uZmlybWVkLCBuYS5ybT1UKSwNCiAgICAgICAgICAgIGRlYXRocyA9IHN1bShkZWF0aHMsIG5hLnJtPVQpLA0KICAgICAgICAgICAgcmVjb3ZlcmVkID0gc3VtKHJlY292ZXJlZCwgbmEucm09VCkpDQoNCmRhdGEuc2FycyAlPD4lIHJiaW5kKHdvcmxkLnNhcnMpDQpkYXRhLnNhcnMgJTw+JSBtdXRhdGUoY3VycmVudC5jb25maXJtZWQgPSBjb25maXJtZWQgLSBkZWF0aHMgLSByZWNvdmVyZWQpDQojVmlldyh3b3JsZC5zYXJzKSANCmRhdGEuc2Fycw0KYGBgDQoNCg0KYGBge3J9DQojcmF0ZQ0KZGF0YS5zYXJzICU8PiUgYXJyYW5nZShjb3VudHJ5LCBkYXRlKQ0KbiA8LSBucm93KGRhdGEuc2FycykNCmRheTEuc2FycyA8LSBtaW4oZGF0YS5zYXJzJGRhdGUpDQpkYXRhLnNhcnMgJTw+JSBtdXRhdGUobmV3LmNvbmZpcm1lZCA9IGlmZWxzZShkYXRlID09IGRheTEuc2FycywgTkEsIGNvbmZpcm1lZCAtIGxhZyhjb25maXJtZWQsIG49MSkpLA0KICAgICAgICAgICAgICAgICBuZXcuZGVhdGhzID0gaWZlbHNlKGRhdGUgPT0gZGF5MS5zYXJzLCBOQSwgZGVhdGhzIC0gbGFnKGRlYXRocywgbj0xKSksDQogICAgICAgICAgICAgICAgIG5ldy5yZWNvdmVyZWQgPSBpZmVsc2UoZGF0ZSA9PSBkYXkxLnNhcnMsIE5BLCByZWNvdmVyZWQgLSBsYWcocmVjb3ZlcmVkLCBuPTEpKSkNCmRhdGEuc2FycyAlPD4lIG11dGF0ZShuZXcuY29uZmlybWVkID0gaWZlbHNlKG5ldy5jb25maXJtZWQgPCAwLCAwLCBuZXcuY29uZmlybWVkKSwNCiAgICAgICAgICAgICAgICAgbmV3LmRlYXRocyA9IGlmZWxzZShuZXcuZGVhdGhzIDwgMCwgMCwgbmV3LmRlYXRocyksDQogICAgICAgICAgICAgICAgIG5ldy5yZWNvdmVyZWQgPSBpZmVsc2UobmV3LnJlY292ZXJlZCA8IDAsIDAsIG5ldy5yZWNvdmVyZWQpKQ0KIyMgZGVhdGggcmF0ZSBiYXNlZCBvbiB0b3RhbCBkZWF0aHMgYW5kIHJlY292ZXJlZCBjYXNlcw0KZGF0YS5zYXJzICU8PiUgbXV0YXRlKHJhdGUudXBwZXIgPSAoMTAwICogZGVhdGhzIC8gKGRlYXRocyArIHJlY292ZXJlZCkpICU+JSByb3VuZCgxKSwNCiAgICAgICAgICAgICAgICAgcmF0ZS51cHBlciA9IGlmZWxzZShpcy5uYW4ocmF0ZS51cHBlciksIDAsIHJhdGUudXBwZXIpKQ0KDQojIyBsb3dlciBib3VuZDogZGVhdGggcmF0ZSBiYXNlZCBvbiB0b3RhbCBjb25maXJtZWQgY2FzZXMNCmRhdGEuc2FycyAlPD4lIG11dGF0ZShyYXRlLmxvd2VyID0gKDEwMCAqIGRlYXRocyAvIGNvbmZpcm1lZCkgJT4lIHJvdW5kKDEpLA0KICAgICAgICAgICAgICAgICByYXRlLmxvd2VyID0gaWZlbHNlKGlzLm5hbihyYXRlLmxvd2VyKSwgMCwgcmF0ZS5sb3dlcikpDQoNCiMjIGRlYXRoIHJhdGUgYmFzZWQgb24gdGhlIG51bWJlciBvZiBkZWF0aC9yZWNvdmVyZWQgb24gZXZlcnkgc2luZ2xlIGRheQ0KZGF0YS5zYXJzICU8PiUgbXV0YXRlKHJhdGUuZGFpbHkgPSAoMTAwICogbmV3LmRlYXRocyAvIChuZXcuZGVhdGhzICsgbmV3LnJlY292ZXJlZCkpICU+JSByb3VuZCgxKSwNCiAgICAgICAgICAgICAgICAgcmF0ZS5kYWlseSA9IGlmZWxzZShpcy5uYW4ocmF0ZS5kYWlseSksIDAsIHJhdGUuZGFpbHkpKQ0KDQojVmlldyhkYXRhLnNhcnMpDQpgYGANCg0KYGBge3J9DQojIyBjb252ZXJ0IGZyb20gd2lkZSB0byBsb25nIGZvcm1hdA0KZGF0YS5zYXJzLmxvbmcgPC0gZGF0YS5zYXJzICU+JQ0KICBzZWxlY3QoYyhjb3VudHJ5LCBkYXRlLCBjb25maXJtZWQsIGN1cnJlbnQuY29uZmlybWVkLCByZWNvdmVyZWQsIGRlYXRocykpICU+JQ0KICBnYXRoZXIoa2V5PXR5cGUsIHZhbHVlPWNvdW50LCAtYyhjb3VudHJ5LCBkYXRlKSkNCiMjIHNldCBmYWN0b3IgbGV2ZWxzIHRvIHNob3cgdGhlbSBpbiBhIGRlc2lyYWJsZSBvcmRlcg0KZGF0YS5zYXJzLmxvbmcgJTw+JSBtdXRhdGUodHlwZT1yZWNvZGVfZmFjdG9yKHR5cGUsIGNvbmZpcm1lZD0nVG90YWwgQ29uZmlybWVkJywNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY3VycmVudC5jb25maXJtZWQ9J0N1cnJlbnQgQ29uZmlybWVkJywNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcmVjb3ZlcmVkPSdSZWNvdmVyZWQnLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkZWF0aHM9J0RlYXRocycpKQ0KI1ZpZXcoZGF0YS5zYXJzLmxvbmcpDQoNCiMgV29ybGQgc2FycycgbG9uZyBkYXRhIA0Kd29ybGQuc2Fycy5sb25nIDwtIGRhdGEuc2Fycy5sb25nICU+JQ0KICBmaWx0ZXIoY291bnRyeSA9PSAiV29ybGQiKQ0KI1ZpZXcod29ybGQuc2Fycy5sb25nKQ0KDQpnIDwtIGdncGxvdCh3b3JsZC5zYXJzLmxvbmcsIGFlcyhkYXRlLCBjb3VudCwgY29sb3IgPSB0eXBlKSkgKw0KICBnZW9tX2xpbmUoKSsNCiAgbGFicyh0aXRsZSA9ICJOdW1iZXIgb2YgQ2FzZXMgV29ybGR3aWRlOiBTQVJzIikrDQogIHhsYWIoIiIpKw0KICB5bGFiKCIiKQ0KZw0KDQpnbHkuZyA8LSBnZ3Bsb3RseShnKQ0KZ2x5LmcNCmdseS5wbG90X3dvcmxkLmxpbmUNCmBgYA0KDQpgYGB7cn0NCmRmX3NhcnNfbGFzdGRhdGUgPC0gZGF0YS5zYXJzICU+JQ0KICBmaWx0ZXIoZGF0ZSA9PSBtYXguZGF0ZS5zKQ0KDQpkZl9zYXJzX2xhc3RkYXRlDQpgYGANCg0KDQpgYGB7cn0NCiMjIEN1cnJlbnQgQ29uZmlybWVkIENhc2VzDQpkYXRhLnNhcnMud29ybGQgPC0gZGF0YS5zYXJzICU+JSBmaWx0ZXIoY291bnRyeT09J1dvcmxkJykNCiNWaWV3KGRhdGEuc2Fycy53b3JsZCkNCm4gPC0gbnJvdyhkYXRhLnNhcnMud29ybGQpDQojVmlldyhkYXRhLnNhcnMud29ybGQpDQpwbG90X3NhcnMuY3VycmNvbmYgPC0gZ2dwbG90KGRhdGEuc2Fycy53b3JsZCwgYWVzKHg9ZGF0ZSwgeT1jdXJyZW50LmNvbmZpcm1lZCkpICsNCiAgZ2VvbV9wb2ludCgpICsgZ2VvbV9zbW9vdGgoKSArDQogIHhsYWIoJycpICsgeWxhYignQ291bnQnKSArIGxhYnModGl0bGU9J0N1cnJlbnQgQ29uZmlybWVkIENhc2VzJykgKw0KICB0aGVtZShheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKQ0KcGxvdF9zYXJzLm5ld2NvbmYgPC0gZ2dwbG90KGRhdGEuc2Fycy53b3JsZCwgYWVzKHg9ZGF0ZSwgeT1uZXcuY29uZmlybWVkKSkgKw0KICBnZW9tX3BvaW50KCkgKyBnZW9tX3Ntb290aCgpICsNCiAgeGxhYignJykgKyB5bGFiKCdDb3VudCcpICsgbGFicyh0aXRsZT0nRGFpbHkgTmV3IENvbmZpcm1lZCBDYXNlcycpICsNCiAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdD0xKSkNCiMjIHNob3cgdHdvIHBsb3RzIHNpZGUgYnkgc2lkZQ0KZ3JpZC5hcnJhbmdlKHBsb3Rfc2Fycy5jdXJyY29uZiwgcGxvdF9zYXJzLm5ld2NvbmYsIG5jb2w9MikNCmBgYA0KDQpgYGB7cn0NCiMjIGEgc2NhdHRlciBwbG90IHdpdGggYSBzbW9vdGhlZCBsaW5lIGFuZCB2ZXJ0aWNhbCB4LWF4aXMgbGFiZWxzDQpwbG90X3NhcnMuYWNjZGVhdGhzIDwtIGdncGxvdChkYXRhLnNhcnMud29ybGQsIGFlcyh4PWRhdGUsIHk9ZGVhdGhzKSkgKw0KICBnZW9tX3BvaW50KCkgKyBnZW9tX3Ntb290aCgpICsNCiAgeGxhYignJykgKyB5bGFiKCdDb3VudCcpICsgbGFicyh0aXRsZT0nQWNjdW11bGF0aXZlIERlYXRocycpICsNCiAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdD0xKSkNCnBsb3Rfc2Fycy5hY2NyZWNvdiA8LSBnZ3Bsb3QoZGF0YS5zYXJzLndvcmxkLCBhZXMoeD1kYXRlLCB5PXJlY292ZXJlZCkpICsNCiAgZ2VvbV9wb2ludCgpICsgZ2VvbV9zbW9vdGgoKSArDQogIHhsYWIoJycpICsgeWxhYignQ291bnQnKSArIGxhYnModGl0bGU9J0FjY3VtdWxhdGl2ZSBSZWNvdmVyZWQgQ2FzZXMnKSArDQogIHRoZW1lKGF4aXMudGV4dC54PWVsZW1lbnRfdGV4dChhbmdsZT00NSwgaGp1c3Q9MSkpDQpwbG90X3NhcnMubmV3ZGVhdGhzIDwtIGdncGxvdChkYXRhLnNhcnMud29ybGQsIGFlcyh4PWRhdGUsIHk9bmV3LmRlYXRocykpICsNCiAgZ2VvbV9wb2ludCgpICsgZ2VvbV9zbW9vdGgoKSArDQogIHhsYWIoJycpICsgeWxhYignQ291bnQnKSArIGxhYnModGl0bGU9J05ldyBEZWF0aHMnKSArDQogIHRoZW1lKGF4aXMudGV4dC54PWVsZW1lbnRfdGV4dChhbmdsZT00NSwgaGp1c3Q9MSkpDQpwbG90X3NhcnMubmV3cmVjb3YgPC0gZ2dwbG90KGRhdGEuc2Fycy53b3JsZCwgYWVzKHg9ZGF0ZSwgeT1uZXcucmVjb3ZlcmVkKSkgKw0KICBnZW9tX3BvaW50KCkgKyBnZW9tX3Ntb290aCgpICsNCiAgeGxhYignJykgKyB5bGFiKCdDb3VudCcpICsgbGFicyh0aXRsZT0nTmV3IFJlY292ZXJlZCBDYXNlcycpICsNCiAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdD0xKSkNCiMjIHNob3cgZm91ciBwbG90cyB0b2dldGhlciwgd2l0aCAyIHBsb3RzIGluIGVhY2ggcm93DQpncmlkLmFycmFuZ2UocGxvdF9zYXJzLmFjY2RlYXRocywgcGxvdF9zYXJzLmFjY3JlY292LCBwbG90X3NhcnMubmV3ZGVhdGhzLCBwbG90X3NhcnMubmV3cmVjb3YsIG5yb3c9MikNCmBgYA0KDQoNCmBgYHtyfQ0KI1RvcCAyMCB3aXRoIGdkcA0KZGF0YS5sb25nR0RQIDwtIGRmX2dkcDIwMTkgJT4lIGdhdGhlcihrZXk9eWVhciwgdmFsdWU9R0RQLCAtYyhjb3VudHJ5KSkNCmRhdGEudG9wIDwtIGRhdGEubGF0ZXN0ICU+JSBmaWx0ZXIoIWNvdW50cnkgJWluJSBjKCdXb3JsZCcsICdPdGhlcnMnKSkNCmRhdGEudG9wIDwtIGhlYWQoZGF0YS50b3AsMjApDQojVmlldyhkYXRhLmxhdGVzdCkNCmRhdGEuZ2RwIDwtIGZpbHRlcihkYXRhLmxvbmdHRFAseWVhcj09JzIwMjAnKQ0KZGZfc2Fyc19sYXN0ZGF0ZV9jb25maXJtZWQgPC0gZGZfc2Fyc19sYXN0ZGF0ZSAlPiUNCiAgc2VsZWN0KCJjb3VudHJ5IiwgImNvbmZpcm1lZCIpICU+JQ0KICByZW5hbWUoc2FycyA9ICJjb25maXJtZWQiKQ0KI1ZpZXcoZGZfc2Fyc19sYXN0ZGF0ZV9jb25maXJtZWQpDQojbWVyZ2UNCm1lcmdjb3VudHJ5ID0gZnVuY3Rpb24oZGF0YTEsZGF0YTIpew0KICBkYXRhIDwtIG1lcmdlKHggPSBkYXRhMSwgeSA9IGRhdGEyLCBieSA9ICJjb3VudHJ5IiwgYWxsLnggPSBUUlVFKSANCiAgcmV0dXJuKGRhdGEpDQp9DQpkYXRhLnRvcC53b3JsZCA8LSBtZXJnZSh4ID0gZGF0YS50b3AsIHkgPSBkZl9nZHAyMDE5LCBieSA9ICJjb3VudHJ5IiwgYWxsLnggPSBUUlVFKSAlPiUgDQogIHNlbGVjdCgtYyhjb2RlLHJhbmssbmV3LmNvbmZpcm1lZCxuZXcuZGVhdGhzLGN1cnJlbnQuY29uZmlybWVkLHBvcHVsYXRpb24pKSAlPiUgDQogIHJlbmFtZShHRFA9IkdEUCAobWlsbGlvbnMgb2YgVVMgZG9sbGFycykiKQ0KDQpkYXRhLnRvcC53b3JsZCA8LSBtZXJnZSh4ID0gZGF0YS50b3Aud29ybGQsIHkgPSBkZl9oZWFsdCwgYnkgPSAiY291bnRyeSIsIGFsbC54ID0gVFJVRSkgJT4lDQogIHJlbmFtZShoZWFsdGhjYXJlPSJoZWFsdGhDYXJlSW5kZXgiKQ0KI2RhdGEudG9wLndvcmxkIDwtIG1lcmdjb3VudHJ5KGRhdGEudG9wLndvcmxkLCBkZl90ZW1wKQ0KDQpkYXRhLnRvcC53b3JsZCA8LSBtZXJnZSh4ID0gZGF0YS50b3Aud29ybGQsIHkgPSBkZl9wb3JuaHViLCBieSA9ICJjb3VudHJ5IiwgYWxsLnggPSBUUlVFKSAlPiUNCiAgcmVuYW1lKFBvcm5odWIgPSAiUG9ybmh1YkluZGV4KCUpIikNCg0KZGF0YS50b3Aud29ybGQgPC0gbWVyZ2UoeCA9IGRhdGEudG9wLndvcmxkLCB5ID0gZGZfc2Fyc19sYXN0ZGF0ZV9jb25maXJtZWQsIGJ5ID0gImNvdW50cnkiLCBhbGwueCA9IFRSVUUpIA0KDQoNCmRhdGEudG9wLndvcmxkIDwtIG1lcmdjb3VudHJ5KGRhdGEudG9wLndvcmxkLCBkZl90ZW1wKQ0KaW5kZXggPC0gaXMubmEoZGF0YS50b3Aud29ybGQpDQpkYXRhLnRvcC53b3JsZFtpbmRleF0gPC0gMA0KZGF0YS50b3Aud29ybGQNClZpZXcoZGF0YS50b3Aud29ybGQpDQoNCg0Kbm9ybWFsaXplID0gZnVuY3Rpb24oZGF0YSl7DQogICNyZXR1cm4gKChkYXRhIC0gbWluKGRhdGEsbmEucm0gPSBUUlVFKSkvKG1heChkYXRhLG5hLnJtID0gVFJVRSkgLSBtaW4oZGF0YSxuYS5ybSA9IFRSVUUpKSkNCiAgeiA8LSBzY2FsZShkYXRhKTsNCiAgdGFuaCh6LzIpDQp9DQoNCm5vcm1fZGF0YSA9IGFzLmRhdGEuZnJhbWUoYXBwbHkoZGF0YS50b3Aud29ybGRbLDI6MTNdLDIsbm9ybWFsaXplKSkNCmNvcnJfZGF0YSA8LSBub3JtX2RhdGENCm5vcm1fZGF0YSRjb3VudHJ5IDwtIGMoIkFyZ2VudGluYSIsIkJhbmdsYWRlc2giLCJCcmF6aWwiLCJDaGlsZSIsIkNvbG9tYmlhIiwiRnJhbmNlIiwiR2VybWFueSIsIkluZGlhIiwiSXJhbiIsIkl0YWx5IiwiTWV4aWNvIiwiUGFraXN0YW4iLCJQZXJ1IiwiUnVzc2lhIiwic2F1ZGkgQXJhYmlhIiwiU291dGggQWZyaWNhIiwiU3BhaW4iLCJUdXJrZXkiLCJVbml0ZWQgS2luZ2RvbSIsIlVTIikNCm5vcm1fZGF0YSA8LSBub3JtX2RhdGFbLGMobmNvbChub3JtX2RhdGEpLDE6KG5jb2wobm9ybV9kYXRhKS0xKSldDQpyb3duYW1lcyhub3JtX2RhdGEpIDwtIG5vcm1fZGF0YVssMV0NCm5vcm1fZGF0YVssMV0gPC0gTlVMTA0KDQojaGVhdG1hcChhcy5tYXRyaXgobm9ybV9kYXRhWywgLTFdKSxzY2FsZT0iY29sdW1uIixjb2w9IGNvbG9yUmFtcFBhbGV0dGUoYnJld2VyLnBhbCg4LCAiQmx1ZXMiKSkoMjUpKQ0KI2xlZ2VuZCh4PSJib3R0b21yaWdodCIsIGxlZ2VuZD1jKCItMSIsICIwIiwgIjEiKSwgDQojICAgICBmaWxsPWNvbG9yUmFtcFBhbGV0dGUoYnJld2VyLnBhbCg4LCAiQmx1ZXMiKSkoMykpDQoNCiMgaW5zdGFsbC5wYWNrYWdlcygiZ3Bsb3RzIikNCmxpYnJhcnkoImdwbG90cyIpDQpub3JtX2RhdGFfcGxvdCA8LSBzZWxlY3Qobm9ybV9kYXRhLCJjb25maXJtLnJhdGUiLCJkZWF0aC5yYXRlIiwicmVjb3Zlci5yYXRlIiwiaGVhbHRoY2FyZSIsIlBvcm5odWIiLCJHRFAiLCJhdmdfdGVtcCIsICJzYXJzIikNCm5vcm1fZGF0YV9wbG90DQpoZWF0bWFwLjIoYXMubWF0cml4KG5vcm1fZGF0YV9wbG90KSwNCiAgICAgICAgICAgICAgICAgICAgc2NhbGU9Im5vbmUiLCANCiAgICAgICAgICAgICAgICAgICAgY29sID0gY29sb3JSYW1wUGFsZXR0ZShjKCIjNkQ5RUMxIiwid2hpdGUiLCIjRTQ2NzI2IikpKG4gPSAyMDApLA0KICAgICAgICAgICAgICAgICAgICBtYXJnaW5zPWMoMTAsNiksdHJhY2U9ImNvbHVtbiIpDQoNCiNub3JtX2RhdGFfcGxvdCA8LSBzZWxlY3Qobm9ybV9kYXRhLCJjb3VudHJ5IiwiY29uZmlybS5yYXRlIiwiZGVhdGgucmF0ZSIsInJlY292ZXIucmF0ZSIsImhlYWx0aGNhcmUiLCJQb3JuaHViIiwiR0RQIiwiYXZnX3RlbXAiLCAic2FycyIpDQojbm9ybV9kYXRhX3Bsb3QgJTw+JSBnYXRoZXIoa2V5PXR5cGUsIHZhbHVlPWNvdW50LCAtYyhjb3VudHJ5KSkNCiNsZXZlbF9vcmRlciA8LSBmYWN0b3Iobm9ybV9kYXRhX3Bsb3QkdHlwZSwgDQogICAgICAgICAgICAgICAgICAgICAgbGV2ZWwgPSAjYygic2FycyIsIlBvcm5odWIiLCJHRFAiLCJhdmdfdGVtcCIsImhlYWx0aGNhcmUiLCJyZWNvdmVyLnJhdGUiLCJkZWF0aC5yYXRlIiwiY29uZmlybS5yYXRlIikpDQpWaWV3KGxldmVsX29yZGVyKQ0KZ2dwbG90KGRhdGEgPSBub3JtX2RhdGFfcGxvdCwgYWVzKHg9Y291bnRyeSwgeT1sZXZlbF9vcmRlciwgZmlsbD1jb3VudCkpICsgDQogIGdlb21fdGlsZSgpICsNCiAgc2NhbGVfZmlsbF9ncmFkaWVudChsb3cgPSAicGluayIsIGhpZ2ggPSAiYmx1ZSIpICsNCiAgeGxhYigiIikgKw0KICB5bGFiKCIiKSArDQogIHRoZW1lX2J3KCkgKw0KICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLHZqdXN0ID0gMSkpKw0KICB0aGVtZSgNCiAgICBheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksDQogICAgYXhpcy50aWNrcyA9IGVsZW1lbnRfYmxhbmsoKSwNCiAgICBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpLA0KICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksDQogICAgcGFuZWwuYm9yZGVyID0gZWxlbWVudF9ibGFuaygpLA0KICAgIHBhbmVsLmJhY2tncm91bmQgPSBlbGVtZW50X2JsYW5rKCksDQogICAgI2xlZ2VuZC5wb3NpdGlvbiA9ICJub25lIg0KICApDQoNCmBgYA0KDQpgYGB7cn0NCiNjb3JyZWxhdGlvbg0KY29ycl9kYXRhICU8PiUgc2VsZWN0KGMoR0RQLGNvbmZpcm0ucmF0ZSxkZWF0aC5yYXRlLHJlY292ZXIucmF0ZSxoZWFsdGhjYXJlLGF2Z190ZW1wLFBvcm5odWIsIHNhcnMpKQ0KaGVhZChjb3JyX2RhdGEpDQpjb3IoY29ycl9kYXRhKQ0KcGxvdGNvcnIgPC0gZ2djb3JycGxvdChjb3IoY29ycl9kYXRhKSxoYy5vcmRlciA9IFRSVUUsDQogICAgICAgICAgIG91dGxpbmUuY29sb3IgPSAid2hpdGUiLA0KICAgICAgICAgICBjb2xvcnMgPSBjKCIjNkQ5RUMxIiwid2hpdGUiLCIjRTQ2NzI2IiksDQogICAgICAgICAgIGxhYiA9IFRSVUUpDQoNCnBsb3Rjb3JyDQpgYGANCg0KYGBge3J9DQojIyBEYXRhIGluIFVTDQoNCmRmX3VzIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIGNvdmlkdXMiKSkNCmRmX3VzIDwtIGFzLmRhdGEuZnJhbWUoZGZfdXMpIA0KZGZfdXMgPC0gc2VsZWN0KGRmX3VzLCBkYXRlLCBzdGF0ZSwgY2FzZXMsIGRlYXRocykNCg0KZGZfdXNfcG9wIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIGRhdGFfcG9wdWxhdGlvbiIpKSANCmRmX3VzX3BvcCA8LSBhcy5kYXRhLmZyYW1lKGRmX3VzX3BvcCkNCg0KZGZfdXNfZ2VuZGVyIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIGRhdGFfZ2VuZGVyIikpIA0KZGZfdXNfZ2VuZGVyIDwtIGFzLmRhdGEuZnJhbWUoZGZfdXNfZ2VuZGVyKQ0KDQpkZl91c19ldGhuaWMgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gZGF0YV9ldGhuaWMiKSkgDQpkZl91c19ldGhuaWMgPC0gYXMuZGF0YS5mcmFtZShkZl91c19ldGhuaWMpDQoNCmRmX3VzX2xvY2tkb3duIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIGRhdGFfbG9ja2Rvd24iKSkgDQpkZl91c19sb2NrZG93biA8LSBhcy5kYXRhLmZyYW1lKGRmX3VzX2xvY2tkb3duKQ0KDQpkZl91c19oZWFsdGggPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gZGF0YV9oZWFsdGgiKSkgDQpkZl91c19oZWFsdGggPC0gYXMuZGF0YS5mcmFtZShkZl91c19oZWFsdGgpDQoNCmRmX3VzX3Rlc3RpbmcgPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gZGF0YV90ZXN0aW5nIikpIA0KZGZfdXNfdGVzdGluZyA8LSBhcy5kYXRhLmZyYW1lKGRmX3VzX3Rlc3RpbmcpDQoNCmBgYA0KDQpgYGB7cn0NCmRmX3VzIDwtIGRmX3VzICAlPiUNCiAgbXV0YXRlKGRhdGUgPSBkYXRlICU+JSBtZHkoKSkgJT4lDQogIHJlbmFtZSgiY29uZmlybWVkIiA9ICJjYXNlcyIpDQojZGF0YS51cw0KZGF0ZXMudXMgPC0gZGZfdXNbLDFdDQpyYW5nZShkYXRlcy51cykNCm1pbi5kYXRlLnVzIDwtIG1pbihkYXRlcy51cykNCm1heC5kYXRlLnVzIDwtIG1heChkYXRlcy51cykNCm1pbi5kYXRlLnR4dC51cyA8LSBtaW4uZGF0ZS51cyAlPiUgZm9ybWF0KCclZCAlYiAlWScpDQptYXguZGF0ZS50eHQudXMgPC0gbWF4LmRhdGUudXMgJT4lIGZvcm1hdCgnJWQgJWIgJVknKQ0KZGF5MS51cyA8LSBtaW4oZGZfdXMkZGF0ZSkNCg0KZGF0YS51cy50b3RhbCA8LSBkZl91cyAlPiUgZ3JvdXBfYnkoZGF0ZSkgJT4lDQogIHN1bW1hcmlzZShzdGF0ZT0nVVMnLA0KICAgICAgICAgICAgY29uZmlybWVkID0gc3VtKGNvbmZpcm1lZCwgbmEucm09VCksDQogICAgICAgICAgICBkZWF0aHMgPSBzdW0oZGVhdGhzLCBuYS5ybT1UKSkNCiNWaWV3KGRhdGEudXMudG90YWwpDQoNCmRhdGEudXMgPC0gZGZfdXMNCmRhdGEudXMgJTw+JSByYmluZChkYXRhLnVzLnRvdGFsKQ0KVmlldyhkYXRhLnVzKQ0KDQpkYXRhLnVzLmxvbmcgPC0gZGF0YS51cyAlPiUgDQogIGdhdGhlcihrZXkgPSB0eXBlLCB2YWx1ZSA9IGNvdW50LCAtYyhkYXRlLCBzdGF0ZSkpIA0KVmlldyhkYXRhLnVzLmxvbmcpDQoNCnVzLnRvdGFsIDwtIGRhdGEudXMudG90YWwgJT4lDQogIG11dGF0ZShuZXcuY29uZmlybWVkID0gaWZlbHNlKGRhdGUgPT0gZGF5MSwgMCwgY29uZmlybWVkIC0gbGFnKGNvbmZpcm1lZCwgbj0xKSksDQogICAgICAgICAgICAgICAgbmV3LmRlYXRocyA9IGlmZWxzZShkYXRlID09IGRheTEsIDAsIGRlYXRocyAtIGxhZyhkZWF0aHMsIG49MSkpKQ0KDQp1cy50b3RhbCAlPD4lIG11dGF0ZShuZXcuY29uZmlybWVkID0gaWZlbHNlKG5ldy5jb25maXJtZWQgPCAwLCAwLCBuZXcuY29uZmlybWVkKSwNCiAgICAgICAgICAgICAgICAgbmV3LmRlYXRocyA9IGlmZWxzZShuZXcuZGVhdGhzIDwgMCwgMCwgbmV3LmRlYXRocykpDQp1cy50b3RhbA0KDQp1cy50b3RhbC5sb25nIDwtIHVzLnRvdGFsICU+JSANCiAgZ2F0aGVyKGtleSA9IHR5cGUsIHZhbHVlID0gY291bnQsIC1jKGRhdGUsIHN0YXRlKSkNCg0KI1ZpZXcodXMudG90YWwubG9uZykNCmBgYA0KDQoNCmBgYHtyfQ0KcGxvdF91cy5jYXNlcyA8LSBkYXRhLnVzLmxvbmcgJT4lIGZpbHRlcihzdGF0ZSA9PSAiVVMiKSAlPiUNCiAgZ2dwbG90KGFlcyh4ID0gZGF0ZSwgeSA9IGNvdW50KSkgKw0KICBnZW9tX2FyZWEoYWVzKGZpbGw9dHlwZSksIGFscGhhPTAuNSkgKw0KIyAgbGFicyh0aXRsZSA9IHBhc3RlMCgiQ3VtdWxhdGl2ZSBjYXNlcyBpbiBVUyA6ICIsIG1pbi5kYXRlLnR4dC51cywgJy0nLCBtYXguZGF0ZS50eHQudXMsICIgKFN0YWNrLCBMb2cgU2NhbGUpIikpICsNCiAgc2NhbGVfeV9jb250aW51b3VzKHRyYW5zPSdsb2cxMCcpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoJ3JlZCcsICdibGFjaycpKSsNCiAgeWxhYigiIikgKw0KICB0aGVtZShheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9OCkpDQojVmlldyh1cy50b3RhbC5sb25nICkNCnBsb3RfdXMubmV3Y29uZiA8LSB1cy50b3RhbC5sb25nICU+JQ0KICBmaWx0ZXIodHlwZSAlaW4lIGMoIm5ldy5jb25maXJtZWQiLCJuZXcuZGVhdGhzIikpICU+JQ0KICBnZ3Bsb3QoYWVzKHggPSBkYXRlLCB5ID0gY291bnQsIGNvbG9yID0gdHlwZSkpICsgDQogIGdlb21fbGluZSgpICsgDQogICNsYWJzKHRpdGxlID0gcGFzdGUwKCJEYWlseSBjb25maXJtZWQgY2FzZXMgaW4gVVMgOiAiLCBtaW4uZGF0ZS50eHQudXMsICctJywgbWF4LmRhdGUudHh0LnVzKSkgKw0KICB4bGFiKCJEYXRlIikgKw0KICB5bGFiKCJDb25maXJtZWQgY2FzZXMiKQ0KDQpnbHkucGxvdF91cy5jYXNlcyA8LSBnZ3Bsb3RseShwbG90X3VzLmNhc2VzKQ0KZ2x5LnBsb3RfdXMuY2FzZXMgPC0gZ2x5LnBsb3RfdXMuY2FzZXMgJT4lIA0KICBsYXlvdXQodGl0bGUgPSBwYXN0ZTAoIjxiPkN1bXVsYXRpdmUgY2FzZXMgaW4gVVMgOiAiLCBtaW4uZGF0ZS50eHQudXMsICctJywgbWF4LmRhdGUudHh0LnVzLCAiIChTdGFjaywgTG9nIFNjYWxlKSIpLCANCiAgICAgICAgIGZvbnQ9bGlzdChzaXplID0gMTApKQ0KZ2x5LnBsb3RfdXMubmV3Y29uZiA8LSBnZ3Bsb3RseShwbG90X3VzLm5ld2NvbmYpICU+JSANCiAgbGF5b3V0KHRpdGxlID0gcGFzdGUwKCI8Yj5EYWlseSBjb25maXJtZWQgY2FzZXMgaW4gVVMgOiAiLCBtaW4uZGF0ZS50eHQudXMsICctJywgbWF4LmRhdGUudHh0LnVzKSwgDQogICAgICAgICBmb250PWxpc3Qoc2l6ZSA9IDEwKSkNCg0KZ2x5LnBsb3RfdXMuY2FzZXMgIA0KZ2x5LnBsb3RfdXMubmV3Y29uZg0KDQpgYGANCg0KYGBge3J9DQpkYXRhLnVzICU8PiUgbXV0YXRlKG5ldy5jb25maXJtZWQgPSBpZmVsc2UoZGF0ZSA9PSBkYXkxLCBOQSwgY29uZmlybWVkIC0gbGFnKGNvbmZpcm1lZCwgbj0xKSksDQogICAgICAgICAgICAgICAgIG5ldy5kZWF0aHMgPSBpZmVsc2UoZGF0ZSA9PSBkYXkxLCBOQSwgZGVhdGhzIC0gbGFnKGRlYXRocywgbj0xKSkpDQoNCmRhdGEudXMgJTw+JSBtdXRhdGUobmV3LmNvbmZpcm1lZCA9IGlmZWxzZShuZXcuY29uZmlybWVkIDwgMCwgMCwgbmV3LmNvbmZpcm1lZCksDQogICAgICAgICAgICAgICAgIG5ldy5kZWF0aHMgPSBpZmVsc2UobmV3LmRlYXRocyA8IDAsIDAsIG5ldy5kZWF0aHMpKQ0KVmlldyhkYXRhLnVzKQ0KDQojZGF0YS51cy5kYWlseSA8LSBkYXRhLnVzICU+JSBmaWx0ZXIoc3RhdGUgPT0gIlVTIikNCiNWaWV3KGRhdGEudXMuZGFpbHkpDQoNCmRhdGEudXMucG9wIDwtIGRmX3VzX3BvcCAlPiUgc2VsZWN0KFN0YXRlLCBQb3B1bGF0aW9uKSAlPiUNCiAgcmVuYW1lKCJzdGF0ZSIgPSAiU3RhdGUiKSANCiNkYXRhLnVzLnBvcA0KDQpkYXRhLnVzLmdlbmRlciA8LSBkZl91c19nZW5kZXIgJT4lIHNlbGVjdChTdGF0ZSwgTWFsZSwgRmVtYWxlKSAlPiUNCiAgcmVuYW1lKCJzdGF0ZSIgPSAiU3RhdGUiKSANCmRhdGEudXMuZ2VuZGVyDQoNCmRhdGEudXMubGF0ZXN0IDwtIGRhdGEudXMgJT4lDQogIGZpbHRlcihkYXRlID09IG1heC5kYXRlLnVzKSAlPiUgDQogIG1lcmdlKGRhdGEudXMucG9wLCBieSA9ICJzdGF0ZSIsIGFsbC54ID0gVCkgJT4lDQogIG1lcmdlKGRhdGEudXMuZ2VuZGVyLCBieSA9ICJzdGF0ZSIsIGFsbC54ID0gVCkNCiNWaWV3KGRhdGEudXMubGF0ZXN0KQ0KDQpkYXRhLnVzLmxhdGVzdCRQb3B1bGF0aW9uW2RhdGEudXMubGF0ZXN0JHN0YXRlID09ICJVUyJdIDwtIHN1bShkYXRhLnVzLnBvcCRQb3B1bGF0aW9uKQ0KZGF0YS51cy5sYXRlc3QgJTw+JSBtdXRhdGUocmFua2luZyA9IGRlbnNlX3JhbmsoZGVzYyhjb25maXJtZWQpKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbmZpcm1lZC5yYXRlID0gKDEwMCAqIGNvbmZpcm1lZCAvIFBvcHVsYXRpb24pICU+JSByb3VuZCgyKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGRlYXRoLnJhdGUgPSAoMTAwICogZGVhdGhzIC8gY29uZmlybWVkKSAlPiUgcm91bmQoMiksDQogICAgICAgICAgICAgICAgICAgICAgICAgICBNYWxlLmNvbmZpcm1lZCA9IChNYWxlICogY29uZmlybWVkKSAlPiUgcm91bmQoMCksDQogICAgICAgICAgICAgICAgICAgICAgICAgICBGZW1hbGUuY29uZmlybWVkID0gRmVtYWxlICogY29uZmlybWVkICU+JSByb3VuZCgwKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgIE1hbGUuZGVhdGhzID0gTWFsZSAqIGRlYXRocywNCiAgICAgICAgICAgICAgICAgICAgICAgICAgIEZlbWFsZS5kZWF0aHMgPSBGZW1hbGUgKiBkZWF0aHMpICU+JQ0KICBhcnJhbmdlKHJhbmtpbmcpIA0KDQp0b3AudXMgPC0gZGF0YS51cy5sYXRlc3RbLDFdDQp0b3AudXMNCg0KI1ZpZXcodG9wLnVzKQ0KDQpkYXRhLnVzLmxhdGVzdC5zaG93IDwtIGRhdGEudXMubGF0ZXN0ICU+JSBzZWxlY3Qoc3RhdGUsZGF0ZSxjb25maXJtZWQsZGVhdGhzLCBjb25maXJtZWQucmF0ZSwgZGVhdGgucmF0ZSkNCiNWaWV3KGRhdGEudXMubGF0ZXN0LnNob3cpDQpgYGANCg0KYGBge3J9DQojIExpc3Qgb2YgdG9wIDIwIHN0YXRlDQprIDwtIDIwDQpkYXRhLnVzLnRvcCA8LSBkYXRhLnVzLmxhdGVzdCAlPiUNCiAgZmlsdGVyKHJhbmtpbmcgPD0gaysxKSAlPiUgDQogIGFycmFuZ2UocmFua2luZykgDQojVmlldyhkYXRhLnVzLnRvcCkNCg0KdXMuc3RhdGUudG9wIDwtIGRhdGEudXMudG9wICU+JSBwdWxsKHN0YXRlKSAlPiUgYXMuY2hhcmFjdGVyKCkNCnVzLnN0YXRlLnRvcCAgJT4lIHNldGRpZmYoJ1VTJykgJT4lIHByaW50KCkNCg0KIyBjb25maXJtZWQgcmF0ZSAmIGRlYXRoIHJhdGUgb2YgdG9wIDIwIHN0YXRlDQpnLnJhdGUgPC0gZGF0YS51cy5sYXRlc3QgJT4lIGZpbHRlcihzdGF0ZSAlaW4lIHVzLnN0YXRlLnRvcCAmIHN0YXRlICE9ICJVUyIpICU+JQ0KICBzZWxlY3Qoc3RhdGUsIGNvbmZpcm1lZC5yYXRlLCBkZWF0aC5yYXRlLCByYW5raW5nKSAlPiUNCiAgZ2F0aGVyKGtleSA9IFR5cGUsIHZhbHVlID0gUGVyY2VudCwgLWMoc3RhdGUsIHJhbmtpbmcpKSAlPiUNCiAgZ2dwbG90KGFlcyh4PXJlb3JkZXIoc3RhdGUsIC1kZXNjKHJhbmtpbmcpKSwgeT1QZXJjZW50LCBmaWxsID0gUGVyY2VudCkpICsNCiAgZ2VvbV9iYXIoc3RhdD0naWRlbnRpdHknKSArDQogIHNjYWxlX2ZpbGxfZ3JhZGllbnQobG93ID0gIiNlYmJjNjIiLCBoaWdoID0gIiNiNDIwMDYiKSArDQogIGdlb21fdGV4dChhZXMobGFiZWw9UGVyY2VudCwgeT1QZXJjZW50KSwgc2l6ZT0zLCB2anVzdD0wKSArDQogIHhsYWIoJycpICsgeWxhYignJykgKw0KICBsYWJzKHRpdGxlPXBhc3RlMCgnQ29uZmlybWVkIFJhdGUgJiBEZWF0aCBSYXRlIG9mIFRvcCAyMCBTdGF0ZSBpbiBVUycpKSArDQogICNzY2FsZV9maWxsX2NvbnRpbnVvdXMobmFtZT0nU3RhdGUnLCBsYWJlbHM9YWVzKFBlcmNlbnQpKSArDQogIHRoZW1lKGxlZ2VuZC50aXRsZT1lbGVtZW50X2JsYW5rKCksDQogICAgICAgIGxlZ2VuZC5wb3NpdGlvbj0nbm9uZScsDQogICAgICAgIHBsb3QudGl0bGU9ZWxlbWVudF90ZXh0KHNpemU9MTMpLA0KICAgICAgICBheGlzLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9OCksDQogICAgICAgIGF4aXMudGV4dC54PWVsZW1lbnRfdGV4dChhbmdsZT00NSwgaGp1c3Q9MSkpICsNCiAgZmFjZXRfd3JhcCh+VHlwZSwgbmNvbD0xLCBzY2FsZXM9J2ZyZWVfeScpIA0KDQpnLnJhdGUNCmBgYA0KDQpgYGB7cn0NCnVzLmNvbmZpcm1lZC5udW0gPC0gZGF0YS51cy50b3AkY29uZmlybWVkW2RhdGEudXMudG9wJHN0YXRlID09ICJVUyJdICU+JSBhcy5udW1lcmljKCkNCnVzLmNvbmZpcm1lZC5udW0NCg0KZGF0YS51cy50b3AgJTw+JSBtdXRhdGUoY29uZmlybWVkLnBlci51cyA9IChjb25maXJtZWQgKiAxMDAgLyB1cy5jb25maXJtZWQubnVtKSAlPiUgcm91bmQoMSkpDQpkYXRhLnVzLnRvcA0KDQpnLnVzLnRvcDEgPC0gZGF0YS51cy50b3AgJT4lDQogIGZpbHRlcihzdGF0ZSAhPSAiVVMiKSAlPiUNCiAgc2VsZWN0KHN0YXRlLCBjb25maXJtZWQsIGNvbmZpcm1lZC5wZXIudXMsIHJhbmtpbmcpICU+JQ0KICBnYXRoZXIoa2V5ID0gVHlwZSwgdmFsdWUgPSBjb3VudCwgLWMoc3RhdGUsIHJhbmtpbmcsIGNvbmZpcm1lZC5wZXIudXMpKSAlPiUNCiAgYXJyYW5nZShyYW5raW5nKSAlPiUNCiAgDQogIGdncGxvdChhZXMoeCA9IHJlb3JkZXIoc3RhdGUsIHJhbmtpbmcpLCB5ID0gY291bnQsIGZpbGwgPSBjb3VudCkpICsgDQogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiKSArDQogICNsYWJzKHRpdGxlID0gIjIwIHN0YXRlIGluIFVTIHdpdGggbW9zdCBjb25maXJtZWQgY2FzZXMiKSArDQogIHNjYWxlX2NvbG9yX2dyYWRpZW50KGxvdyA9ICIjOTNEQkZGIiwgaGlnaCA9ICIjRkY3NzcxIikgKw0KICB4bGFiKCIiKSArIA0KICB5bGFiKCJDb25maXJtZWQgQ2FzZXMiKSArDQogIGdlb21fdGV4dChhZXMobGFiZWw9cGFzdGUwKGNvbmZpcm1lZC5wZXIudXMsICIlIikpLCBzaXplPTMsIHZqdXN0PS0wLjUpICsNCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCx2anVzdCA9IDAuNSkpKw0KICB0aGVtZSgNCiAgICBheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksDQogICAgYXhpcy50aWNrcyA9IGVsZW1lbnRfYmxhbmsoKSwNCiAgICBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpLA0KICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksDQogICAgcGFuZWwuYm9yZGVyID0gZWxlbWVudF9ibGFuaygpLA0KICAgIHBhbmVsLmJhY2tncm91bmQgPSBlbGVtZW50X2JsYW5rKCksDQogICAgI2xlZ2VuZC5wb3NpdGlvbiA9ICJub25lIg0KICApDQoNCmcudXMudG9wMQ0KDQpnbHkudXMudG9wMSA8LSBnZ3Bsb3RseShnLnVzLnRvcDEpDQpnbHkudXMudG9wMQ0KYGBgDQoNCmBgYHtyfQ0KZGF0YS51cy5sYXRlc3QuY29yciA8LSBkZl91cyAlPiVmaWx0ZXIoZGF0ZT09bWF4LmRhdGUudXMpDQojVmlldyhkYXRhLnVzLmxhdGVzdC5jb3JyKQ0KDQpkYXRhLnVzLmxhdGVzdC5jb3JyICU8PiUgbXV0YXRlKHJhbmtpbmcgPSBkZW5zZV9yYW5rKGRlc2MoY29uZmlybWVkKSkpICU+JSBhcnJhbmdlKHJhbmtpbmcpDQp0b3AudXMgPC0gZGF0YS51cy5sYXRlc3QuY29yclssMl0NCiN0b3AudXMNCg0KIyBMaXN0IG9mIHRvcCAyMCBzdGF0ZQ0KayA8LSAxOQ0KZGF0YS51cy5sYXRlc3QudG9wIDwtIGRhdGEudXMubGF0ZXN0LmNvcnIgJT4lDQogICAgZmlsdGVyKHJhbmtpbmcgPD0gaysxKSAlPiUgDQogICAgYXJyYW5nZShyYW5raW5nKSANCiNWaWV3KGRhdGEudXMubGF0ZXN0LnRvcCkNCnVzLnN0YXRlLnRvcCA8LSBkYXRhLnVzLmxhdGVzdC50b3AgJT4lIHB1bGwoc3RhdGUpICU+JSBhcy5jaGFyYWN0ZXIoKQ0KdXMuc3RhdGUudG9wICAlPiUgc2V0ZGlmZignVVMnKSAlPiUgcHJpbnQoKQ0KI1ZpZXcodXMuc3RhdGUudG9wKQ0KYGBgDQoNCg0KDQpgYGB7cn0NCiNnZW5kZXIgaW4gdXMNCmRmX3VzX2dlbmRlciA8LSB0YmwobXlfZGIsIHNxbCgic2VsZWN0ICogZnJvbSBkYXRhX2dlbmRlciIpKQ0KZGZfdXNfZ2VuZGVyIDwtIGFzLmRhdGEuZnJhbWUoZGZfdXNfZ2VuZGVyKQ0KZGZfdXNfZ2VuZGVyIDwtIHNlbGVjdChkZl91c19nZW5kZXIsYygiU3RhdGUiLCJNYWxlIiwiRmVtYWxlIikpDQpkZl91c19nZW5kZXIgPC0gcmVuYW1lKGRmX3VzX2dlbmRlciwic3RhdGUiPSJTdGF0ZSIpDQpkZl91c19nZW5kZXINCg0KI3BvcHVsYXRpb24gaW4gdXMNCmRmX3VzX3BvcCA8LSB0YmwobXlfZGIsIHNxbCgic2VsZWN0ICogZnJvbSBkYXRhX3BvcHVsYXRpb24iKSkNCmRmX3VzX3BvcCA8LSBhcy5kYXRhLmZyYW1lKGRmX3VzX3BvcCkNCmRmX3VzX3BvcCA8LSBzZWxlY3QoZGZfdXNfcG9wLGMoIlN0YXRlIiwiUG9wdWxhdGlvbiIpKQ0KZGZfdXNfcG9wIDwtIHJlbmFtZShkZl91c19wb3AsInN0YXRlIj0iU3RhdGUiKQ0KZGZfdXNfcG9wDQoNCiNsb2NrZG93biBpbiB1cw0KZGZfdXNfbG9ja2Rvd24gPC0gdGJsKG15X2RiLCBzcWwoInNlbGVjdCAqIGZyb20gZGF0YV9sb2NrZG93biIpKQ0KZGZfdXNfbG9ja2Rvd24gPC0gYXMuZGF0YS5mcmFtZShkZl91c19sb2NrZG93bikNCmRmX3VzX2xvY2tkb3duIDwtIHNlbGVjdChkZl91c19sb2NrZG93bixjKCJTdGF0ZSIsIkRheSBsb2NrZG93biIpKQ0KZGZfdXNfbG9ja2Rvd24gPC0gcmVuYW1lKGRmX3VzX2xvY2tkb3duLCJzdGF0ZSI9IlN0YXRlIikNCmRmX3VzX2xvY2tkb3duDQoNCiNHRFAgaW4gdXMNCmRmX3VzX2dkcCA8LSB0YmwobXlfZGIsIHNxbCgic2VsZWN0ICogZnJvbSB1c19nZHAiKSkNCmRmX3VzX2dkcCAgPC0gYXMuZGF0YS5mcmFtZShkZl91c19nZHAgKQ0KZGZfdXNfZ2RwICA8LSBzZWxlY3QoZGZfdXNfZ2RwICxjKCJTdGF0ZSIsIkdEUHMiKSkNCmRmX3VzX2dkcCAgPC0gcmVuYW1lKGRmX3VzX2dkcCAsInN0YXRlIj0iU3RhdGUiKQ0KZGZfdXNfZ2RwIA0KDQojaG9tZWxlc3MgaW4gdXMNCmRmX3VzX2hvbWVsZXNzIDwtIHRibChteV9kYiwgc3FsKCJzZWxlY3QgKiBmcm9tIHVzX2hvbWVsZXNzIikpDQpkZl91c19ob21lbGVzcyA8LSBhcy5kYXRhLmZyYW1lKGRmX3VzX2hvbWVsZXNzKQ0KZGZfdXNfaG9tZWxlc3MgPC0gc2VsZWN0KGRmX3VzX2hvbWVsZXNzLGMoIlN0YXRlIiwiSG9tZWxlc3MiKSkNCmRmX3VzX2hvbWVsZXNzIDwtIHJlbmFtZShkZl91c19ob21lbGVzcywic3RhdGUiPSJTdGF0ZSIpDQpkZl91c19ob21lbGVzcw0KDQojVmlldyhkZl91cykNCg0KDQpgYGANCmBgYHtyfQ0KI21lcmdlDQptZXJnY291bnRyeSA9IGZ1bmN0aW9uKGRhdGExLGRhdGEyKXsNCiAgZGF0YSA8LSBtZXJnZSh4ID0gZGF0YTEsIHkgPSBkYXRhMiwgYnkgPSAic3RhdGUiLCBhbGwueCA9IFRSVUUpIA0KICByZXR1cm4oZGF0YSkNCn0NCg0KZGF0YS50b3Auc3RhdGUgPC0gZGF0YS51cy5sYXRlc3QudG9wICU+JSBmaWx0ZXIoc3RhdGUgIT0gIlVTIikgJT4lDQogIHNlbGVjdChzdGF0ZSwgY29uZmlybWVkLCBkZWF0aHMpDQojVmlldyhkYXRhLnRvcC5zdGF0ZSkNCg0KZGZfQWxsdXMgPC0gbWVyZ2NvdW50cnkoZGF0YS50b3Auc3RhdGUsIGRmX3VzX2dlbmRlcikNCg0KZGZfQWxsdXMgPC0gbWVyZ2NvdW50cnkoZGZfQWxsdXMsIGRmX3VzX3BvcCkgDQoNCmRmX0FsbHVzIDwtIG1lcmdjb3VudHJ5KGRmX0FsbHVzLCBkZl91c19sb2NrZG93bikgDQoNCmRmX0FsbHVzIDwtIG1lcmdjb3VudHJ5KGRmX0FsbHVzLCBkZl91c19nZHApDQoNCmRmX0FsbHVzIDwtIG1lcmdjb3VudHJ5KGRmX0FsbHVzLCBkZl91c19ob21lbGVzcykNCg0KI1ZpZXcoZGZfQWxsdXMpDQpgYGANCg0KYGBge3J9DQppbmRleCA8LSBpcy5uYShkZl9BbGx1cykNCmRmX0FsbHVzW2luZGV4XSA8LSAwDQojVmlldyhkZl9BbGx1cykNCg0Kbm9ybWFsaXplID0gZnVuY3Rpb24oZGF0YSl7DQogICNyZXR1cm4gKChkYXRhIC0gbWluKGRhdGEsbmEucm0gPSBUUlVFKSkvKG1heChkYXRhLG5hLnJtID0gVFJVRSkgLSBtaW4oZGF0YSxuYS5ybSA9IFRSVUUpKSkNCiAgeiA8LSBzY2FsZShkYXRhKTsNCiAgdGFuaCh6LzIpDQp9DQpBbGx1cyA9IGFzLmRhdGEuZnJhbWUoYXBwbHkoZGZfQWxsdXNbLDI6OV0sMixub3JtYWxpemUpKQ0KY29ycl9kYXRhVVMgPC0gQWxsdXMgDQoNCiNWaWV3KGNvcnJfZGF0YVVTKQ0KYGBgDQoNCmBgYHtyfQ0KQWxsdXMkc3RhdGUgPC0gYyhkZl9BbGx1cyRzdGF0ZSkNCiNWaWV3KEFsbHVzKQ0KDQpjb3JyX2RhdGFVUyA8LSByZW5hbWUoY29ycl9kYXRhVVMgLCJEYXlsb2NrZG93biIgPSAiRGF5IGxvY2tkb3duIikNCmBgYA0KDQpgYGB7cn0NCkFsbHVzX3Bsb3QgPC0gc2VsZWN0KEFsbHVzLCJzdGF0ZSIsImNvbmZpcm1lZCIsImRlYXRocyIsIk1hbGUiLCJGZW1hbGUiLCJQb3B1bGF0aW9uIiwiRGF5IGxvY2tkb3duIiwiR0RQcyIsIkhvbWVsZXNzIikNCkFsbHVzX3Bsb3QgJTw+JSBnYXRoZXIoa2V5PXR5cGUsIHZhbHVlPWNvdW50LCAtYyhzdGF0ZSkpDQpsZXZlbF9vcmRlciA8LSBmYWN0b3IoQWxsdXNfcGxvdCR0eXBlLCANCiAgICAgICAgICAgICAgICAgICAgICBsZXZlbCA9IGMoImNvbmZpcm1lZCIsImRlYXRocyIsIk1hbGUiLCJGZW1hbGUiLCJQb3B1bGF0aW9uIiwiRGF5IGxvY2tkb3duIiwiR0RQcyIsIkhvbWVsZXNzIikpDQoNCkFsbHVzX0hlYXQgPC0gQWxsdXNbLGMobmNvbChBbGx1cyksMToobmNvbChBbGx1cyktMSkpXQ0Kcm93bmFtZXMoQWxsdXNfSGVhdCkgPC0gQWxsdXNfSGVhdFssMV0NCkFsbHVzX0hlYXRbLDFdIDwtIE5VTEwNCmhlYXRtYXAuMihhcy5tYXRyaXgoQWxsdXNfSGVhdCksDQogICAgICAgICAgICAgICAgICAgIHNjYWxlPSJub25lIiwgDQogICAgICAgICAgICAgICAgICAgIGNvbCA9IGNvbG9yUmFtcFBhbGV0dGUoYygiIzZEOUVDMSIsIndoaXRlIiwiI0U0NjcyNiIpKShuID0gMjAwKSwNCiAgICAgICAgICAgICAgICAgICAgbWFyZ2lucz1jKDEwLDYpLHRyYWNlPSJjb2x1bW4iKQ0KDQpoZWF0bWFwdXMgPC0gZ2dwbG90KGRhdGEgPSBBbGx1c19wbG90LCBhZXMoeD1zdGF0ZSwgeT1sZXZlbF9vcmRlciwgZmlsbD1jb3VudCkpICsgDQogIGdlb21fdGlsZSgpICsNCiAgc2NhbGVfZmlsbF9ncmFkaWVudChsb3cgPSAicGluayIsIGhpZ2ggPSAiYmx1ZSIpICsNCiAgeGxhYigiIikgKw0KICB5bGFiKCIiKSArDQogIHRoZW1lX2J3KCkgKw0KICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLHZqdXN0ID0gMSkpKw0KICB0aGVtZSgNCiAgICBheGlzLmxpbmUgPSBlbGVtZW50X2JsYW5rKCksDQogICAgYXhpcy50aWNrcyA9IGVsZW1lbnRfYmxhbmsoKSwNCiAgICBwYW5lbC5ncmlkLm1pbm9yID0gZWxlbWVudF9ibGFuaygpLA0KICAgIHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksDQogICAgcGFuZWwuYm9yZGVyID0gZWxlbWVudF9ibGFuaygpLA0KICAgIHBhbmVsLmJhY2tncm91bmQgPSBlbGVtZW50X2JsYW5rKCksDQogICAgI2xlZ2VuZC5wb3NpdGlvbiA9ICJub25lIg0KICApICtsYWJzKHRpdGxlPSdUaGUgaGVhdG1hcCBvZiBDT1ZJRC0xOSBpbmZlY3Rpb25zIGluIFVTJykNCg0KaGVhdG1hcHVzDQpgYGANCg0KDQpgYGB7cn0NCiNjb3JyZWxhdGlvbg0KY29ycl9kYXRhVVMgJTw+JSBzZWxlY3QoYyhjb25maXJtZWQsZGVhdGhzLE1hbGUsRmVtYWxlLFBvcHVsYXRpb24sRGF5bG9ja2Rvd24sR0RQcyxIb21lbGVzcykpDQpoZWFkKGNvcnJfZGF0YVVTKQ0KY29yKGNvcnJfZGF0YVVTKQ0KZ2djb3JycGxvdChjb3IoY29ycl9kYXRhVVMpLGhjLm9yZGVyID0gVFJVRSwNCiAgICAgICAgICAgb3V0bGluZS5jb2xvciA9ICJ3aGl0ZSIsDQogICAgICAgICAgIGNvbG9ycyA9IGMoIiM2RDlFQzEiLCJ3aGl0ZSIsIiNFNDY3MjYiKSwNCiAgICAgICAgICAgbGFiID0gVFJVRSkNCg0KDQoNCmBgYA0KDQpgYGB7cn0NCiNDb3ZpZF9UaGFpbGFuZA0KZGZfdGhhaSA8LSB0YmwobXlfZGIsIHNxbCgic2VsZWN0ICogZnJvbSBjb3ZpZDE5X3RoYWlsYW5kIikpDQpkZl90aGFpIDwtIGFzLmRhdGEuZnJhbWUoZGZfdGhhaSkNCiNWaWV3KGRmX3RoYWkpDQpgYGANCg0KYGBge3J9DQojY2xlYW4gQ292aWRfVGhhaWxhbmQNCmRhdGVzLnRoIDwtIGRmX3RoYWlbLDJdJT4lIG1keSgpDQpyYW5nZShkYXRlcy50aCkNCm1pbi5kYXRlLnRoIDwtIG1pbihkYXRlcy50aCkNCm1heC5kYXRlLnRoIDwtIG1heChkYXRlcy50aCkNCm1pbi5kYXRlLnR4dC50aCA8LSBtaW4uZGF0ZS50aCAlPiUgZm9ybWF0KCclZCAlYiAlWScpDQptYXguZGF0ZS50eHQudGggPC0gbWF4LmRhdGUudGggJT4lIGZvcm1hdCgnJWQgJWIgJVknKQ0KYGBgDQoNCmBgYHtyfQ0KZGZfdGhhaSRhbm5vdW5jZV9kYXRlIDwtIG1keShkZl90aGFpJGFubm91bmNlX2RhdGUpDQpkZl90aGFpJG5vdGlmaWNhdGlvbl9kYXRlIDwtIG1keShkZl90aGFpJG5vdGlmaWNhdGlvbl9kYXRlKQ0KZGZfdGhhaQ0KYGBgDQoNCmBgYHtyfQ0KZGZfdGhhaSA8LSBkZl90aGFpICU+JSBzZWxlY3QoIU5vLikgJT4lIHNlbGVjdCghbm90aWZpY2F0aW9uX2RhdGUpICU+JSANCiAgZ3JvdXBfYnkoYW5ub3VuY2VfZGF0ZSkNCmRmX3RoYWkNClZpZXcoZGZfdGhhaSkNCmBgYA0KDQpgYGB7cn0NCiMgVG90YWwgY29uZmlybWVkIGNhc2VzIGluIFRoYWlsYW5kDQpkYXRhLnRoYWkuY291bnQgPC0gZGZfdGhhaSAlPiUNCiAgc2VsZWN0KGFubm91bmNlX2RhdGUpICU+JQ0KICBzdW1tYXJpc2UoY29uZmlybWVkID0gbigpKSAlPiUgYXMuZGF0YS5mcmFtZSgpDQpWaWV3KGRhdGEudGhhaS5jb3VudCkNCmRhdGEudGhhaS5jb3VudCRjdW11bGF0aXZlX2NvbmZpcm1lZCA8LSBjdW1zdW0oZGF0YS50aGFpLmNvdW50WywgMl0pDQpgYGANCmBgYHtyfQ0KZGF0YS50aGFpLmNvdW50IDwtIGRhdGEudGhhaS5jb3VudCAlPiUgcmVuYW1lKG5ldy5jb25maXJtZWQ9Y29uZmlybWVkKSAlPiUgcmVuYW1lKGNvbmZpcm1lZD1jdW11bGF0aXZlX2NvbmZpcm1lZCkNCnRoYWkudG90YWwubG9uZyA8LSBkYXRhLnRoYWkuY291bnQgJT4lIA0KICBnYXRoZXIoa2V5ID0gdHlwZSwgdmFsdWUgPSBjb3VudCwgLWMoYW5ub3VuY2VfZGF0ZSkpDQp0aGFpLnRvdGFsLmxvbmcNCmBgYA0KDQpgYGB7cn0NCg0KcGxvdF90aGFpbGFuZC5jb25mIDwtIHRoYWkudG90YWwubG9uZyAlPiUNCiAgZmlsdGVyKHR5cGUgJWluJSBjKCJuZXcuY29uZmlybWVkIikpICU+JQ0KICBnZ3Bsb3QoYWVzKHggPSBhbm5vdW5jZV9kYXRlLCB5ID0gY291bnQsIGNvbG9yID0gdHlwZSkpICsgDQogIGdlb21fbGluZSgpICsgDQogIHhsYWIoIkRhdGUiKSArDQogIHlsYWIoImNhc2VzIikNCg0KZ2x5LnBsb3RfdGhhaS5jYXNlcyA8LSBnZ3Bsb3RseShwbG90X3RoYWlsYW5kLmNvbmYpDQpnbHkucGxvdF90aGFpLmNhc2VzIDwtIGdseS5wbG90X3RoYWkuY2FzZXMgJT4lIA0KICBsYXlvdXQodGl0bGUgPSBwYXN0ZTAoIjxiPkRhaWx5IFRoYWkgQ29uZmlybWVkIENhc2VzIGluIFRoYWlsYW5kIDogSmFuIDIwMjAgLSBKYW4gMjAyMSIpLCANCiAgICAgICAgIGZvbnQ9bGlzdChzaXplID0gMTApKQ0KDQpnbHkucGxvdF90aGFpLmNhc2VzDQpgYGANCmBgYHtyfQ0KcGxvdF90aGFpbGFuZC5jdW11bGF0aXZlIDwtIHRoYWkudG90YWwubG9uZyAlPiUgDQogIGZpbHRlcih0eXBlICVpbiUgYygiY29uZmlybWVkIikpICU+JQ0KICBnZ3Bsb3QoYWVzKHggPSBhbm5vdW5jZV9kYXRlLCB5ID0gY291bnQpKSArDQogIGdlb21fYXJlYShhZXMoZmlsbD10eXBlKSwgYWxwaGE9MC41KSArDQoNCiAgc2NhbGVfeV9jb250aW51b3VzKHRyYW5zPSdsb2cxMCcpICsNCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWMoJ3JlZCcpKSsNCiAgeWxhYigiIikgKw0KICB0aGVtZShheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9OCkpDQoNCmdseS5wbG90X3RoYWkuY3VtdWxhdGl2ZSA8LSBnZ3Bsb3RseShwbG90X3RoYWlsYW5kLmN1bXVsYXRpdmUpDQpnbHkucGxvdF90aGFpLmN1bXVsYXRpdmUgPC0gZ2x5LnBsb3RfdGhhaS5jdW11bGF0aXZlICU+JSANCiAgbGF5b3V0KHRpdGxlID0gcGFzdGUwKCI8Yj5DdW11bGF0aXZlIENhc2VzIGluIFRoYWlsYW5kIDogSmFuIDIwMjAgLSBKYW4gMjAyMSAoU3RhY2ssTG9nIHNjYWxlKSIpLCANCiAgICAgICAgIGZvbnQ9bGlzdChzaXplID0gMTApKQ0KDQpnbHkucGxvdF90aGFpLmN1bXVsYXRpdmUNCmBgYA0KDQoNCg0KYGBge3J9DQojIyBUaGFpIENvbmZpcm1lZCBDYXNlcyAoSmFuIDIwMjAgLSBKYW4gMjAyMQ0KcGxvdDEgPC0gZ2dwbG90KGRhdGEudGhhaS5jb3VudCwgYWVzKHg9YW5ub3VuY2VfZGF0ZSwgeT1jdW11bGF0aXZlX2NvbmZpcm1lZCkpICsNCiAgZ2VvbV9wb2ludCgpICsgZ2VvbV9zbW9vdGgoKSArDQogIHhsYWIoIiAiKSArIHlsYWIoIkNvdW50IikgKyBsYWJzKHRpdGxlPSdUaGFpIEN1bXVsYXRpdmUgQ29uZmlybWVkIENhc2VzIChKYW4gMjAyMCAtIEphbiAyMDIxKScpICsNCiAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdD0xKSkNCnBsb3QyIDwtIGdncGxvdChkYXRhLnRoYWkuY291bnQsIGFlcyh4PWFubm91bmNlX2RhdGUsIHk9Y29tZmlybWVkKSkgKw0KICBnZW9tX3BvaW50KCkgKyBnZW9tX3Ntb290aCgpICsNCiAgICB4bGFiKCIgIikgKyB5bGFiKCJDb3VudCIpKyBsYWJzKHRpdGxlPSdUaGFpIENvbmZpcm1lZCBDYXNlcyAoSmFuIDIwMjAgLSBKYW4gMjAyMSBsb2cgc2NhbGUpJykgKw0KICAgdGhlbWUoYXhpcy50ZXh0Lng9ZWxlbWVudF90ZXh0KGFuZ2xlPTQ1LCBoanVzdD0xKSkrc2NhbGVfeV9jb250aW51b3VzKHRyYW5zPSdsb2cxMCcpDQojIyBzaG93IHR3byBwbG90cyBzaWRlIGJ5IHNpZGUNCmdyaWQuYXJyYW5nZShwbG90MSwgcGxvdDIsIG5jb2w9MSkNCmBgYA0KDQoNCg0KYGBge3J9DQojIENvbmZpcm1lZCBjYXNlcyBkaXZpZGVkIGJ5IHNleCAoZ2VuZGVyKQ0KZGZfdGhhaSRzZXhbZGZfdGhhaSRzZXggPT0gIiJdIDwtICJVbmtub3duIg0KZGF0YS50aGFpLmdlbmRlciA8LSBkZl90aGFpICU+JQ0KICBncm91cF9ieShzZXgpICU+JQ0KICBzdW1tYXJpc2UoY291bnQgPSBuKCkpICU+JQ0KICBtdXRhdGUocGVyY2VudCA9IChjb3VudCAvIHN1bShjb3VudCkgKiAxMDApICU+JSByb3VuZCgyKSkgJT4lDQogIGZpbHRlcihwZXJjZW50PjEpJT4lDQogICNtdXRhdGUocG9zID0gY3Vtc3VtKHBlcmNlbnQpIC0gMC41KnBlcmNlbnQpICU+JQ0KICBhcnJhbmdlKGRlc2MocGVyY2VudCkpDQoNCmRhdGEudGhhaS5nZW5kZXINCmBgYA0KDQpgYGB7cn0NCmRhdGEudGhhaS5nZW5kZXIkc2V4IDwtIGZhY3RvcihkYXRhLnRoYWkuZ2VuZGVyJHNleCwgbGV2ZWxzID0gYXMuY2hhcmFjdGVyKGRhdGEudGhhaS5nZW5kZXIkc2V4KSkNCmRhdGEudGhhaS5nZW5kZXIkc2V4DQpnLnRoLmdlbmRlciA8LSBkYXRhLnRoYWkuZ2VuZGVyICU+JSANCiAgZ2dwbG90KGFlcyh4ID0gIiIsIHkgPSBwZXJjZW50LCBmaWxsID0gc2V4KSkgKw0KICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5Iiwgd2lkdGggPSAxKSArDQogIGNvb3JkX3BvbGFyKCJ5IikgKw0KICB0aGVtZV92b2lkKCkgKw0KICBsYWJzKHRpdGxlPSdHZW5kZXIgb2YgVGhhaSBDb25maXJtZWQgQ2FzZXMgKEphbiAyMDIwIC0gSmFuIDIwMjEpJykrDQogIGdlb21fdGV4dChhZXMobGFiZWwgPSBwYXN0ZTAocGVyY2VudCwgIiUiKSksIGNvbG9yID0gIndoaXRlIiwgc2l6ZSA9IDUsIHBvc2l0aW9uID0gcG9zaXRpb25fc3RhY2sodmp1c3QgPSAwLjUpKSArDQogIGd1aWRlcyhmaWxsID0gZ3VpZGVfbGVnZW5kKHJldmVyc2UgPSBUUlVFKSkgDQpnLnRoLmdlbmRlcg0KYGBgDQoNCmBgYHtyfQ0KIyBDb25maXJtZWQgY2FzZXMgZGl2aWRlZCBieSBwcm92aW5jZV9vZl9vbnNldA0KZGF0YS50aGFpLm9uc2V0IDwtIGRmX3RoYWkgJT4lDQogICAgZ3JvdXBfYnkocHJvdmluY2Vfb2Zfb25zZXQpICU+JQ0KICAgIHN1bW1hcmlzZShjb3VudCA9IG4oKSkgJT4lDQogICAgYXJyYW5nZShkZXNjKGNvdW50KSklPiUNCiAgICByZW5hbWUob25zZXQ9ImNvdW50IiklPiUNCiAgICByZW5hbWUocHJvdmluY2U9InByb3ZpbmNlX29mX29uc2V0IikNCmRhdGEudGhhaS5vbnNldCRwcm92aW5jZVtkYXRhLnRoYWkub25zZXQkcHJvdmluY2UgPT0gIiJdIDwtICJVbmtub3duIg0KZGF0YS50aGFpLm9uc2V0DQpgYGANCg0KYGBge3J9DQojIENvbmZpcm1lZCBjYXNlcyBkaXZpZGVkIGJ5IHByb3ZpbmNlX29mX2lzb2xhdGlvbiANCmRhdGEudGhhaS5pc29sYXRpb24gPC0gZGZfdGhhaSAlPiUNCiAgICBncm91cF9ieShwcm92aW5jZV9vZl9pc29sYXRpb24pICU+JQ0KICAgIHN1bW1hcmlzZShjb3VudCA9IG4oKSkgJT4lDQogICAgYXJyYW5nZShkZXNjKGNvdW50KSklPiUNCiAgICByZW5hbWUoaXNvbGF0aW9uID0iY291bnQiKSU+JQ0KICAgIHJlbmFtZShwcm92aW5jZT0icHJvdmluY2Vfb2ZfaXNvbGF0aW9uIikNCg0KZGF0YS50aGFpLmlzb2xhdGlvbiRwcm92aW5jZVtkYXRhLnRoYWkuaXNvbGF0aW9uJHByb3ZpbmNlID09ICIiXSA8LSAiVW5rbm93biINCmRhdGEudGhhaS5pc29sYXRpb24NCmBgYA0KDQpgYGB7cn0NCmRhdGEucHJvdmluY2UgPC0gbWVyZ2UoeCA9IGRhdGEudGhhaS5vbnNldCwgeSA9IGRhdGEudGhhaS5pc29sYXRpb24sIGJ5ID0gInByb3ZpbmNlIiwgYWxsLnggPSBUUlVFKSANCmRhdGEucHJvdmluY2UgPC0gZGF0YS5wcm92aW5jZSU+JWZpbHRlcihvbnNldD4xMjUpDQpkYXRhLnByb3ZpbmNlDQpgYGANCg0KYGBge3J9DQpkYXRhLnByb3ZpbmNlLmxvbmcgPC0gZGF0YS5wcm92aW5jZSAlPiUNCiAgICBnYXRoZXIoa2V5PXR5cGUsIHZhbHVlPWNvdW50LCAtYyhwcm92aW5jZSkpDQpkYXRhLnByb3ZpbmNlLmxvbmcNCmBgYA0KDQpgYGB7cn0NCnByb3ZpbmNlLnRoYWkgPC0gZ2dwbG90KGRhdGEucHJvdmluY2UubG9uZyxhZXMoeD1wcm92aW5jZSx5PWNvdW50KSkrDQogICAgZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIscG9zaXRpb24gPSAiZG9kZ2UiLGFlcyhmaWxsPXR5cGUpKSsNCiAgICBsYWJzKHRpdGxlID0gIlByb3ZpbmNlIG9mIG9uc2V0IGFuZCBQcm92aW5jZSBvZiBpc29sYXRpb24gaW4gVGhhaWxhbmQiKSsNCiAgICB0aGVtZShsZWdlbmQudGl0bGU9ZWxlbWVudF9ibGFuaygpLA0KICAgICAgICAgICAgICAgICAgI2xlZ2VuZC5wb3NpdGlvbj0nbm9uZScsDQogICAgICAgICAgICAgICAgICBwbG90LnRpdGxlPWVsZW1lbnRfdGV4dChzaXplPTEyKSwNCiAgICAgICAgICAgICAgICAgIGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT0xMCksDQogICAgICAgICAgICAgICAgICBheGlzLnRleHQueD1lbGVtZW50X3RleHQoYW5nbGU9NDUsIGhqdXN0PTEpKSsgeGxhYignJykgKyB5bGFiKCdDb3VudCcpDQpnbHkucHJvdmluY2UudGhhaSA8LSBnZ3Bsb3RseShwcm92aW5jZS50aGFpKQ0KZ2x5LnByb3ZpbmNlLnRoYWkNCmBgYA0KDQpgYGB7cn0NCiMgQ29uZmlybWVkIGNhc2VzIGRpdmlkZWQgYnkgcmlzaw0KZGF0YS50aGFpLnJpc2sgPC0gZGZfdGhhaSAlPiUNCiAgZ3JvdXBfYnkocmlzaykgJT4lDQogIHN1bW1hcmlzZShjb3VudCA9IG4oKSkgJT4lDQogIGFycmFuZ2UoZGVzYyhjb3VudCkpDQpkYXRhLnRoYWkucmlzaw0KYGBgDQoNCmBgYHtyfQ0KIyBDb25maXJtZWQgY2FzZXMgZGl2aWRlZCBieSByaXNrDQpkYXRhLnRoYWkucmlzayA8LSBkZl90aGFpICU+JQ0KICBncm91cF9ieShyaXNrKSAlPiUNCiAgc3VtbWFyaXNlKGNvdW50ID0gbigpKSAlPiUNCiAgbXV0YXRlKHBlcmNlbnQgPSAoY291bnQgLyBzdW0oY291bnQpICogMTAwKSAlPiUgcm91bmQoMikpICU+JQ0KICBmaWx0ZXIocGVyY2VudD40KSAlPiUNCiAgYXJyYW5nZShkZXNjKHBlcmNlbnQpKSANCmRhdGEudGhhaS5yaXNrDQpgYGANCg0KYGBge3J9DQpkYXRhLnRoYWkucmlzayRyaXNrW2RhdGEudGhhaS5yaXNrJHJpc2sgPT0gIkMiXSA8LSAiQ2xvc2UgY29udGFjdCB3aXRoIHRoZSBwYXRpZW50Ig0KZGF0YS50aGFpLnJpc2skcmlza1tkYXRhLnRoYWkucmlzayRyaXNrID09ICJDR3VzdERyIFNhbXV0IFNhR0dBRiJdIDwtICJDbHVzdGVyIFNhbXV0IFNha2hvbiINCmRhdGEudGhhaS5yaXNrJHJpc2tbZGF0YS50aGFpLnJpc2skcmlzayA9PSAiRiJdIDwtICJTdGF0ZSBRdWFyYW50aW5lIg0KZGF0YS50aGFpLnJpc2skcmlza1tkYXRhLnRoYWkucmlzayRyaXNrID09ICJHIl0gPC0gIkdvIHRvIGEgY3Jvd2RlZCBwbGFjZSINCmRhdGEudGhhaS5yaXNrJHJpc2tbZGF0YS50aGFpLnJpc2skcmlzayA9PSAiQSJdIDwtICJQZW9wbGUgdHJhdmVsbGluZyBmcm9tIGFicm9hZCINCmRhdGEudGhhaS5yaXNrJHJpc2tbZGF0YS50aGFpLnJpc2skcmlzayA9PSAiVUZHRkF3RiJdIDwtICJVbmtub3duIg0KZGF0YS50aGFpLnJpc2skcmlza1tkYXRhLnRoYWkucmlzayRyaXNrID09ICJEIl0gPC0gIkNhcmVlciBhdCByaXNrIg0KZGF0YS50aGFpLnJpc2skcmlza1tkYXRhLnRoYWkucmlzayRyaXNrID09ICJDR3VzdERyIFJheUFGZyJdIDwtICJDbHVzdGVyIFJheW9uZyINCmRhdGEudGhhaS5yaXNrJHJpc2tbZGF0YS50aGFpLnJpc2skcmlzayA9PSAiQ0d1c3REciBQYXR0YXlhIENhc0NGQSJdIDwtICJDbHVzdGVyIFBhdHRheWEgQ2FzaW5vIg0KZGF0YS50aGFpLnJpc2skcmlza1tkYXRhLnRoYWkucmlzayRyaXNrID09ICJDR3VzdERyIEFGZyBUR0FGZyJdIDwtICJDbHVzdGVyIEFuZyBUaG9uZyINCmRhdGEudGhhaS5yaXNrDQpgYGANCg0KYGBge3J9DQpkYXRhLnRoYWkucmlzayRyaXNrIDwtIGZhY3RvcihkYXRhLnRoYWkucmlzayRyaXNrLCBsZXZlbHMgPSBhcy5jaGFyYWN0ZXIoZGF0YS50aGFpLnJpc2skcmlzaykpDQpkYXRhLnRoYWkucmlzayRyaXNrDQpnLnRoLnJpc2sgPC0gZGF0YS50aGFpLnJpc2sgJT4lIA0KICBnZ3Bsb3QoYWVzKHggPSAiIiwgeSA9IHBlcmNlbnQsIGZpbGwgPSByaXNrKSkgKw0KICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5Iiwgd2lkdGggPSAwLjUpICsNCiAgY29vcmRfcG9sYXIoInkiKSArDQogIHRoZW1lX3ZvaWQoKSArDQogIGxhYnModGl0bGU9J1Jpc2sgb2YgVGhhaSBDb25maXJtZWQgQ2FzZXMoSmFuIDIwMjAgLSBKYW4gMjAyMSknKSsNCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbCA9IHBhc3RlMChwZXJjZW50LCAiJSIpKSwgY29sb3IgPSAiQmxhY2siLCBzaXplID0gMywgcG9zaXRpb24gPSBwb3NpdGlvbl9zdGFjayh2anVzdCA9IDAuNSkpICsNCiAgZ3VpZGVzKGZpbGwgPSBndWlkZV9sZWdlbmQocmV2ZXJzZSA9IFRSVUUpKSANCmcudGgucmlzaw0KYGBgDQoNCmBgYHtyfQ0KIyBDb25maXJtZWQgY2FzZXMgZGl2aWRlZCBieSBhZ2UNCmRhdGEudGhhaS5hZ2UgPC0gZGZfdGhhaSAlPiUNCiAgZ3JvdXBfYnkoYWdlLHNleCkgJT4lIA0KICBmaWx0ZXIoYWdlICE9IDApJT4lDQogIHN1bW1hcmlzZShjb3VudCA9IG4oKSkgJT4lDQogIGFycmFuZ2UoZGVzYyhjb3VudCkpDQpkYXRhLnRoYWkuYWdlDQpgYGANCg0KYGBge3J9DQpkYXRhLnRoYWkuYWdlIDwtIGRhdGEudGhhaS5hZ2UjICU+JSBmaWx0ZXIoc2V4ICE9ICJVbmtub3duIikNCmdncGxvdChkYXRhLnRoYWkuYWdlLGFlcyh4PWFnZSx5PWNvdW50LGZpbGw9c2V4KSkrZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIpKw0KICAgICAgbGFicyh0aXRsZT0nQWdlIG9mIFRoYWkgQ29uZmlybWVkIENhc2VzIChTdGFjayknKStndWlkZXMoZmlsbD1ndWlkZV9sZWdlbmQocmV2ZXJzZSA9IFQpKQ0KYGBgDQoNCmBgYHtyfQ0KIyBDb25maXJtZWQgY2FzZXMgZGl2aWRlZCBieSBuYXRpb25hbGl0eQ0KZGF0YS50aGFpLm5hdGlvbmFsaXR5IDwtIGRmX3RoYWkgJT4lDQogIGdyb3VwX2J5KG5hdGlvbmFsaXR5KSAlPiUNCiAgc3VtbWFyaXNlKGNvdW50ID0gbigpKSAlPiUNCiAgZmlsdGVyKGNvdW50ID4gMTEpJT4lDQogIGFycmFuZ2UoZGVzYyhjb3VudCkpDQpkYXRhLnRoYWkubmF0aW9uYWxpdHkkbmF0aW9uYWxpdHlbZGF0YS50aGFpLm5hdGlvbmFsaXR5JG5hdGlvbmFsaXR5ID09ICI/Pz8/Pz8/PyJdIDwtICJVbmtub3duIg0KZGF0YS50aGFpLm5hdGlvbmFsaXR5JG5hdGlvbmFsaXR5W2RhdGEudGhhaS5uYXRpb25hbGl0eSRuYXRpb25hbGl0eSA9PSAiIl0gPC0gIlVua25vd24iDQpkYXRhLnRoYWkubmF0aW9uYWxpdHkkbmF0aW9uYWxpdHlbZGF0YS50aGFpLm5hdGlvbmFsaXR5JG5hdGlvbmFsaXR5ID09ICJCdXJtYSJdIDwtICJNeWFubWFyIg0KZGF0YS50aGFpLm5hdGlvbmFsaXR5DQpgYGANCg0KYGBge3J9DQpnZ3Bsb3QoZGF0YS50aGFpLm5hdGlvbmFsaXR5KSsgDQogIGdlb21fYmFyKGFlcyh4PW5hdGlvbmFsaXR5LHk9Y291bnQsZmlsbCA9IGNvdW50KSxzdGF0ID0gImlkZW50aXR5IixzaG93LmxlZ2VuZCA9IEYpKw0KICAgIGdndGl0bGUoJ05hdGlvbmFsaXR5IG9mIFRoYWkgQ29uZmlybWVkIENhc2VzIChsb2cgc2NhbGUpJykrDQogICAgdGhlbWUocGxvdC50aXRsZSA9IGVsZW1lbnRfdGV4dChoanVzdCA9IDAuNSkpKw0KICAgIHRoZW1lKGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTApLA0KICAgICAgICAgICAgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTApLA0KICAgICAgICAgICAgcGxvdC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTApKSsNCiAgICB4bGFiKCIiKSt5bGFiKCJOdW1iZXIgb2YgTmF0aW9uYWxpdHkiKSsNCiAgIHNjYWxlX3lfY29udGludW91cyh0cmFucz0nbG9nMTAnKStjb29yZF9mbGlwKCkrc2NhbGVfZmlsbF9ncmFkaWVudChsb3c9ImJsdWUiLGhpZ2ggPSAicmVkIikNCmBgYA0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0K